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Abstract: Correlation functions in Liouville theory are meromorphic functions of the 
Liouville momenta, as is shown explicitly by the DOZZ formula for the three-point 
function on S 2 . In a certain physical region, where a real classical solution exists, the 
semiclassical limit of the DOZZ formula is known to agree with what one would expect 
from the action of the classical solution. In this paper, we ask what happens outside 
of this physical region. Perhaps surprisingly we find that, while in some range of the 
Liouville momenta the semiclassical limit is associated to complex saddle points, in 
general Liouville's equations do not have enough complex-valued solutions to account 
for the semiclassical behavior. For a full picture, we either must include "solutions" 
of Liouville's equations in which the Liouville field is multivalued (as well as being 
complex- valued) , or else we can reformulate Liouville theory as a Chern-Simons theory 
in three dimensions, in which the requisite solutions exist in a more conventional sense. 
We also study the case of "timelike" Liouville theory, where we show that a proposal 
of Al. B. Zamolodchikov for the exact three-point function on S 2 can be computed by 
the original Liouville path integral evaluated on a new integration cycle. 
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1 Introduction 

Quantum Liouville theory has been studied extensively since it was first introduced 
by Polyakov several decades ago in the context of non-critical string theory [1]. Since 
then it has been invoked as a model for higher-dimensional Euclidean gravity, as a non- 
compact conformal field theory, and as a dilaton background in string theory. Among 
more recent developments, Liouville theory has been found [2] to have a connection to 
four-dimensional gauge theories with extended supersymmetry and has emerged as an 
important component of speculative holographic duals of de Sitter space and the multi- 
verse [3-5]. In many of these applications the Liouville objects of interest are evaluated 
at complex values of their parameters. The goal of this paper is to understand to what 
extent these analytically continued objects are computed by appropriately continued 
versions of the Liouville path integral. 

Liouville theory has been studied from many points of view, but the essential point 
for studying the question of analytic continuation is that several nontrivial quantities 
are, remarkably, exactly computable. A basic case is the expectation value on a two- 
sphere S 2 of a product of three primary fields of Liouville momentum oij, i — 1,2, 3: 



For this correlation function, there is an exact formula - known as the DOZZ formula 
[6, 7]. The existence of such exact formulas makes it possible to probe questions that 
might otherwise be out of reach. We will exploit this in studying the semiclassical limit 
of Liouville theory in the present paper. 

Liouville theory is conveniently parametrized in terms of a coupling constant b, 
with the central charge being c = 1 + QQ 2 , where Q = b + b~ x . For a semiclassical 
limit, we take b — > 0, giving two interesting choices for the Liouville momenta. We 
can consider "heavy" operators, ojj = rji/b, with rji fixed as b — > 0. The insertion of 
a heavy operator changes the saddle points which dominate the functional integral. 
Alternatively, we can consider "light" operators, «j = boi, where now we keep <7j fixed 
as b — > 0. Light operators do not affect a saddle point; they just give us functions that 
have to be evaluated at a particular saddle point. We will consider both cases in the 
present paper. 
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A real saddle point in the Liouville path integral is simply a real solution of the 
classical equations of motion 

— d a d a (j) + QTZ + 87r6yue 26 * = 0. (1.2) 

Such a solution is a critical point 1 of the classical action S. Path integrals are most 
simple if they are dominated by a real saddle point. For the Liouville correlation 
function of three heavy fields on S 2 , there is a real saddle point that dominates the 
semiclassical limit of the path integral if and only if the rji are real, less than 1/2, and 
obey YliVi > 1- These inequalities, which define what we will call the physical region, 
were described in [8] and the explicit solution was described and its action computed 
in [7]. The action evaluated at the classical solution is of the form S = G(r]i, rj 2 , rj 3 )/b 2 
where the function G can be found in (4.26). In [7], it was shown that, in the physical 
region, the weak coupling limit of the three-point function of heavy fields is 

/j|exp(277^/6)(^)y ~ exp(-S) = exp (—G(r)i,T)2,r) 3 )/b 2 + 0(1)) , (1.3) 

as one would expect. 

1.1 Analytic Continuation 

What happens when we leave the physical region? There is no problem in continuing 
the left hand side of (1.3) beyond the physical region. Indeed, the DOZZ formula shows 
that the left hand side of the three-point function (1.1) is, for fixed b, a meromorphic 
function of the variables c^, and in particular one can analytically continue the rji to 
arbitrary complex values. Similarly, the DOZZ formula is analytic in b 2 for b 2 > 0. 

What happens on the right hand side of eqn. (1.3)? If continued outside the 
physical region, the function G extends to a multivalued function of complex variables 
7]i. This multivaluedness takes a very simple form. There are branch points at integer 
values of the rji or of simple sums and differences such as rji + rj2 — rj 3 . The monodromy 
around these branch points takes the form 

G -> G + 2ni + ^ m^A , (1.4) 

where n and the rrii are integers and the mi are either all even or all odd. 

There actually is one specific region of complex rji - the case that Re rji = 1/2 and 
Im?7j > 0, so that the external states are normalizable states in the sense of [8] - in 

^^When heavy operators are present, they add delta function terms in (1.2) and also make contri- 
butions to the action S. In this introduction, we will omit such details. 
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which a semiclassical interpretation of the DOZZ formula is available [9, 10] in terms 
of real singular solutions of Liouville's equations that have a natural interpretation 
in terms of hyperbolic geometry. The action of these singular solutions is given by a 
particular branch of the multivalued function G, for the values of the r)i in question. 
(The solutions themselves are given by an analytic continuation of those constructed in 
[7] and thus are a special case of the solutions we discuss later.) In the present paper, 
we aim to interpret the DOZZ formula semiclassically for arbitrary complex rji. 

Our investigation started with the question of how to interpret the multivaluedness 
(1.4). The most obvious interpretation is that the branches of G might correspond to 
the actions of complex solutions of the Liouville equation. Outside the physical region, 
the correlation function of a product of heavy fields does not have a real saddle point, 
but one might hope that it would have one or more complex saddle points. 

The most obvious notion of a complex saddle point is simply a complex-valued 
solution of the classical Liouville equation (1.2). Such a solution is a critical point of 
the Liouville action S, interpreted now as a holomorphic function of a complex Liouville 
field cf). Optimistically, one would think that Liouville theory for the case of three heavy 
operators on S 2 has complex-valued solutions parametrized by the integers n and rrii 
that appeared in eqn. (1.4). Then one would hope that for any given values of the rji, 
the path integral could be expressed as a sum of contributions from the complex saddle 
points. Which saddle points must be included (and with what coefficients) would in 
general depend on the rji, as Stokes phenomena may intrude. 

To appreciate the analytic continuation of path integrals, one needs to know that 
to a given saddle point one can associate, in principle, much more than a perturbative 
expansion. The basic machinery of complex saddle points and Stokes phenomena says 
the following. 2 Let S be the set of complex saddle points; these are also known as 
critical points of the complexified action. To each p e S, one associates an integration 
cycle 3 C p in the complexified path integral. Roughly speaking, C p is defined by steepest 
descent, starting at the critical point p and descending by gradient flow with respect 
to the "Morse function" h = —Re S. The C p are well-defined for generic values of the 
parameters; in our case, the parameters are b and the rji. The definition of C p fails if 
two critical points p and p have the same value of Im S and unequal values of Re S. 
In this case, the difference S p — S p is either positive or negative (we write S p for the 
value of S at p and similarly h p for the value of h); if for example it is positive, then 

2 See for example section 2 of [11] for an elementary explanation, much more detailed and precise 
than we can offer here. Some familiarity with this material is necessary for a full appreciation of the 
present paper. 

3 An integration "cycle" is simply the multi-dimensional analog of an integration "contour." For 
simplicity, we assume that the critical points are isolated and nondegenerate. 
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the jumping of integration cycles takes the form 

C p ^C p + nC p , (1.5) 

for some integer n. 

Any integration cycle C on which the path integral converges can always be ex- 
pressed in terms of the C p : 

C = ]Ta%, (1.6) 

pes 

with some coefficients a p . In particular - assuming that the machinery of critical points 
and Stokes walls applies to Liouville theory, which is the hypothesis that we set out 
to test in the present paper - the integration cycle for the Liouville path integral must 
have such an expansion. The subtlety is that the coefficients in this expansion are not 
easy to understand, since one expects them to jump in crossing Stokes walls. However, 
there is one place where the expansion (1.6) takes a simple form. In the physical region, 
one expects Liouville theory to be defined by an integral over the ordinary space of real 
fields. On the other hand, in the physical region, there is a unique critical point po 
corresponding to a real solution. Starting at a real value of <fi and performing gradient 
flow with respect to h, <fi remains real. (When is real, h is just the ordinary real 
Liouville action.) So C po is just the space of real fields as long as the rji are in the 
physical region. In the physical region, the expansion (1.6) collapses therefore to 

C = C P0 . (1.7) 

In principle - if the machinery we are describing does apply to Liouville theory - the 
expansion (1.6) can be understood for any values of the rji and b by starting in the 
physical region and then varying the parameters at will, taking Stokes phenomena into 
account. 

If one knows the coefficients in the expansion (1.6) for some given values of the 
parameters, then to determine the small b asymptotics of the path integral 

Z = J D<j> exp(-S) (1.8) 

is straightforward. One has Z = ^2 p a p Z p , with Z p = j c Dip exp(-S). On the other 
hand, the cycle C p was defined so that along this cycle, h = —Re S is maximal at the 
critical point p. So for small b, 

Z p ~ exp(S p ). (1.9) 

The asymptotic behavior of Z is thus given by the contributions of those critical points 
that maximize h p = Re(—S p ), subject to the condition that a p ^ 0. 
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At this point, we can actually understand more explicitly why (1.7) must hold in 
the physical region. A look back to (1.4) shows that as long as b and the rji are real, 
all critical points have the same value of Re S. So all critical points with a p ^ are 
equally important for small b in the physical region. Thus, the computation of [7] 
showing that in the physical region the Liouville three-point function is dominated by 
the contribution of the real critical point also shows that in this region, all other critical 
points have a p = 0. 

1.1.1 The Minisuperspace Approximation 

The Gamma function gives a practice case for some of these ideas. (For a previous 
analysis along similar lines to what we will explain, see [12]. For previous mathematical 
work, see [13, 14].) The most familiar integral representation of the Gamma function 
is 

poo 

T(z)= / d*t* _1 exp(-t), Rez>0. (1.10) 
Jo 

A change of variables t = converts this to 

/oo 
d<\> exp(-S), (1.11) 
-oo 

where the "action" is 

S=-z<f> + e+. (1.12) 

This integral is sometimes called the minisuperspace approximation [8, 15, 16] to Li- 
ouville theory, as it is the result of a truncation of the Liouville path integral to the 
constant mode of <fi (and a rescaling of <fi to replace e 2b< ^ by e^). 

If z is real and positive, the action S has a unique real critical point at = log z, 
and this is actually the absolute minimum of S (on the real axis). We call this critical 
point pq. Gradient flow from p keeps real, so the corresponding integration cycle C Po 
is simply the real axis. If z is not real but Re z > 0, then C Po , defined by gradient flow 
from p , is not simply the real axis, but is equivalent to it modulo Cauchy's Residue 
theorem. The original integral (1.10) or (1.11) can be approximated for z — > oo in the 
half-space Rez > by an expansion near the critical point p$, at which the value of 
the action is S = —z log z + z. The contribution of this critical point leads to Stirling's 
formula T(z) ~ exp(zlogz — z + 0(logz)), Rez > 0. 

The Gamma function can be analytically continued beyond the half-space to a 
meromorphic function of z, defined in the whole complex plane with poles at non- 
positive integers. This is analogous to the fact that the exact Liouville three-point 
function ( 1 . 1 ) is a meromorphic function of the a» , even when we vary them to a region 
in which the path integral over real <fi does not converge. The analytic continuation of 
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the Gamma function to negative Kez can be exhibited by deforming the integration 
contour in (1.11) into the complex <fr plane as z varies. To understand the behavior 
of the integral for Re z < 0, it helps to express this integral in terms of contributions 
of critical points. The complex critical points of S are easily determined; they are the 
points p n with 

<f) n = log z + 2irin, nGZ. (1-13) 

For Kez > 0, the integration contour defining the Gamma function is simply C po , but 
for negative Re z, the integration contour has a more elaborate expansion J2 n eT ^°™> 
where T is determined in Appendix C. Once one determines T, the analog of Stirling's 
formula for Re z < is immediate. 

The essential point is that the integration contour in the definition of the Gamma 
function can be chosen to vary smoothly as z varies, but its expression as a sum of 
critical point contours C Pn jumps in crossing the Stokes walls at Re z = 0. The present 
paper is an attempt to understand to what extent the machinery just sketched actually 
applies to the full Liouville theory, not just the minisuperspace approximation. 

1.2 Complex Solutions Of Liouville Theory 

The result of our investigation has not been as simple as we originally hoped. The 
classical Liouville equations do not have enough complex critical points to account for 
the multivaluedness (1.4), at least not in a straightforward sense. 

As soon as one allows the Liouville field <fi to be complex, one meets the fact that 
the classical Liouville equations are invariant under 

(p + ikw/b, keZ. (1.14) 

This assertion, which extends what we just described in the minisuperspace approxi- 
mation, actually accounts for part of the multivaluedness (1.4). The shift (1.14) gives 

G^G + 27rik(l-J2 iV i). 

This is all we get by considering complex solutions of the Liouville equations in a 
simple way. For example, in the physical region, even if we allow to be complex, the 
most general solution of Liouville's equations is the one described in [7], modulo a shift 
(1.14). We prove this in section 3 by adapting standard arguments about Liouville's 
equations in a simple way. 

Outside of the physical region, the solutions of the complex Liouville equations are 
no more numerous. One can try to find some complex solutions by directly generalizing 
the formulas of [7] to complex parameters. For a certain (difficult to characterize) range 
of the parameters rji and b, this procedure works and gives, again, the unique solution 
of the complex Liouville equations, modulo a transformation (1.14). In other ranges of 
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the parameters, the formulas of [7] do not generalize and one can then argue that the 
complex Liouville equations have no solutions at all. 

The way that the formulas of [7] fail to generalize is instructive. 4 In general, when 
one extends these formulas to complex values of rji and b, branch points appear in the 
solution and (ft is not singlevalued. (The quantities such as exp(260) that appear in 
the classical Liouville equations remain singlevalued.) Singlevaluedness of cf) places a 
serious constraint on the range of parameters for which a complex critical point exists. 

We will show that, after taking into account the symmetry (1.14), ordinary, single- 
valued complex solutions of Liouville's equations suffice for understanding the semiclas- 
sical asymptotics of the Liouville two-point function, and also for understanding the 
semiclassical asymptotics of the three-point function in a somewhat larger region than 
considered in [7]. In particular we will see that the old "fixed-area" prescription for 
computing correlators outside the physical region can be replaced by the machinery of 
complex saddlepoints, which makes the previously-subtle question of locality manifest. 
But for general values of the rji, there are not enough singlevalued complex solutions 
to account for the asymptotics of the three-point function. 

What then are we to make of the semiclassical limit outside of the region where 
solutions exist? Rather surprisingly, we have found that allowing ourselves to use the 
multivalued "solutions" just mentioned in the semiclassical expansion enables us to 
account for the asymptotics of the DOZZ formula throughout the full analytic contin- 
uation in the rji. There is a simple prescription for how to evaluate the action of these 
"solutions", and which has them as stationary points. This prescription agrees with 
the conventional Liouville action on singlevalued solutions and produces its analytic 
continuation when evaluated on the multivalued "solutions". 5 In particular once is 
multivalued, to evaluate the action one must pick a branch of <fi at each insertion point 
of a heavy vertex operator exp(2r]j0/6); allowing all possible choices, one does indeed 
recover the multivaluedness expected in (1.4). We do not have a clear rationale for why 
this is allowed. For one thing if we do not insist on expanding on cycles attached to 
critical points as in (1.6) then it seems clear that for any value of rji we can always find 
an integration cycle that passes only through single-valued field configurations simply 
by arbitrarily deforming the original cycle in the physical region in a manner that pre- 
serves the convergence as we continue in rji. It is only when we try to deform this cycle 
in such a way that the semiclassical expansion is transparent that we apparently need 

This was anticipated by A. B. Zamolodchikov, whom we thank for discussions. 
5 As explained in section 5.1 the prescription is essentially to re-express the Liouville field in terms 
of the "physical metric" gtj = e 2b ^gij, which is always single valued. The branch points of (f> become 
isolated divergences of the metric that have a very specific form, and which turn out to be integrable 
if a "principal value" regularization is used. 
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to consider these exotic integration cycles attached to multivalued "solutions" . 

We also attempt to probe the classical solutions that contribute to the three-point 
function of heavy primary fields by inserting a fourth light primary field. This is not 
expected to significantly modify the critical points contributing to the path integral, 
but should enable one to "measure" or observe those critical points. If the light primary 
field is "degenerate" in the sense of [17], then one can obtain a very concrete formula 
for the four-point function, and this formula supports the idea that the three-point 
function is dominated by the multivalued classical solution. When the light operator 
is non-degenerate the situation is more subtle, a naive use of the multivalued solution 
suggests an unusual singularity in the Liouville four-point function which we are able 
to prove does not exist. We speculate as to the source of this discrepancy, but we have 
been unable to give a clear picture of how it is resolved. 

1.3 Liouville Theory And Chern-Simons Theory 

Since it somewhat strains the credulity to believe that the Liouville path integral 
should be expanded around a multivalued classical solution, we have also looked for 
another interpretation. Virasoro conformal blocks in two dimensions have a relation to 
Chern-Simons theory in three-dimensions with gauge group 5X(2,C) that was identi- 
fied long ago [18]. An aspect of this relation is that quantization of Teichmuller space 
[19], which is an ingredient in SL(2, C) Chern-Simons theory, can be used to describe 
Virasoro conformal blocks [20]. Since Liouville theory can be constructed by gluing 
together Virasoro conformal blocks for left- and right-movers, it should also have an 
expression in terms of Chern-Simons theory; one hopes to express Liouville theory on 
a Riemann surface £ in terms of Chern-Simons on E x /, where J is a unit interval. 
The boundary conditions required at the ends of / are those of [18]. These bound- 
ary conditions have recently been reconsidered and the relation between Liouville and 
Chern-Simons theory developed in more detail [21]. 

Given these facts, instead of looking for complex solutions of Liouville theory on E, 
we can look for complex solutions of SL(2, C) Chern-Simons theory on £ x / with the 
appropriate boundary conditions. Here we find a simpler story than was summarized in 
section 1.2. Solutions are precisely parametrized by the integers n and rrii of eqn. (1.4) 
and the action depends on those parameters in precisely the expected fashion. So a 
possible interpretation of the results of the present paper is that if one wishes to apply 
the machinery of complex saddle points and integration cycles to Liouville theory in a 
conventional way, one should use the Chern-Simons description. Possibly this reflects 
the fact that the gradient flow equations of complex Chern-Simons theory are elliptic 
(as analyzed in [11]); this is not so for complexified Liouville theory. 
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1.4 Timelike Liouville Theory 

As an application of these ideas, we will consider the case of what we will call 
timelike Liouville theory, or Liouville theory with large negative central charge. This 
is the case that b is small and imaginary, so that b 2 < 0. If b is imaginary, then the 
exponential term exp(260) in the Liouville action is of course no longer real. One can 
compensate for this by taking <ft — > i<f), but then the kinetic energy of the Liouville field 
becomes negative, and the Liouville field becomes timelike in the sense of string theory. 
From that point of view, ordinary Liouville theory, in which the kinetic energy of has 
the usual sign, might be called spacelike Liouville theory. We will use that terminology 
occasionally. Timelike Liouville theory has possible applications in quantum cosmology 
[3-5] , and also as the worldsheet description of closed string tachyon condensation [22] . 

It was shown in [23] that the DOZZ formula, when analytically continued in b, 
has a natural boundary of holomorphy on the imaginary axis. On the other hand, it 
was also shown that the Ward identities that lead to the DOZZ formula have a second 
solution - which we will call the timelike DOZZ formula - that is well-behaved on the 
imaginary b axis, but runs into trouble if analytically continued to real b. If b is neither 
real nor imaginary, the two formulas are both well-behaved but different. The timelike 
DOZZ formula has also been independently derived as a possible "matter" theory to be 
coupled to spacelike Liouville in [24] and further studied in [25, 26]. Its first appearance 
seems to be as equation (4.5) in [27], where it appeared as an intermediate step in a 
proposal for the c = 1 limit of Liouville. 6 

From the perspective of the present paper, with all fields and parameters poten- 
tially continued to complex values, timelike Liouville theory and ordinary or spacelike 
Liouville theory are the same theory, possibly with different integration cycles. We will 

6 In [27] it was argued that the timelike DOZZ formula should be multiplied by various nonanalytic 
factors in order for it to describe timelike Liouville. This proposal seems to work only when b = ip/q 
with p, q £ Z, and does not allow continuation to generic a. These modifications are allowed because 
for these special values of b the uniqueness argument for the timelike DOZZ formula breaks down. 
Some interesting applications of timelike Liouville seem to require generic values of b and a, for example 
in coupling to a "matter" CFT with generic c > 25, so we are interested in describing a theory that 
works for generic imaginary b. From the point of view of this paper analyticity in b is also more 
natural to consider; since the integrand of the path integral is an analytic function of its parameters 
the integral should be analytic as well. We will see below that evaluating the Liouville path integral 
in the timelike regime does not produce any nonanalytic factors, so they would have to be put in by 
hand. Schomerus's justification of the extra factors involves wanting the three point function to reduce 
to the two point function, which may be an appropriate requirement in a theory that works precisely 
at b = ip/q (and may be related to Virasoro minimal models). The timelikc DOZZ formula does not 
have this property, but we suggest an alternative interpretation in section 7.4 that does not require 
the new factors. 
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investigate this question and show that the timelike DOZZ formula can indeed come 
from the same path integral that gives the ordinary or spacelike DOZZ formula, with 
an extra factor that represents a change in the integration cycle. 

It was shown in [22, 23, 27] that timelike Liouville theory does not at first seem 
to have all of the usual properties of a conformal field theory; this issue was discussed 
further in [28] but not resolved. The simplest problem in interpreting the timelike DOZZ 
formula in terms of conformal field theory is that naively it appears that the two-point 
function is not diagonal in the conformal dimensions. Our path integral interpretation 
of the timelike Liouville correlators sheds some light on this question; we will argue 
that the two-point function is indeed diagonal and conjecture that the problems which 
have been identified have to do with the existence of new degenerate fields that do 
not decouple in the conventional way. This is possible because of the intrinsically 
nonunitary nature of timelike Liouville. We have not been able to answer the more 
subtle question of which states to factorize correlators on. For a minisuperspace analysis 
of this problem, see [29]. 

1.5 Outline 

An outline of this paper is as follows. In section 2, we review Liouville theory. In 
section 3, we study complex solutions of Liouville's equations on the sphere with heavy 
operators. In section 4, show that the analytic continuation of the DOZZ formula in a 
restricted region can be interpreted in terms of complex classical solutions. In section 
5, we study the full analytic continuation and confront the issue of the nonexistence of 
nonsingular solutions. We then use a fourth light primary field to probe the classical 
configurations contributing to the three-point function of heavy primaries, confirming 
our explanation of the DOZZ analytic continuation in terms of singular "solutions". In 
section 6, we reinterpret the question of complex classical solutions in terms of Chern- 
Simons theory. In section 7, we consider timelike Liouville theory. In section 8 we give 
a brief summary of our results and suggest directions for future work. Finally, in a 
series of appendices, we describe a variety of useful technical results. 

The length of the paper is partly the result of an attempt to keep it self-contained. 
We have written out fairly detailed accounts of a variety of results that are known but 
are relatively hard to extract from the literature. This is especially so in section 2 and 
in some of the appendices. Casual readers are welcome to skip to the conclusion, which 
contains the highlights in compact form. 
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2 Review Of Liouville Theory 



We begin with an overview of Liouville theory. The goal is to present and motivate 
all the existing results that we will need in following sections; there are no new results 
here. Some relatively modern reviews on Liouville theory are [30, 31]; a much older 
one is [8]. Our conventions are mostly those of [7]. 

2.1 Action, Boundary Condition, and Equation of Motion 

The Liouville action, obtained for example by gauge fixing a generic conformal field 
theory coupled to two-dimensional gravity [32], is 

Sl = hJ [ d ^ d ^9 ab + W + 4vr/ie 2 ^] . (2.1) 

Here Q — b + |, and the exponential operator is defined in a renormalization scheme 
using g to measure distances. The metric g is referred to as the "reference" metric (1Z 
is its scalar curvature), while the quantity g a b = e< ■ g a b is referred to as the "physical" 
metric. Since we are viewing Liouville theory as a complete theory in and of itself, the 
"physical" metric is no more physical than the reference one, but it is extremely useful 
for semiclassical intuition so we will often discuss it in what follows. This theory is 
invariant, except for a c-number anomaly, under conformal transformations: 



z = w(z) 
{z' } z') = <p(z,z) - ^log 



dw 



dz 



2 



(2.2) 



Here we use a complex coordinate z — ^ + it; 2 , and w(z) is any locally holomorphic 
function. Under these transformations the renormalized exponential operators have 
conformal weights 

A(e 2 ^) = A(e 2 ^) = a(Q - a), (2.3) 
as we explain in section 2. 2. 7 The stress tensor is 

T(z) = -(d<P) 2 + Qd 2 <P, (2.4) 

and the central charge of the conformal algebra is 

c L = 1 + 6Q 2 = 1 + 6(6 + b- 1 ) 2 . (2.5) 



7 In the terminology that we adopt, the scaling dimension of an operator is A + A, which is twice 
the weight for a scalar operator. 
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We will study this theory on a two-sphere. It is convenient to take the reference 
metric to be the flat metric ds 2 = dz dz, with 

<p = -2Qlogr + 0(1) as r -» oo, r = \z\, (2.6) 

which ensures that the physical metric is a smooth metric on S 2 . This ensures that <fi 
is nonsingular at infinity with respect to (2.2). The intuition for the condition (2.6) is 
that there is an operator insertion at infinity representing the curvature of S 2 , which 
has been suppressed in taking the reference metric to be flat. 

Though the use of a flat reference metric is convenient, with this choice there is 
some subtlety in computing the action; one must regulate the region of integration and 
introduce boundary terms. Following [7], we let D be a disk of radius R, and define 
the action as the large R limit of 

S L = ^~ [ d\ [d a <pd a <P + 47r^e 2 ^] +Q<( (j)d6 + 2Q 2 log R. (2.7) 
4?r Jd 7r J dD 

The last two terms ensure finiteness of the action and also invariance under (2. 2). 8 

The semiclassical limit b — > is conveniently studied with a rescaled field C = 2b(f>, 
in terms of which the action becomes 

h2Sh= lib / d ^[ d ^cd a 4>c + l67ifib 2 e^]+^ ^4> c d9 + 2\ogR + 0(b 2 ), (2.10) 
and the boundary condition becomes 

<f) c {z,z) = -2\og{zz) + 0{1) as \z\ -> oo. (2.11) 
The equation of motion following from this action is 

dd4> c = 2tt//&V c . (2.12) 

If we now define A = irfib 2 to be fixed as b — > 0, then C will have a fixed limit for 
b — > 0. 9 Since g a b = e^ c 5 a b, the physical metric has a good limit as well. The equation 

8 Onc way to interpret them is to note that if we begin with the original Liouville action (2.1) with 
round reference metric 

d S 2 = ^-A_(dr 2 + rW), (2.8) 

then the field redefinition 

^0-Qlog^^y (2.9) 

produces exactly the action (2.7) up to a finite field-independent constant. Rather than trying to keep 
track of this, we will just take the action (2.7) as our starting point. 

intuitively this choice of scaling ensures that the radius of curvature A -1 / 2 of the physical metric 
is large in units of the "microscopic" scale 
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of motion is equivalent to the condition of constant negative curvature of g a b, and this 
is the source of the classical relationship between Liouville's equation (2.12) and the 
uniformization of Riemann surfaces. 

2.2 Conformal Primary Operators and Semiclassical Correlators 

Because of the unusual nature of the transformation (2.2), we can guess that it will 
be exponentials of that transform with definite conformal weights. Classically we see 
that 

= ff^ ff^j ^ gactf-rt (2.13) 

so that classically V a = e 2a ^ is a primary conformal operator with conformal weights 
A = A = aQ [17]. a is called the Liouville momentum. Quantum mechanically, the 
conformal weights of these operators are modified. In free field theory, normal ordering 
contributes a well-known additional term —a 2 to each weight. In Liouville theory, 
the quantum correction is exactly the same, since we can compute the weight of the 
operator V a by considering correlations in a state of our choice. We simply consider 
correlations in a state in which <fi « 0, thus turning off the Liouville interactions 
and reducing the computation of operator scaling to the free field case. 10 So V a has 
conformal weight a(Q — a), as in (2.3). 

In this subsection we will discuss the properties of these operators and their correla- 
tors in more detail at the semiclassical level, in particular seeing how this factor emerges 
in the formula for A. In the following subsections we will review the exact construction 
of Liouville theory that confirms this expression for A beyond the semiclassical regime. 

We will now consider correlation functions of primary fields, 

(v ai ( Zl ,zi)---Va n (zn,Zn)^ = J V<f> c e~ 8 - f[ exp ^MplM ) . (2.14) 

We would like to approximate this path integral using the method of steepest descent 
for small b, but to do so we must decide how the aj's scale with b. The action (2.10) 
scales like b~ 2 , so for an operator to have a nontrivial effect on the saddle points we 
must choose its Liouville momentum a to scale like b~ l . Thus if we want an operator 
to affect the saddle point, we take a = rj/b and keep rj fixed for b — > 0. This gives 
what is conventionally called a "heavy" Liouville primary field. Asymptotically such 
a field has A = 77(1 — i])/b 2 for b — > 0. One can also define "light" operators with 
a = ba, where a is kept fixed for b — > 0. Light operators have fixed dimensions in the 

10 We will see below that this argument requires Re a < Q/2, since otherwise the backreaction of 
the operator will prevent <f> « near the operator. 
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semiclassical limit. Insertion of such an operator has no effect on the saddle point C , 
and to lowest order in b can be approximated by a ^-independent factor of e ai ^ c( - z ' z \ 

Semiclassically the insertion of a heavy operator has the effect of adding an addi- 
tional delta function term to the action, leading to a new equation of motion: 

dd(j) c = 2tt/i6V c - 2tt ^ 2 (£ - &) (2-15) 

i 

Let us assume that in the vicinity of one of the operator insertions we may ignore the 
exponential term. This equation then becomes Poisson's equation: 11 

V 2 c = -87T7^ 2 (£-6)- (2.16) 

This has the solution 

4> c (z,z) = C - Ar]i log \z - Zi\, (2.17) 
so we find that in a neighborhood of a heavy operator we have 

(f) c (z,z) = -Ar]i\og\z - Zi\ + 0(1) asz^-Zi. (2.18) 

We also find that that the physical metric in this region has the form: 

ds 2 = ^-(dr 2 + r 2 d6 2 ) (2.19) 

We can insert this solution back into the equation of motion to check whether the 
exponential is indeed subleading. We find that this is the case if and only if 

Refa) < l -. (2.20) 

If this inequality is not satisfied, then the interactions affect the behaviour of the field 
arbitarily close to the operator. In [8], this was interpreted as the non-existence of 
local operators with Re(^) > |, and the condition that "good" Liouville operators have 
Re (r/) < | is referred to as the Seiberg bound. The modern interpretation of this result, 
as we will see in the following section, is that both a and Q — a correspond to the same 
quantum operator, with a nontrivial rescaling: 

V Q - a = R{a)V a . (2.21) 

R(a) is referred to as the reflection coefficient, for reasons explained in [8, 16]. Either a 
or Q — a will always obey the Seiberg bound, and we can always choose that one when 

11 Note the the convention that Add = V 2 . 
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studying the semiclassical limit. We will thus focus only on semiclassical solutions for 
which all operators have Re(r/j) < \. 

We will in general be interested in complex values of rji, so the metric (2.19) will 
be complex and thus not admit a simple geometric interpretation. For the next few 
paragraphs, however, we will assume that 77, is real to enable us to develop some useful 
intuition. We first observe that since r/j < |, we can do a simple change of variables to 
find a metric 

ds 2 = dr' 2 + r' 2 d6' 2 , (2.22) 

where the coordinate ranges are r' e [0, 00) and & G [0, (1 — 2f]i)2ir). Thus we can 
interpret the effect of the operator as producing a conical singularity in the physical 
metric, with a conical deficit for < rji < | and a conical surplus for rji < 0. Finding 
real saddle points in the presence of heavy operators with real 77's is thus equivalent 
to finding metrics of constant negative curvature on the sphere punctured by conical 
singularities of various strength. 

An interesting additional constraint comes from the Gauss-Bonnet theorem. The 
integrated curvature on a sphere must be positive to produce a positive Euler character, 
so for a metric of constant negative curvature to exist on a punctured sphere the 
punctures must introduce sufficient positive curvature to cancel the negative curvature 
from the rest of the sphere. By integrating equation (2.15) and using the boundary 
condition (2.11) we find a real solution C can exist only if 

5^7fc>l (2.23) 

i 

This inequality along with the Seiberg bound leads to interesting constraints on 
Liouville momenta. In particular for the case of three heavy operators on S 2 , the 
inequalities together imply < rji < \. Unless we satisfy these inequalities, there is no 
real saddle point for the Liouville path integral, even if the rji are all real. The Gauss- 
Bonnet constraint (2.23) also implies that there is no real saddle point for a product of 
light fields on S 2 ; this case amounts to setting all 77, to zero. In particular, there is no 
real saddle point for the Liouville partition function on S 2 . This has traditionally been 
dealt with by fixing the area (calculated in the physical metric) and then attempting 
to integrate over the area; the fixed area path integral has a real saddle point. We will 
develop an alternative based on complex saddle points. 

More generally, if the 77's are complex, then as we mentioned above a saddle point 
cf) c will in general be complex and there is no reason to impose (2.23). 

So far we have not encountered the renormalization issues mentioned at the begin- 
ning of the section. But if we try to evaluate the action (2.10) on a solution obeying 
(2.18), then we find that both the kinetic term and the source term contributed by 
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the heavy operator are divergent. 12 To handle this, again following [7], we perform 
the action integral only over the part of the disk D that excludes a disk di of radius e 
about each of the heavy operators. We then introduce "semiclassically renormalized" 
operators 



Vi 



2n 



Vm(zi,Zi) « e'W exp (t <p c d6 . (2.24) 



ddi 



It is easy to check that this operator multiplied by the usual integrand of the path 
integral (the exponential of minus the action) has a finite limit as e — > when evaluated 

on a solution obeying (2.18). The prefactor e f> 2 in (2.24) contributes a term —2r]f/b 2 
to the scaling dimension of the operator V^/j; this is a contribution of —f] 2 /b 2 to both 
A and A, consistent with the quantum shift —a 2 of the operator weights. We can thus 
incorporate the effects of all the heavy operators by introducing a modified action: 

b 2 S L =-i- / d 2 i (d a <j> c d a (j) c + 16Ae* c ) + ^-(f <f) c d6 + 2 log R 

1D7T JD-Uidi Zn JdD 

" E f^T / tcMi + 2 V 2 loge) (2.25) 



2tt 

The equations of motion for this action automatically include both Liouville's equation 
(2.12) and the boundary conditions (2.11) and (2.18). The final semiclassical expression 
for the expectation value of a product of heavy and light primary fields is 

m 



Vnx(z!,zx) ■ ■ ■ VW(z n ,z n )H ffl (xi,xi) • • • V bam (x mi x m )\ w e Sl[< ^ ] JJe 



i=l 



(2.26) 



Here there are n heavy operators and m light operators, and cj)^ is the solution of (2.15) 
obeying the correct boundary conditions. In this formula effects that are 0(b°) in the 
exponent have been kept only if they depend on the positions or conformal weights of 
the light operators. We will do light operator computations in sections (4.3, 5.4, 5.5, 
7.2.3), and we will be more careful about these corrections there. If there is more than 
one solution, and we will find that in general there will be, then the right hand side of 
(2.26) will include a sum (or integral) over these saddlepoints. 

2.3 DOZZ Formula 

In two-dimensional conformal field theory, the expectation value of a product of 
three primary operators on S 2 is determined up to a constant by conformal symmetry 

12 The exponential term is finite since we are assuming (2.20). 
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[17]. We saw in the previous section that the operators V a are primaries of weight 
A = a(Q — a), so their three-point function must be of the form 



(V ai (Z! , Z\ ) V a2 (z 2 , z 2 ) V a3 (z 3 , z 3 ) ) 



C(aci,a 2 ,a 3 ) 



Here Zij 



j ■ 



U l2 |2(A 1 +A 2 -A 3 )|^ l3 |2(A 1 +A3-A 2 )| Z23 |2(A 2 +A3-A 1 ) ' 

(2.27) 

The function C{a,x, 02,03) is the main dynamical data of any 
CFT. In a CFT with only finitely many primaries, matrix elements of C are often 
called structure constants, but this terminology does not seem entirely felicitous when 
C depends on continuous variables. The DOZZ formula is an analytic expression for C 
in Liouville theory [6, 7]. This proposal satisfies all the expected conditions in Liouville 
theory, and is the unique function that does so; in particular, it is the unique solution 
of recursion relations that were derived in [30, 33]. Knowing (7(0:1,02,03), along with 
rules for a sewing construction of higher order amplitudes can be viewed as an exact 
construction of the quantum Liouville theory. 
The DOZZ formula reads: 

r , 2 l (Q~T. a i)/ b 

C(ai, a 2 , a 3 ) = Tx^(b 2 )b 2 ^ 2 



ToT 6 (2a 1 )T b (2a2)T b (2a 3 ) 



T b (ai + a 2 + a 3 - Q)Y b (cti + a 2 - a 3 )T b (a 2 + a 3 - «i)T fe (ai + a 3 - a 2 ) 
Here T^(x) is an entire function of x defined (for real and positive b) by 



(2.28) 



logT 6 (x) 



00 dt 

T 



(Q/2 



x?e- 1 



sinh 2 ((Q/2-x 



sinh ^ sinh ^ 



< Re(x) < Q. 



(2.29) 

Though this integral representation is limited to the strip < Re(x) < Q, T b (x) has an 
analytic continuation to an entire function of x. This follows from recursion relations 
that are explained in Appendix A, along with other properties of T b . Tq is defined 
as |x=o, and 7(0;) = r(i— L) • ■'■ n ^ e f°ll° w i n g section we will discuss some of the 
motivation for this formula, but for the moment we will just make three observations: 

(1) This expression obeys C(Q — ati, a 2 , a 3 ) = R(ai)C(ai, a 2l a 3 ) with 

(2a-o)/6 T b (2ai - Q) 



R{a) 



n f i 1 (b 2 )b 2 - 2bI 
[n^(b 2 )f a - Q)/b 



T(2aj 



7(2a/6- 1 - l/b 2 ) 1 (2ba-b 2 )' 



(2.30) 



this justifies the reflection formula (2.21). To derive this, one uses the reflection 
symmetry T b (Q — x) = T b (x) and also the recursion relations for T b . 
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(2) The entire expression (2.28) is almost invariant under b — > r, and it becomes so 
if we also send \i — > Jl, with 

7r/l7(l/6 2 ) = [7r/i 7 (& 2 )] (2.31) 

This is a weak-strong duality, in the sense that if /i scales like b~ 2 to produce good 
semiclassical saddle points with finite curvature as b — > 0, then Jib 2 = Ji/b 2 will be 
extremely singular in the same limit so the dual picture will not be semiclassical. 

(3) C(ai, a 2 , a 3 ) as defined in (2.28) is a meromorphic function of the a^, with the 
only poles coming from the zeros of the T^'s in the denominator. In particular it 
is completely well-behaved in regions where the inequalities (2.23) and (2.20) are 
violated. That said, the integral representation of Yj, is only valid in the strip 
< Re(x) < Q, and in the semiclassical limit, for four of the T;,'s in (2.28), the 
boundary of the strip is precisely where the inequality (2.23) or (2.20) breaks 
down. This can lead to a change in the nature of the semiclassical limit. In 
particular when all three a's are real and obey the Seiberg and Gauss-Bonnet 
inequalities, all seven T;,'s can be evaluated by the integral (2.29). This is not an 
accident; in particular, we will argue below that analytically continuing past the 
line Re(r/i + r] 2 + r] 3 ) = 1 corresponds to crossing a Stokes line in the Liouville 
path integral; the number of contributing saddle points increases as we do so. 

2.4 Four-Point Functions and Degenerate Operators 

We will for the most part be studying the semiclassical limit of the DOZZ formula, 
but we will find it extremely helpful to also consider certain four-point functions. 13 In 
two-dimensional CFT, the four-point function on S 2 is the first correlation function 
whose position dependence is not completely determined by conformal symmetry. It 
is strongly constrained, but unfortunately there is much freedom in how to apply the 
constraint and there do not seem to be standard conventions in the literature. We will 
define: 

^ V ai {z x , z x ) V a2 (z 2 , z 2 )V a3 (z 3 , z 3 )V a4 (z 4 , z 4 ) \ 

_ I |2(A 4 -Ai-A 2 -A 3 )| |2(A 2 +A 3 -A 1 -A 4 )| r |-4A 2 i |2(A 1 +A2-A 3 -A 4 ) f y ^\ 

(2.32) 

13 The material discussed here is mostly not required until the final two parts of section 5, so the 
reader who is unfamiliar with the CFT techniques of [17] may wish to stop after equation (2.38) and 
postpone the rest. 
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with the harmonic ratio y defined as: 

y = (2.33) 

^13^24 

This parametrization is chosen so that the limit — > oo, z 3 — y 1, z 2 — > y, and Z\ — > 
is simple: 



lim \z 4 \ 4A ^(V ai (0,0)V a2 (y,y)Va 3 (lA)V a4 (z 4 ^)) = G 12U (y,y) (2.34) 
Using radial quantization as in [17], we can write this as 

Gi234(y,V) = (a^\V a3 (l, l)K 2 (y,y)|a;i). (2.35) 

We can also write C as 

C(ai, «2, a 3 ) = (a 3 |V a2 (l, l)|«i). (2.36) 

In a conventional two-dimensional CFT, these two equations are the starting point for 
the conformal bootstrap program [17]. In this program, one expresses the four-point 
function (2.35) in terms of products of three points functions in two different ways, 
either by inserting a complete set of states between the fields V a3 (l, 1) and V a2 (y,y) 
in (2.35), or by using the operator product expansion to replace the product of those 
two fields with a single field. In Liouville, the situation is more subtle since a is a 
continuous label with complex values and it is not immediately clear what is meant 
by a complete set of states. Similarly, in making the operator product expansion, 
one expands the product of two fields in terms of a complete set of fields, and it is 
again not clear how to formulate this. This problem was solved by Seiberg [8], who 
argued using minisuperspace that the states with a = ® + iP are indeed delta- function 
normalizable for real and positive P, and moreover that these states along with their 
Virasoro descendants are a complete basis of normalizable states. One can check the 
first of these assertions directly from the DOZZ formula by demonstrating that 14 

limC(Q/2 + iPt, e, Q/2 + iP 2 ) = 2ti5{P 1 - P 2 )G(Q/2 + iP x ), (2.37) 

with the two-point normalization G(a) given by 

G(a) = ^ = 1 [n^(b 2 )] iQ - 2a)/h j(2a/b - 1 - l/6 2 ) 7 (26a - 6 2 ). (2.38) 

14 In showing this, one uses the fact that the numerator of the DOZZ formula has a simple zero for 
e — > 0, while the denominator has a double zero for e — > and Pi — P 2 —> 0. One encounters the 
relation lim e ^ e/((Pi - P 2 ) 2 + e 2 ) = tt6(P 1 - P 2 ). 
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Seiberg also argued semiclassically that the state V a2 (y, y)\aii) with both a's real and 
less than Q/2 is normalizable if and only if on + a 2 > This follows from the Gauss- 
Bonnett constraint. If we assume that ot\ and a 2 are in this range, then we can expand 
this state in terms of the normalizable states \Q/2 + iP, k, k). Here \Q/2 + iP, k, k) is 
a shorthand notation for Vg/2+ip(0, 0)\v ac) and its Virasoro descendants. Similarly if 
a 3 + a 4 > y the state (a^V^l, 1) is also normalizable, and we can evaluate (2.32) 
by inserting a complete set of normalizable states. Using (2.21), (2.36), and (2.37) this 
leads to 

1 r°° dP 

G 12 3*(y,y) = - / —C(a 1 ,a 2 ,Q/2 + iP)C(a 3 ,a i ,Q/2-tP) 

x F 12M {A h A P ,y)7- 12 3 4 (A i , A P ,y). (2.39) 

Here Ap = P 2 + Q 2 /4, and the function J-1234 is the familiar Virasoro conformal block 
[17], expressible as 

^ M (A,. A P , y) = ± ft? ^|±^ °> yK (2.40) 

^ C(a 3 ,a4,g/2 + « J P) 

The sum over k is heuristic; it really represents a sum over the full conformal family 
descended from Vq/2+%p- The power of y for a given term is given by the level of the 
descendant being considered, so for example L_iL_ 2 |Q/2 + iP) contributes at order 
y 3 . i s defined in [17], it appears here in the expansion of V^ 2 |a!i) via 



V a2 {y,y)\ai) 



poo jp 

/ —C(a u a 2 , Q/2 + iP)R{Q/2 + iP) [^(Ap-a^a,) 
Jo 27T 

oo 

E PuPn y k f\Q/2 + ^ k, k). (2.41) 



■DO 

X 

k,k=0 



Both and the conformal block itself are universal building blocks for two-dimensional 
CFT's, and conformal invariance completely determines how they depend on the con- 
formal weights and central charge. 

We can then define the general four-point function away from the specified region of 
oti, a 2 by analytic continuation of (2.39). As observed in [7], this analytic continuation 
changes the form of the sum over states. The reason is that as we continue in the a's, 
the various poles of the C's can cross the contour of integration and begin to contribute 
discrete terms in addition to the integral in (2.39). 

One final tool that will be useful for us is the computation of correlators that 
include degenerate fields. A degenerate field in 2D CFT is a primary operator whose 
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descendants form a short representation of the Virasoro algebra, and this implies that 
correlation functions involving the degenerate field obey a certain differential equation 
[17]. Such short representations of the Virasoro algebra can arise only for certain values 
of the conformal dimension. In Liouville theory the degenerate fields have 



where n and m are nonnegative integers [33]. In particular we see that there are 
both light and heavy degenerate fields. We will be especially interested in the light 
degenerate field V_b/2, so we observe here that the differential equation its correlator 
obeys is 



Here A is the conformal weight of the field V^b/2- An identical equation holds for 
correlators involving V_±, with A now being the weight of V_±. For example, by 

26 2b 

applying this equation to the three-point function (V-b/2V ai V a2 ) and using also the fact 
that it must take the form (2.27), one may show that this three-point function vanishes 
unless a 2 = Oil ± 6/2. (This relation is known as the degenerate fusion rule.) We 
can check that the DOZZ formula indeed vanishes if we set a% = —6/2 and consider 
generic ai, a 2 , but there is in important subtlety in that if we simultaneously set a 2 = 
q;i±6/2 and 03 = —6/2 the value of C(ai, a 2 , 013) is indeterminate. (The numerator and 
denominator both vanish.) The lesson is that correlators with degenerate fields cannot 
always be simply obtained by specializing generic correlators to particular values. 

One can actually obtain a good limit for the four-point function with a degenerate 
operator from the integral expression (2.39) [30]. The evaluation is subtle in that there 
are poles of C(a\, a 2 , Q/2 + iP) that cross the contour as we continue a 2 — > —6/2, and 
in particular there are two separate pairs of poles that merge as a 2 — > — 6/2 into double 
poles at the "allowed" intermediate channels a(P) = oil ± 6/2. If we are careful to 
perform the integral with a 2 = — 6/2 + e and then take e — > 0, we find that the formula 
for the four-point function involving the light degenerate field V-b/2 simplifies into a 
discrete formula of the usual type [33]: 




(2.43) 



G(y, y) = 2_j C± i2C34±^i234(Ai, A±, y)J r i234(A i , A±, y). 



(2.44) 



± 
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Here we have taken a 2 = —b/2, and ± corresponds to the operator V ai ±b/2- The raised 
index ± is defined using the two-point function (2.38), so: 

C± 12 = C(ai ± b/2, at, -b/2)R(a 1 ± b/2) = C{a x , -b/2, Q - a x =p b/2). (2.45) 

As just discussed the value of the structure constant on the right cannot be determined 
unambiguously from the DOZZ formula, but the contour manipulation of the four-point 
function gives 

C(a 1 ,-b/2,Q-a 1 Tb/2) = lim [lim e C{a l , -b/2 + 5,Q - a^b/2 + e - 5)] . (2.46) 

<5— >0 U-»0 J 

This definition agrees with a Coulomb gas computation in free field theory [7]. 15 Ex- 
plicitly, from the DOZZ formula we find 

7T/1 



C 



12 



7(-6 2 )7(2ai6) 7 (2 + b 2 - 2ba^ 
C~ 12 = 1 (2.47) 

As shown in [17], by applying the differential equation (2.43) to (2.44) we can 
actually determine J-1234 in terms of a hypergeometric function. This involves using 
SL(2, C) invariance to transform the partial differential equation (2.43) into an ordinary 
differential equation, which turns out to be hypergeometric. 16 The analysis is standard 
and somewhat lengthy, so we will only present the result: 



JWAi, A±, y) = y a * (1 - yfF(A T , B T , C T , y), (2.48) 



with: 



A± = (a 1 ±6/2)(Q-a 1 =F&/ 2 ) 

1 3b 2 
A = 1 

2 4 

a T = A± - A - Ai 
(3 = A - A - A 3 
A T = tK^ - Q/2) + b(a 3 + a* - 6) - 1/2 
B T = Tb{a x - Q/2) + b(a 3 - a 4 ) + 1/2 
C T = l T b(2a 1 -Q). 



15 That computation is based on the observation that for the ctfj's occuring in this structure constant, 
the power of fi appearing in the correlator is either zero or one. This suggests computing the correlator 
by treating the Liouville potential as a perturbation of free field theory and then computing the 
appropriate perturbative contribution to produce the desired power of \l. 

16 Hypergeometric functions will appear repeatedly in our analysis, so in Appendix B we present a 
self-contained introduction. 
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Using this expression and formula (B.15) from the Appendix, Teschner showed that 
(2.44) will be singlevalued only if the structure constant obeys a recursion relation [33]: 

C(a 3 ,« 4 ,«i + 6/2) _ 7(-6 2 ) 
C(«3, «4, a% — b/2) -nji 

7(2a 1 6) 7 (26« 1 - b 2 )-f(b(a 3 + a A - a x ) - b 2 /2) 



x 



7(6(ai + a 4 - a 3 ) - 6 2 /2)7(6(«i + a 3 - a 4 ) - & 2 / 2 )7( & («i + "s + « 4 ) - 1 - 36 2 /2) 
The reader can check that the DOZZ formula indeed obeys this recursion relation. In 
fact, Teschner ran the logic the other way: by combining this recursion relation with 
a similar one derived from the four-point function with degenerate operator V_j_, he 
showed that the DOZZ formula is the unique structure constant that allows both four- 
point functions to be singlevalued. In this version of the logic, C u is determined by 
the Coulomb gas computation rather than the limit (2.46) of the DOZZ formula. This 
at last justifies the DOZZ formula (2.28). 

3 Complex Solutions of Liouville's Equation 

In this section we will describe the most general complex-valued solutions of Liouville's 
equation on S 2 with two or three heavy operators present. The solutions we will present 
are simple extensions of the real solutions given for real 77's in [7]. We will emphasize the 
new features that emerge once complex 77's are allowed and also establish the uniqueness 
of the solutions. One interesting issue that will appear for the three-point function is 
that for many regions of the parameters 771,772,773, there are no nonsingular solutions 
of Liouville's equation with the desired properties, not even complex-valued ones. We 
will determine the analytic forms of the singularities that appear and comment on their 
genericity. 

3.1 General Form of Complex Solutions 

We will first determine the local form of a solution Liouville's equation with flat 
reference metric: 

dd(j) c = 2Ae<K (3.1) 

We have defined A = Tifib 2 , which we hold fixed for b — > to produce a nontrivial 
semiclassical limit. It will be very convenient to parametrize <p c in terms of 

e Mzrz) = YtA^ (3-2) 



which gives equation of motion 



ddf = j(dfdf-l). (3.3) 
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There is a classic device [34] that allows the transformation of this partial differential 
equation into two ordinary differential equations, using the fact that the stress tensor 
(2.4) is holomorphic. In particular, the holomorphic and antiholomorphic components 
of the stress tensor are proportional to W — —d 2 f/f and W = —d f / f respectively. 
We thus have: 

d 2 f + W(z)f = (3.4) 
d 2 f + W(z)f = (3.5) 

with W and W holomorphic. In these equations, we may treat z and z independently, so 
we must be able to write / locally as a sum of the two linearly independent holomorphic 
solutions of the W(z) equation with coefficients depending only on z: 

f = u(z)u(z) — v(z)v(z) (3.6) 

Plugging this ansatz into the W equation, we see that u and v are anti-holomorphic 
solutions of that equation. Going back to the original Liouville equation, we find: 

(udv — vdu){udv — vdu) = 1. (3-7) 

The first factor is a constant since it is the Wronskian evaluated on two solutions of the 
W(z) equation, and similarly the second factor is constant. Both Wronskian factors 
must be nonzero to satisfy this equation, so u and v are indeed linearly independent, 
and similarly u and v. So each pair gives a basis of the two linearly independent 
holomorphic or antiholomorphic solutions of the appropriate equation. We thus arrive 
at a general form for any complex solution of Liouville's equation, valid locally as long 
as the reference metric is ds 2 = dz <g> dz: 

<* 1 1 , , 

\{u{z)u{z)-v{z)v{z)) 2 ' K ' 



with u and v obeying 
and u and v obeying 



d 2 g + W(z)g = (3.9) 



d g + W(z)g = 0. (3.10) 



The representation in (3.8) is not quite unique; one can make an arbitrary invertible 
linear transformation of the pair ( J , with a compensating linear transformation of 



To specify a particular solution, we need to choose W and W and also a basis for 
the solutions of (3.9), (3.10). These choices are constrained by the boundary conditions, 
in particular (2.11) and (2.18). If this problem is undetermined then there are moduli 
to be integrated over, while if it is overdetermined there is no solution. 

In the following subsections, we we will show what happens explicitly in the special 
cases of two and three heavy operators on the sphere. But we first make some general 
comments valid for any number of such operators. The presence of heavy operators 
requires the solution C to be singular at specific points Zi. In terms of /, we need 

f(z,z) ~ \z - Zi\ 2r)l as z -> (3-H) 

Looking at the form (3.8), there are two possible sources of these singularities. The 
first is that at least one of u, v , u, v is singular. The second is that all four functions are 
nonzero but uu — vv = 0, because of a cancellation between the two terms. Assuming 
that this cancellation happens at a place where none of u, v, u, v are singular, we can 
expand 

/ ~ A(z-z )+B(z-z ) + O(\z-z \ 2 ) asz^zo- (3.12) 

Inserting this into (3.3), we find that AB = 1 and thus A and B are both nonzero. It 
thus cannot produce the desired behavior (3.11). 

We will have more to say about this type of singularity later, but for now we will 
focus on singularities that occur because some of the functions are singular. In order 
to produce the behavior (3.11) from singularities of the individual functions, u and v 
must behave as linear combinations of (z — Zi) Vi and {z — Zi) Vi for z — > z^, with similar 
behavior for u, v. To get this behavior, W and W must have double poles at z = z iy 
with suitably adjusted coefficients. A double pole of W or W in a differential equation 
of the form (3.4) is called a regular singular point. A double pole is the expected 
behavior of the stress tensor at a point with insertion of a primary field. 

Moreover, for the solution to be regular at the point at infinity on S 2 , we need 
(2.11), which translates into 

f(z,z) ~ \z\ 2 as \z\ oo. (3.13) 

To achieve this, the two holomorphic solutions of the differential equation (d 2 + W)g = 
should behave as 1 and z, respectively, near z = oo. Asking for this equation to have a 
solution of the form a^ + ao + a-i-z -1 + . . . with arbitrary a_i and a and no logarithms 
in the expansion implies that W vanishes for z — > oo at least as fast as \jz A . This is 
also the expected behavior of the stress tensor in the presence of finitely many operator 
insertions on IR 2 . Given this behavior, the differential equation again has a regular 
singular point at z = oo. 
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We do not want additional singularities in W or W as they would lack a physical 
interpretation. To be more precise, a pole in W leads to a delta function or derivative 
of a delta function in dW. Liouville's equation implies that dW = 0, and a delta 
function correction to that equation implies the existence of a delta function source 
term in Liouville's equation - that is, an operator insertion of some kind. 

Thus for a finite number of operator insertions, W and W have only finitely many 
poles, all of at most second order. In particular, W and W are rational functions. The 
parameters of these rational functions must be adjusted to achieve the desired behavior 
near operator insertions and at infinity. We now study this problem in the special cases 
of two or three heavy operators. 

3.2 Two-Point Solutions 

Specializing to the case of two operators, W should have two double poles and 
should vanish as 1/ z A for z — > oo; W should be similar. Up to constant multiples, these 
functions are determined by the positions of the poles: 



W{z) 



w(l - w)z\ 2 
(z - z x ) 2 {z - z 2 ) 2 



We have picked a convenient parametrization of the constants. In this case, the ODE's 
can be solved in terms of elementary functions. A particular basis of solutions is 



9i{z) 
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~Zl] 


ri 


[z 




Z2Y- 


miz) 
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-zi) 
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9i(z) 
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Z2) 1 - 


92{z) 


= (z 


-zi. 


) 1_ 




[z 


-Z2) 



(3.15) 

It remains to determine w and w in terms of rji and f] 2 and to write the it's and v's in 
terms of this basis. In doing this we need to make sure that (3.11) is satisfied, and also 
that the product of the Wronskians obeys (3.7). Up to trivial redefinitions, the result 
is that we must have rji = r] 2 = w = w = rj, also having 

u(z)= gi (z) (3.16) 
v(z) = g 2 (z) 
u{z) = Kgi(z) 

h{z) 



viz 



k(1 - 2riy\z 12 \ 
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This leads to the solution 
1 



A 



k\z - z^lz - z 2 \ 2 ~ 2 v - \ {1 _ 2v ]^ zu? \z - Zl \ 2 -^\z - z 2 \^ 



2 ■ 



(3.17) 



The criterion r\\ = r] 2 is expected, since in conformal field theory, the two-point function 
for operators of distinct conformal weights always vanishes, k is an arbitrary complex 
number, but it is slightly constrained if we impose as a final condition that / be 
nonvanishing away from the operator insertions. The denominator in (3.17) can vanish 
only if k lies on a certain real curve i in the complex plane [I is simply the real axis if 
r) is real). Omitting the curve I from the complex k plane, and taking into account the 
fact that the sign of k is irrelevant, we get a moduli space of solutions that has complex 
dimension one and that as a complex manifold is a copy of the upper half-plane H. 
Returning now to the general form (3.17), we will make two comments: 

(i) Suppose that r\ is real. To avoid singularities, we cannot have k be real, but we 
can instead choose it to be purely imaginary. e^ c will then be real but negative 
definite, and <p c will be complex. Nonetheless this situation still has a simple 
geometric interpretation: we can define a new metric —e^ c 5 a b, which is indeed 
a genuine metric on the sphere, and because of the sign change it has positive 
curvature! It has two conical singularities, and for positive rfs it describes the 
intrinsic geometry of an American football. This observation is a special case of a 
general bijection between saddle points of spacelike and timelike Liouville, which 
will be explored later. 

(ii) Eqn. (3.17) gives the most general form of e^ c , but this leaves the possibility of 
adding to C itself an integer multiple of 2iri, as in eqn. (1.14). Thus the moduli 
space of solutions has many components and is isomorphic to H x Z. 

3.3 Three-Point Solutions 

For the case of three heavy operators, the potentials W and W must now be rational 
functions with three double poles. Their behaviour at infinity determines them up to 
quadratic polynomials in the numerator, which we can further restrict by demanding 
the correct singularities of u, v, u, and v at the operator insertions. There will be a 
new challenge, however; while we can easily choose a basis of solutions of (3.9) and 
(3.10) with the desired behavior near any one singular point, it is nontrivial to arrange 
to get the right behavior at all three singular points. 
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Insisting that the residues of the poles in W and W have the correct forms to 
produce (3.11) leads to unique expressions for W and W: 



W{z) 
W(z) 



7/i(l - 771)2:122:13 7/2(1 - r] 2 )z 2 iz 2 3 7/3(1 ~ 773)^31^32 



2: — 



z - z 2 



z - z 3 



1 



7/i(l - 7/i),Zi2 gig 7/2(1 - 7/2)^21 ^23 7/3(1 -7/3)^31^32 



Z — Z\ 



z — z 2 



z- z 3 



(z - zi)(z - z 2 )(z - z 3 ) 

1 

(z- Z!)(z- Z 2 )(z- Z 3 ) 

(3.18) 



With these potentials, the differential equation of interest becomes essentially the hy- 
pergeometric equation, modulo an elementary normalization. So the solutions can be 
expressed in terms of hypergeometric functions, or equivalently, but slightly more ele- 
gantly, in terms of Riemann's P functions. 17 P functions are solutions of a differential 
equation with three regular singularities at specified points, and with no singularity at 
infinity. The equations (3.9) and (3.10) are not quite of this form since they do have 
a regular singular point at infinity, but we can recast them into Riemann's form by 
defining g(z) = (z — z 2 )h(z) and g(z) = (z — z 2 )h(z). One can check that the equations 
obeyed by h and h are special cases of Riemann's equation B.5, with the parameters 18 



a = r/i 

P = -m 
7 = 7/3 



a' 
i 



1 - 7/1 

7/2-1 
1 - 7/3- 



(3.19) 



We now observe that the boundary conditions (3.11) ensure that without loss of gener- 
ality we can choose u, v, u, and v to diagonalize the monodromy about any particular 
singular point, say Z\. Also picking a convenient normalization of these functions, we 
have 



u{z) = (z-z 2 )P r "(x) 
v(z) = (z - z 2 )P 1 ~ Vl (x) 
u(z) = ai {z-zt)P m {x) 
v{z) = a 2 {z-zl)P 1 -^{x) 

17 In Appendix B, we give a self-contained development of the minimum facts we need concerning 
hypergeometric and P-functions. The reader unfamiliar with these functions is encouraged to read 
this appendix now. 

18 The unpleasant asymmetry of the second line follows from the definition of h, but a symmetric 
definition introduces significant complication in the formulas that follow so we will stay with this 
choice. 



(3.20) 
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Here a 1; a 2 are complex numbers to be determined, x = z 23 (z — Z\)/z\${z — z 2 ), and the 
P functions explicitly are related to hypergeometric functions by 

P Vl (x) = x m (l - x) m F(r]i +773 — r] 2l r]i + r] 2 + 7] 3 - l,2rji,x) 
P l ~^{x) = x 1 - 711 ^ - x) l -^F(l - 77! + 772 - 773, 2 - Vl - r)2 - Vs, 2-2r h ,x). (3.21) 

We can determine the product a\a 2 by imposing (3.7); by construction we know that 
udv — vdu and udv — vdu are both constant, so to make sure their product is 1 it 
is enough to demand it in the vicinity of z = z\. This is easy to do using the series 
expansion for the hypergeometric function near x = 0, leading to 



a±a 2 = 





zial 


2 


ki2 
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Z23 


l 2 ( 


1 - 2Vi) 2 



(3.22) 



It is clear from the above formulas that / = uu — vv is singlevalued about z = z\. For 
this to also be true near z 2 , z 3 is a non-trivial constraint, which can be evaluated using 
the connection formulas (B.ll). For example, 

f\z - z 2 \- 2 = ai P m (x)P m (x) - a 2 P 1 -' r,1 {x)P 1 -' ni {x) 



will be singlevalued near z = z 3 , which corresponds to x — 1, only if 

O'XO'ijX, T)3^"rji, 1— ?73 Q>i(L\— rj 1 ^^l— r)i ,1— 773 • (3.23) 

The connection coefficients are given by (B.12), so combining this with (3.22) we 
find 



( Ql )2 _ 1 131 1x11 - 12 ~ '1111x11 X /3 ~ '1111x11 x 12 X /3 ~ ^ (3 24) 





2l3 


2 


7(^71 + »72 - ^3)7(^1 + ^73 — 772)7(771 + ^2 + r] 3 - 1) 


1^12 


21 


^23 


2 7(277i) 2 7(r/ 2 + ??3 - ?7i) 



Thus both a\ and a 2 are determined (up to an irrelevant overall sign) , so the solution 
is completely determined. The reader can check that with the ratio given by (3.23) 
the solution is also singlevalued near z 2 . This is a nontrivial computation using (B.13), 
but it has to work, since the absence of monodromy around Zi, z 3 , and 00 implies that 
there must also be none around z 2 . 

The final form of the solution is thus 



A [ ai P^(x)P^(x) - a 2 P 1 -^(x)P 1 -^(x)f 



(3.25) 



In the end, this is simply the analytic continuation in r]i of the real solution presented 
in [7], but our argument has established its uniqueness. 
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There is still a potential problem with the solution. The coefficients 0,1,0,2 were 
completely determined without any reference to avoiding cancellations between the 
terms in the denominator, and it is not at all clear that the denominator has no zeroes 
for generic 77's. It is difficult to study the existence of such cancellations analytically 
for arbitrary 77's, but we have shown numerically that they indeed happen for generic 
complex 77's. If we assume that such a singularity is present at z — z , then we saw 
above that its analytic form is given by (3.12). 

For real 77's, we can say more. When the 77's are real, the right hand side of (3.24) 
is real so a\ is either real or imaginary. If it is imaginary, then (3.22) shows that ci2 
will also be imaginary and with opposite sign for its imaginary part. Moreover for 
real 77's, P m (x)P m (x) and P 1 ~ Vl (x)P 1 ~ ril (x) are strictly positive. Thus when a\ is 
purely imaginary, both terms in the denominator have the same phase and there can 
be no singularities arising from cancellation. The metric e^ c 5 a b will however be negative 
definite, so this will be a complex saddle point for <ft c . If we start with such r/'s and 
allow them to have small imaginary parts then cancellations do not appear at once, 
but we find numerically that if we allow the imaginary parts to become large enough 
then cancellations in the denominator do occur. 

We can also consider the case that the 77's are real and a\ is also real. 02 will then 
be real and with the same sign as a\, so cancellations are now possible. We learned in 
section 2.2 that real solutions can only occur if certain inequalities (2.20) and (2.23) are 
satisfied. So if the 77's are real but violate the inequalities, the denominator in (3.25) 
definitely vanishes somewhere. On the other hand, if the 77's are real and satisfy the 
inequalities, then a real metric of constant negative curvature corresponding to a real 
solution of Liouville's equations does exist. It can be constructed by gluing together two 
hyperbolic triangles, or in any number of other ways. So in this case, the denominator 
in (3.25) is positive definite away from the operator insertions. 19 This is the region 
studied in [7]. 

We conclude with two remarks about the nature of these singularities near a zero of 
the denominator in the formula for e^ c . We first observe that the singularities naturally 
come in pairs since the denominator of (3.25) is symmetric under exchanging x and x, so 
for example if we choose the Z{ to be real then the solution is symmetric under reflection 
across the real z-axis. We secondly comment on the stability of these singularities: the 
general local expansion (3.12) near a zero involves two complex coefficients A and B. 
When these are of unequal magnitude, the existence of a zero of / is stable under small 
perturbations. This is because one can associate to an isolated zero of the complex 
function / an integer- valued invariant, the winding number. To define it, set / = se 1 ^ 

19 We show this explicitly below in Appendix B.5. 
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where s is a positive function and ip is real. Supposing that / has an isolated zero 
at z = Zo, consider e 1 ^ as a function defined on the circle z = zq + ee l9 , for some 
small positive e and real 6. The winding number is defined as ^ § T d9 dip / d9 , and is 
invariant under small changes in /. (If / is varied so that several zeroes meet, then 
only the sum of the winding numbers is invariant, in general.) In the context of (3.12), 
the winding number is 1 for \A\ > \B\, and —1 for \A\ < \B\, and depends on higher 
terms in the expansion if \A\ = \B\. In the case of a zero of the denominator in (3.25), 
one generically has \A\ ^ \B\ if the 77's are complex, so isolated singularities arising 
by this mechanism are stable against small perturbations. When the 77's are real, the 
behavior near singularities of this type requires more examination. 

4 Analytic Continuation and Stokes Phenomena 

In this section, we use the complex classical solutions constructed in the previous section 
to interpret the analytic continuation first of the two-point function (2.38) and then 
of the three-point function as given by the DOZZ formula (2.28). We will find that 
for the two-point function there is a satisfactory picture in terms of complex saddle 
points, which agrees with and we believe improves on the old fixed-area results in the 
semiclassical approximation. For the three-point function we will find that the situation 
is more subtle; we will be able to "improve" on the fixed-area result here as well, but to 
understand the full analytic continuation we will need to confront the singularities at 
which the denominator of the solution vanishes. For ease of presentation we postpone 
our discussion of those singularities until section 5, and we here focus only on the part 
of the analytic continuation that avoids them. We also include the case of three light 
operators as check at the end of the section. 

4.1 Analytic Continuation of the Two-Point Function 

We saw in section 2.4 that the DOZZ formula implies that the Liouville two-point 
function takes the form 

(V a (z 1 ,z 1 )V a (z 2 ,z 2 )) = 

\z 12 \- MQ - a)2 ^ [^ 7 (6 2 )] (Q " 2a)/b 7 (2a/6 - 1 - l/6 2 ) 7 (26a - b 2 )S(0) (4.1) 

The factor of 5(0) is a shorthand which reflects the continuum normalization of the 
operators with a = j + iP and the fact that we have taken the two fields in (4.1) 
to have the same Liouville momentum. It may seem unphysical to study the analytic 
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continuation of a divergent quantity, but as we will review, the divergence has a simple 
semiclassical origin that is independent of a. 20 

This "exact" result for the two-point function does not come from a real Liouville 
path integral, even if a is real. One can easily show that, for the two-point function on 
S 2 , the path integral over real Liouville fields does not converge. Consider a smooth real 
field configuration that obeys the boundary conditions (2.11) and (2.18). The modified 
action (2.25) will be finite. Now consider adding a large negative real number A0 C to 
(f) c . The kinetic term will be unaffected and the exponential term will become smaller in 
absolute value, but the boundary terms will add an extra term A0 C (1 — 2rj). Recalling 
that we always choose the Seiberg bound to be satisfied, we see that by taking A<p c to 
be large and negative we can thus make the action as negative as we wish. The path 
integral therefore cannot converge as an integral over real c 's [8]. 21 

The original approach to resolve this divergence, proposed in [8], was to restrict the 
path integral only to field configurations obeying J G? 2 £e^ c = A. This clearly avoids the 
divergence. However, if one tries to integrate over A, one would get back the original 
divergence, while on the other hand if one simply keeps A fixed, one would not expect 
to get a local quantum field theory. As an alternative proposal, we claim that (4.1) is 
computed by a local path integral over a complex integration cycle. This is analogous to 
the suggestion [35] of dealing with a somewhat similar divergence in the path integral of 
Einstein gravity by Wick rotating the conformal factor of the metric to complex values. 
To motivate our proposal, we will show that the semiclassical limit of (4.1), with a 
scaling as rj/b, is reproduced by a sum over the complex saddle points with two heavy 
operators that we constructed in section 3.2. We interpret this as suggesting that the 
path integral is evaluated over a cycle that is a sum of cycles attached to complex saddle 
points, as sketched in section 1.1. The requisite sum is an infinite sum, somewhat like 
what one finds for the Gamma function for Kez < 0, as described in Appendix C. 
We will also find that the set of contributing saddle points jumps discontinuously as 
T) crosses the real axis. This again parallels a result for the Gamma function, and we 
interpret it Stokes phenomenon. 

20 Indeed if we were to use the Liouville theory as part of a gravity theory where conformal symmetry 
is gauged, then to compute a two-point function of integrated vertex operators we would partially fix 
the gauge by fixing the positions of the two operators and then divide by the volume of the remaining 
conformal symmetries. This would remove this divergent factor. 

21 This divergence should not be confused with the factor of S(0), which we will see has to do with 
an integral over a noncompact subgroup of SL(2, C). In particular we can make the same argument 
for the three-point function with three real a's and find the same divergence if ai < Q, and since 
the DOZZ formula does not have any <5(0) it is clear that this is a different issue [8]. 
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4.1.1 Evaluation of the Action for Two-Point Solutions 

In computing the action of the two-point solution (3.17), we first need to deal with 
taking the logarithm to get <f> c . The branch cut in the logarithm makes this a nontrivial 
operation. To make the following manipulations simpler, we will here relabel k = ik, 
so the solution becomes 

a = -a^7 \ (4 ' 2) 

We choose k to ensure that the denominator has no zeroes. Since we are imposing the 
Seiberg bound, we have Re(l — 2rj) > 0. There is a sign choice in defining k, so we 
will choose it to have positive real part. In particular note that if r\ is real then we can 
have k be real and positive. Our prescription for taking the logarithm will then be 

4> Cj n{z, z) =in + 2niN — log A — 2 log k 

- 21og (\z - Zl \^\z - z 2 \ 2 ^ + — 1 \z - Zl \ 2 -^\z - z 2 \ 2 ^ 

(4.3) 

The choice of branch for the final logarithm is inessential, in the sense that making a 
different choice would be equivalent to shifting the integer iV in (4.3). We will choose 
the branch such that the final logarithm behaves like —4r] log \z — z± \ + (4rj — 4) log \zi% \ 
near z%. Its value away from z\ is defined by continuity; there is no problem in extending 
this logarithm throughout the z-plane (punctured at z\ and z 2 ). 22 We will have no such 
luck for the three-point function; in that case, zeroes of the logarithm are essential. 

We will see momentarily that to compute the action, we need to know the leading 
behaviour near z\ and z 2 , so we observe that 

4> c ,n{z,z) ->■ -4r?log \z - Zi\ + Q as z -> z h (4.4) 



with 



C x = 2ni [ N + i ) - log A - 2 log k + (4-r? - 4) log \z 12 \ 

1 
2 



C 2 = 2-ni [ N + - ) - log A + 21ogK + 4t] log\z 12 \ +41og(l - 2r]). (4.5) 



22 Because of the boundary conditions (3.11), there cannot be monodromy of this logarithm about 
Zi, z 2 even though its argument vanishes there. 
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To verify 23 that the same integer N appears in both C\ and C 2 , we note that this is 
clear for real 77 and k, since then the final logarithm in (4.3) has no imaginary part; in 
general it then follows by continuity. 

Now to compute the modified action (2.25), we use a very helpful trick from [7]. 
This is to compute dSi/drj when Sl is evaluated on a saddle point. A priori, there 
would be r) dependence both implicitly through the functional form of the saddle point 
and explicitly through the boundary terms in Sl, but the variation of (2.25) with 
respect to C is zero when evaluated on a solution and only the explicit 77-dependence 
matters. We thus have the remarkably simple equation: 

h 2dSL = _ Ci _ C2 

drj 

= -2m(2N + 1) + 2 log A + (4 - 877) log \z 12 \ - 4 log(l - 2rj) (4.6) 

We can thus determine Sl[4> c ,n] up to a constant by integrating this simple function, 
and we can determine the constant by comparing to an explicit evaluation of the action 
when 77 = 0. When 77 is zero, the saddle point (4.3) becomes an SL(2, C) transformation 
of a metric which is just minus the usual round sphere 

4> c = in + 2niN - log A - 21og(l + zz). (4.7) 

For this solution we can evaluate the action (2.25) explicitly, finding b 2 So = 2iri(N + 
|) — log A — 2. Now doing the integral of (4.6) our final result for the action (2.25) with 
nonzero 77 is thus 

b 2 S L =2m(N + 1/2)(1 - 2t/) + (2t/ - 1)A + 4(7/ - r/ 2 ) log \z l2 \ 

+ 2 [(1 - 2t/) log (1 - 2t/) - (1 - 2t/)] . (4.8) 

We can observe immediately that the dependence is consistent with the two-point 
function of a scalar operator of weight (77 — r] 2 )/b 2 . This action is independent of k, so 
when we integrate over it this will produce a divergent factor, which we interpret as 
the factor 5(0) in (4.1). 

Before moving on to the exact expression, we we will observe here that this action 
is multivalued as a function of 77, with a branch point emanating from 77 = 1/2, where 
the original solution (4.2) is not well-defined. Under monodromy around this point, N 
shifts by 2, so all even and likewise all odd values of N are linked by this monodromy. 
Of course, to see the monodromy, we have to consider paths in the 77 plane that violate 
the Seiberg bound Re (77) < |. 

23 We thank X. Dong for a discussion of this point and for suggesting the following line of argument. 
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4.1.2 Comparison with Limit of Exact Two-Point Function 

We now compute the semiclassical asymptotics of (4.1). We can easily find that 



(V«(zi,*i)Va(z2,z 2 )> ~ 5(0)\z 12 \- 4 ^-nW x (i-2r,)/b* 



b 2 



(1-2tj)/6 2 



7 



(2^7 - 1) 
b 2 



(4.9) 

The first three factors obviously match on to the result (4.8) that we found in the previ- 
ous section, but the last two have more subtle semiclassical limits. It is not hard to see 
that the factor involving r y(b 2 ) is asymptotic for small positive b to exp |_ 4 ( 1 ~ 2 ^) lQ g fe 
but to understand the final factor, we need to understand the asymptotics of the F 
function at large complex values of its argument. For real positive arguments, this 
is the well-known Stirling approximation, but for complex arguments, the situation is 
more subtle: 



x-x+0{\ogx) R&{x) > Q 

Re(x) < 0. 



r(aO = { , , , , „ 4.10 

I 1 a;log(-!B)-a;+0(log(-a!)) T> ~( „\ ^ n 



This result can be obtained in a variety of ways; because of the fact (see [8, 15, 16] and 
section 1.1.1) that the integral representation of the Gamma function is a minisuper- 
space approximation to Liouville theory, we present in Appendix C a derivation using 
the machinery of critical points and Stokes lines. Using (4.10), we find 



. ,' (2g-m i 

M £2 J ~ e iw(2 V -l)/b 2 _ e -iir(2 V -l)/b 2 eXP 



( b/ ~ 2 * (log(l — 277) -21og6-l) 



6 2 

(4.11) 

So we can write the semiclassical limit as 

{V a {z u z x )V a {z 2 ,z 2 )) ~ ^OJIzial^-^Af 1 -^ 

X p -^[(l"2»?)log(l-2r,)-(l-2r ? )] 1 , o\ 

e i7r(2 ? ?-l)/6 2 _ e -j7r(2r;-l)/6 2 • l^-^ 

All factors now clearly match (4.8) except for the last. To complete the argument, 
setting y = e m ( 2 v-^)/b 2 ^ we neec i to know that the function l/(y — y^ 1 ) can be expanded 
in two ways: 

_J_ = sr y - {2k+1) = -X> 2fc+1 . ( 4 - 13 ) 
y-y £z U 

One expansion is valid for \y\ > 1 and one for \y\ < 1. So either way, there is a set T 
of integers with 



3 i7r(2?7-l)/& 2 _ e -ra-(2r;-l)/& 2 

NeT 



± ^ e 27ri(V T l/2)(2r,-l)/b 2 ^ (<&.14) 
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T consists of nonnegative integers if Im((2r/ — l)/b 2 > and of nonpositive ones if 
Im ( (2?7 — 1) jb 2 < 0. We have to interpret the line Im ((2r] — l)/b 2 ) = as a Stokes line 
along which the representation of the integration cycle as a sum of cycles associated 
to critical points changes discontinuously. If b is real, the criterion simplifies and only 
depends on the sign of Imr]. The sign in (4.14) has an analog for the Gamma function 
and can be interpreted in terms of the orientations of critical point cycles. 

4.1.3 Relationship to Fixed-Area Results 

We will now briefly discuss how to relate this point of view to the more traditional 
fixed-area technique [8]. For this section we restrict to real a's. We begin by defining 
the fixed-area expectation value for a generic Liouville correlator as 

(V ai ■ ■ ■ V an ) A = (/^)(S^-Q)/6 1 n ., h A V ai ■ ■ ■ V an ). (4.15) 

r ((£< «i - Q)/o) 

Assuming that Re(£ a i: — Q) > 0, an equivalent formula is 

r°° HA 

(V ai ---V an ) = j o —e-» A (V ai ---V an ) A . (4.16) 

With the A dependence of {V ai ■ ■ ■ V an ) A being the simple power of A given on the 
right hand side of (4.15), the A integral in (4.16) can be performed explicitly, leading 
back to (4.15). So far this is just a definition, but comparison of (4.16) to the original 
Liouville path integral suggests an alternate proposal for how to compute the fixed- 
area expectation value: evaluate the Liouville path integral dropping the cosmological 
constant term and explicitly fixing the physical area J d 2 ^e^ A = A. Semiclassically we 
can do this using a Lagrange multiplier, 24 which modifies the equation of motion: 

99<Pa = ^(X> " l)e* A - 2tt $>5 2 (£ - (4-17) 

i 

The point to notice here is that when YliVi < 1> if we define C) jv = in + 2iriN + (f>A and 
A = ~ t ne solutions of this equation are mapped exactly into the complex 

saddle points we have been discussing. One can check explicitly for the semiclassical 
two-point function we just computed that the various factors on the right hand side of 
(4.15) conspire to remove the evidence of the complex saddle points and produce the 
usual fixed-area result [7]: 

(V>(1, 1)^(0,0))^ = 2n6(0)G A ( V /b) » 2 7r5(0)e-^( 1 - 2 ")( lo ^ +log ( 1 -^" 1 ). (4.18) 

24 In eqn. (4.17), we set the Lagrange multiplier to the value at which the equation has a solution. 
To find this value, one integrates over the z-plane, evaluating the integral of the left hand side with 
the help of (2.6). 
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Historically the proposal was to use (4.15) in the other direction, as a way to define 
the Liouville correlator when YliVi < 1j but it was unclear that this would be valid 
beyond the semiclassical approximation. We see now how it emerges naturally from 
the analytic continuation of the Liouville path integral. 

4.2 Analytic Continuation of the Three-Point Function 

We now move on to the three-point function. We will initially focus on two partic- 
ular regions of the parameter space of the variables % i = 1, . . . , 3. In what we will 
call Region I, we require that ^ j Re(?jj) > 1, and that the imaginary parts lm(rji) are 
small enough that the solution (3.25) does not have singularities coming from zeroes of 
the denominator. The inequality Yli^ e (Vi) > 1 is needed to prevent the path integral 
over (f> c from diverging at large negative <p c , as discussed above in the context of the 
two-point function. When the rji are actually real and less than 1/2, we get the physical 
region studied in [7], which is the only range of rji in which Liouville's equation has real 
nonsingular solutions. In this sense the three-point function is a simpler case than the 
two-point function, since in that case no choice of rj allowed a real integration cycle for 
the path integral. 

We will also be interested in the region defined by 



where again the imaginary parts are taken to be small enough that there are no singu- 
larities from zeroes of the denominator. We will refer to this as Region II. Note that if 
the imaginary parts are all zero, we can see from (3.22) and (3.24) that a± and ci2 will 
be purely imaginary in this region and there will be no singularities. The third line of 
(4.19) has not appeared before in our discussion; we call it the triangle inequality. Its 
meaning is not immediately clear. It is automatically satisfied when Y2i R- e (Vi) > 1 an d 
Ke(r]i) < 7j, but when X]jR- e (^) < 1 it becomes a nontrivial additional constraint. 

To get some intuition about this constraint, recall that in Region II with real 77's, 
a 1; a 2 are imaginary. The metric — e^ c 5 a b is thus well defined and has constant positive 
curvature. Since we have taken the rji to be positive, the metric has three conical 
deficits. Such metrics have been studied in both the physics and math literature [36, 37], 
and they can be constructed geometrically in the following way. Suppose that we can 
construct a geodesic triangle on S 2 whose angles are 9i — (1 — 2rji)-n. We can glue 
together two copies of this triangle by sewing the edges together, and since the edges 



< Re(77i) < - 




(4.19) 



< Re(r]i + rjj - rj k ) 



-38- 



are geodesies they have zero extrinsic curvature and the metric will be smooth accross 
the junction. The angular distance around the singular points will be 29 i = (1 — 2^)277, 
and as explained in the discussion of (2.22), this is the expected behavior for a classical 
solution with the insertion of primary fields of Liouville momenta rji/b. So this gives 
a metric of constant positive curvature with the desired three singularities. For this 
construction to work, we need only make sure that a triangle exists with the specified 
angles. First note that because of the positive curvature of S 2 , we must have ^ 9i > ir, 
which gives Yli Vi < 1- We can choose one of vertices of the triangle, say the one labeled 
by T]i, to be the north pole, and then the two legs connected to it must lie in great 
circles passing through both the north and the south pole. If we extend these legs all 
the way down to the south pole, then the area between them is a "diangle," as shown 
in Figure 1. The third leg of the triangle then splits the diangle into two triangles, 
the original one and its complement, labelled A and B respectively in the figure. The 
inequality ^ 6i > rr applied to the complementary triangle then gives r] 2 + rj 3 — i]i > 0, 
so this is the source of the triangle inequality in (4.19). As we approach saturating the 
inequality, the complementary triangle B becomes smaller and smaller and the original 
triangle A degenerates into a diangle. Once the inequality is violated, no metric with 
only the three desired singularities exists. 




Figure 1. Spherical triangles. 



We will now compare the semiclassical actions of the complex saddle points (3.25) in 
these two regions with the semiclassical limit of the DOZZ formula. We will see that in 
Region I only the real saddle point contributes (this is expected for reasons explained in 
section 1.1) while in Region II, similarly to what we found for the two-point function, 
infinitely many contribute. We interpret this change as a Stokes phenomenon; the 
condition Re (771 + r] 2 + r] 3 ) = 1 separating the two regions evidently defines a Stokes 
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wall. In Region II, we will initially assume that all three operators are heavy, but in a 
final subsection we will treat the case that they are light and again find evidence for a 
path integral interpretation of the DOZZ formula. 

4.2.1 Evaluation of the Action for Three-Point Solutions 

To evaluate the action for a saddle point contributing to the three-point function, we 
can again use the trick of differentiating the action with respect to r^. So we need to 
determine the asymptotic behaviour of (3.25) near Z(. We denote as 4> Ct N the solution 
corresponding to (3.25), where the subscript iV labels the possibility of shifting <p c by 
27iiN. We will again have 

Ci tv(z, z) — > — 4r/ log \z — Zi\ + Ci as z — )■ Zi, (4.20) 

and to determine C{ we again need to confront the problem of defining the logarithm 
of /. We will first treat Region I, where we define 

cj) CjN (z,z) = 2niN — log A — 4 log \z — z 2 \ 

- 2 log ( ai P m (x)P m (x) - a 2 P l ~ m {x)P 1 - m (x)) . (4.21) 

The branch in the logarithm is chosen so that using (3.24) and the series expansion of 
P 7 ? 1 , we find 

d = 2rriN- log A - (1 - 2rfr) log 

7(2r/i) 2 7(^2 + m~vi) 

The function <f) c ^(z,z) that we get by continuation away from z\ will be singlevalued 
by the same argument as for the two-point function. To find C 2 and C 3 , we can use the 
connection coefficients B.13, B.12, but it is easier to just permute the indices to find 

C 2 = 2niN- log A - (1 - 27/a) log 



\Zl2 


2 


Zl 3 \ 


2 




223 


2 



1^12 


21 


223 


2 




2l3 


2 



-log fo x 2 7 , x > 4 - 2 3 

7(2r/ 2 )^7(li +13 - 12) 



C 3 = 2niN- log A - (1 - 2t] 3 ) log 



223 


21 


2l3l 


2 




2-12 


2 



_ j 7(l3 + 12 - 1l)7(ll + 13 ~ l2)7(ll + 12 + 13 - 1) / 4 

7(2ls) 2 7(li + 12 - Is) 
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As with the two-point function, we can justify the equality of N in the vicinity of 
different points by observing that we may begin with real 77's obeying YliVi > 1> f° r 
which the argument of the logarithm is real and positive. We then continue to the 
desired value of rj on a path that remains in Region I. As before, by continuity iV 
cannot change. As a check of this claim, we observe that paths within Region I cannot 
activate the branch cuts of the logarithms in these expressions for Cj. Indeed, for any 
set of r^'s which is in Region I, all of the arguments of 7(-) have real part between 
zero and one. 7(-) has no zeros or poles in this strip, so any loop in Region I can be 
contracted to a point without changing the monodromy of the logarithm. Thus there 
is no monodromy. 

To compute the action, we need to integrate 

b 2 f± = -Q, (4.25) 

OTji 

which gives 

b 2 S L = \ J2vi - lj log A + {S 1 + S 2 - S 3 ) log \z 12 \ 2 + (5a + S 3 - 5 2 ) log |^i 3 | 2 

+ (s 2 + s 3 - #i) log \z 23 \ 2 + F(ri l + rj 2 - m) + F (m + m~ m) + F (v2 + m- m) 

+ F{ Vl + V2 + V3-I)- F{2 Vl ) - F(2r] 2 ) - F(2t7 3 ) - F(0) 

+ 2niN(l - ^rji). (4.26) 



Here we have 



»7J 

F(rj) = I log-f(x)dx, (4.27) 



with the contour staying in the strip < Re(x) < 1, and also 

Si = 7]i{\ - 7]i). 

The ^-independent constant was determined in [7] by explicitly evaluating the action 
in the case 'YliVi = 1> with the result b 2 SL = Yli<j ^ViVj^°S \ x i ~ x j\ 2 - The integral 
involved is quite difficult and will not be described here. 25 We can observe immediately 
that the ^-dependence in (4.26) is of the correct form for a conformal three-point 
function. 



25 The condition that J^i Vi = 1 means that the flat SX(2,C) bundle over the thrcc-punctured sphere 
associated to the solution actually has abelian monodromy. This can perhaps be used to evaluate its 
action in the Chcrn-Simons description that we will discuss in section 6. 
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We now evaluate the action for rji's in Region II. The manipulations are similar, 
but we now define the branch so that 



d =2ni N 



1 



log A - (1 - 2?7i)log 



1*12 


21 


Zl3 


2 




2 23 


2 



+ 21og 



log 



7(771 + V2- 773)7(771 + V3- 772)7(771 + V2 + m) 

7(2?7i) 2 7(f?2 + V3 - Vi) 



(4.28) 



We have used j(x — 1) = — ^l^ li?) to make sure that when we take the 77's to be 
real (and in Region II), the only imaginary parts comes from the first term. We can 
again permute to find: 



C 2 =2rri N + 



1 



log A - (1 - 2?7 2 ) log 



1*12 


21 


^23 


2 




Z 13 


2 



+ 2logll-J2vi 



- log 



7(771 + V2- 773)7(7/2 + m- 771)7(771 + m + 773) 



C 3 =2ni N + 



1 



7(2?7 2 ) 2 7(77i + 773 - V2) 



log A - (1 - 277 3 )log 



(4.29) 



\ z 23 


21 


213 


2 




212 


2 



+ 21og ll-^TH 



log 



7(772 + m - ^1)7(771 + 773 - 772)7(771 + V2 + V3) 



7(2^ 3 ) 2 7(?7i + 772 - 773) 
Finally we can again integrate this to find 



(4.30) 



b 2 S, 



(v^Vi ~ 1 j log A + (S 1 + 5 2 - 5 3 ) log |^i 2 | 2 + (*i + S 3 - 8 2 ) log l^] 2 

+ (S 2 + S 3 - Si) log I ^23 1 2 + F(r]i + r? 2 — 773) + Ffa + r] 3 - r] 2 ) + F(r) 2 + r) 3 - 771) 
+ F( Vl + m + m ) - F(2 Vl ) - F(2t7 2 ) - F(2 m ) - F(0) 



+ 2 



+ 27rz(iV + 1/2)(1 - ^2 



(4.31) 



Here we determined the constant by matching to rji = 0, which as we found before gives 
an action 2m(N + 1/2) - log A - 2. 

Before comparing these expressions with the asymptotics of the DOZZ formula 
in Regions I and II, we first comment on their multivaluedness. In order to do this 
we must determine the analytic structure of the function ^(77). It is clear from the 
definition (4.27) that F(j]) has branch points at each integer 77. The form of the branch 
points for 77 = —n with n = 0, 1, 2, . . . is — (77 + n) log(r] + n), while for 77 = m with 
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m = 1,2, . . . it is (77 — m) log(m — 77). We thus find that the monodromy of F(rj) around 
any loop in the 77-plane is 

00 00 
F( V ) -> F( V ) + J^iv ~ m)2mN m - ^(77 + n)2mN n , (4.32) 

m=l n=0 

where N n and iV m count the number of times the loop circles the branch points in 
a counterclockwise direction. Now applying this to (4.26), we see that continuation 
in the 77^ can produce far more branches than can be accounted for by nonsingular 
complex solutions. In particular, the various nonsingular solutions can only account 
for multivaluedness of the form 2iriN(l — Y^iVi)> while continuation around a loop in 
the general 77, parameter space can easily produce shifts of the action by terms such 
as 27iiN(r]i + 772 — 773). There thus seems to be a mismatch between the branches of 
the action (4.26) and the available saddle points. One might be tempted to interpret 
this multivaluedness as indicating the existence of additional solutions, but we showed 
in section 3.3 that there are no more solutions. We will suggest a mechanism for 
explaining this additional multivaluedness in section 5, as part of our discussion of the 
singularities that appear in the case of general rji, and another possible interpretation 
in section 6. The situation however is simpler for continuations that stay in Region I 
and/or Region II. Such a continuation will only activate the branch cuts in F(^2iVi ~ 
1), and this produces the kind of multivaluedness that can be accounted for by the 
known nonsingular solutions. In particular the action (4.31) can be gotten by analytic 
continuation from (4.26) along a path that goes from Region I to Region II, with the 
particular saddle point we land on being determined by the number of times the path 
wraps around 'YjiVi = 1- 

4.2.2 Comparison with Asymptotics of the DOZZ Formula 

We now compute the semiclassical limit of the DOZZ formula (2.28) with three heavy 
operators in Regions I and II. 26 

The semiclassical behavior of the prefactor of the DOZZ formula is clear: 

(Q-E^)/ fe 



A 7 (6 2 )6 



-2b 2 



y exp 




l> 2 



(4.33) 

To study the remaining terms, we need the b — > behaviour of T^ry/fe). In Appendix 
A, we show that 

T b (rj/b) = e ^[^)-(^i/2) 2 iog^(Mo gfe) ] Q < Re(r/) < L (4 34) 



26 This computation was done in [7] in Region I with real r\i, and our computations here are simple 
extensions of that. 
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In Region I, all of the Y^'s have their arguments in the region of validity for this formula, 
so we find that they asymptote to: 



exp 



b 2 



F{2r tl ) + F{2 V2 ) + F(2 m ) + F(0) 



(4.35) 



Combining these two contributions, we find complete agreement with (4.31) with A" = 
0. Thus in Region I, only one saddle point contributes and we can interpret the path 
integral as being evaluated on a single integration cycle passing through it. 

In Region II, the only new feature is that (£\ ot-i ~ Q) is no longer in the region 
where we can apply (4.34). To deal with this, we can use the recursion relation (A. 4) 
to move the argument back to the region where we can use (4.34): 

(4.36) 

Using also (4.10) for the asymptotics of the Gamma function (and hence of 7(x) = 
T(x)/r(l — x)), we finally arrive at 



C(r)i/b) ~ exp 



+ F( ^Vi) + F(rj 1 +7] 2 - m) + F{r]x + r] 3 - r? 2 ) + F(rj 2 + i] 3 - r}i) 



1 



E 



x 



(4.37) 



This is in complete agreement with (4.31), provided that as with the two-point function 
we interpret the final factor as coming from a sum over infinitely many complex saddle 
points. Rather as before, the saddle points that contribute are A" = { — 1, —2, . . . } 
when Im r? - 1) /b 2 ) < and N = {0, 1, 2, ... } when Im rji - l)/b 2 ) > 0. The 
condition l^a((YliVi ~ l)/^ 2 ) = defines a Stokes wall. 
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4.3 Three-Point Function with Light Operators 

So far we have considered only correlators where all operators are heavy. As a final 
check we will compute the semiclassical limit of the DOZZ formula (2.28) with three 
three light operators of Liiouville momenta on = b<Ji, with <Ji fixed for b — > 0, and 
compare it with a semiclassical computation based on equation (2.26). This compution 
is essentially a repackaging of a fixed-area computation outlined in [7]; we include it as 
an additional illustration of the machinery of complex saddle points and also because 
many of the details were omitted in [7] . In section 7, we will also use the same tools to 
do a new check in the context of timelike Liouville, so it is convenient to first present 
them in a more familiar context. 

We begin by computing the asymptotics of the DOZZ formula with three light 
operators; in order to capture the nontrivial effects of the operators we need to compute 
to higher order in b than before. To order b° in the exponent the pref actor not involving 
T/,'s becomes 



A 7 (6 2 )6" 



-2b l 



-£i«0/& 



& -2/6 2 +2£i <7i-4^1/6 2 +l-Ei <Ti e -2 7j5 +0(61og6) 



(4.38) 



Here je is the Euler-Mascheroni constant je = ^^n^oo w2k=i I ~~ logn). To take the 
limits of the functions, we need the the asymptotics of Tb(ab) as b — > 0. This is 
given by equation (A. 10): 



T b (ba) 



cb 1 ' 2 - 



exp 



^log6+^ + 0(6 2 log6) 



(4.39) 



Here C is a constant that will cancel in the final result. This along with (A. 11) is 
sufficient to determine the asymptotics of all parts of the DOZZ formula except for the 
Tfe involving Xli "*- For this one we can use the recursion relation: 

T 6 [(X> i -l)6-l/&] =7(E^- 1 - 1 ^ 2 )^ 3+2 /" 2 - 2 ^-T b [(^a l -l)6] (4.40) 

i i i 

To evaluate the semiclassical limit of this we need the corrections to (4.10). We can 
get these by using the machinery of Appendix C, but we can simplify the discussion 
using Euler's reflection formula r(x)T(l — x) 



7(x - l/b 2 ) 



sin(7rx) ' 
7T 



(4.41) 



r(l -x+ l/6 2 ) 2 sin (tt(x- l/b 2 )) 

The T function appearing on the right hand side of this equation always has positive 
real part as b — > 0, so we can simply include the first subleading terms in Stirling's 
formula to find 

r(l - x + l/b 2 ) = V2^b- 2 / b2+2x - l e- l/i; \l + Oib 2 )). (4.42) 
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This then gives 



7($>-l -1/& 2 ) 



i 



b^ 2 -^^e 2 ' b \l + 0{b 2 )). 



(4.43) 



Combining all these results together we can write: 




r(<7i + Q-2 - o- 3 )r(o-i + 03 - cr 2 )r(er 2 + o- 3 - Q-i)r(q- 1 + q- 2 + or 3 - 1) 

r(2 ( x 1 )r(2 r 2 )r(2 ( T3) 

(4.44) 



We now compare this result to an appropriate refinement of (2.26). There are 
several subtleties to consider. With all operators light the appropriate saddle point 
is the sphere (4.7). As with the saddle point (3.17) for the two-point function with 
heavy operators there will be a moduli space of such solutions, in this case given 
by the quotient SL(2,C)/SU(2), since the subgroup of SX(2,C) that leaves fixed a 
particular round sphere metric is a copy of SU(2). The light operator insertions will 
depend explicitly on these moduli, so we need the general SL(2, C) transformation of 
the saddlepoint (4.7). From (2.2) this is given by 



with a, {3, 7, 5 G C and obeying a5 — f3j = 1. In using (2.26) we will need to integrate 
the right hand side over all such saddlepoints. 

An additional subtlety is that in (2.26) all effects of the operator insertions are 
O(b ) in the exponent. To precisely include all effects of this order, we would need to 
carefully compute the renormalized fluctuation determinant about each saddle point, 
and also include the O(b ) corrections to the action (2.10). Moreover we would need 
the Jacobian in transforming the integral over <ft c into an integral over the parameters 
a, (3, 7, 5. We will include the subleading terms in the action explicitly, but to simplify 
things we will represent the fluctuation determinant and Jacobian as a 6-dependent 
prefactor A(b) which is at most (9(log&) in the exponent. 27 Note that neither of these 
things should be affected by shifting the saddlepoint by 2iri so we expect A(b) to be 
independent of N. It is also independent of <j\ since neither effect has anything to do 

27 We do NOT need to include 0(b 2 ) corrections to the saddlepoint (4.45) even though they are 
present. The reason is that the leading order saddlepoints are stationary points of the leading order 
action, so perturbing the solution at 0(b 2 ) does not affect the action until 0(b 2 ), which is beyond our 
interest. 



(j) C:N {z, z) = 2ni(N + 1/2) - log A - 2 log (\az + /3| 2 + \jz + 5\ 2 ) , 



(4.45) 
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with the operator insertions. With this convention, we can now write a more precise 
version of (2.26) that is appropriate for comparison with (4.44): 



3 



(V bai ( Zl ^i)Vba 2 (z 2 ,z 2 )V ba;i (z 3l z 3 )} « A(6) J2e~ SL[ ^ N] [ d(i(a,P,j,S)f[e at 

N£T J i=l 

(4.46) 

Here T is some set of integers and 

dfi{a, (3, 7, 5) = 45 2 {a5 -07-1) d 2 a d 2 (5 d 2 j d 2 5 

is the invariant measure on SL(2, C) [7]. The integrals over over the full a, /3, . . . planes. 

The G(b°) correction to the action (2.10) is given by ^- § dD (p c dO + 41ogi?. For the 
saddle point (4.45) the leading part was computed above (4.8), and now including the 
subleading term we find 



Sl[4>c,n] — jn 



2ni(N + l/2) -logA-2 + 2ni(N + 1/2) - log A + Q{b 2 ). (4.47) 



The integral over the moduli is quite difficult, we will simplify it some here and then 
relegate the final computation to an appendix. Our technique is identical to that in 
[7]. We first note that 

r 3 

/ dfi(a, (3, 7, 5) Y[ e^^*"*) = A" E * <V i(7V+1/2) £* 



dn(a,P,y,8) 



/ \ 2<T\ / \ 2(72 / \ 2(T3 ' 

(^laz! + /?| 2 + + <5| 2 J (\az 2 + P\ 2 + \ 1 z 2 + 6\ 2 j (\az 3 + (3\ 2 + \yz 3 + 5\ 2 j 

(4.48) 

The position dependence of this integral can be extracted by using its SL(2, C) trans- 
formation properties; changing variables by the transformation which sends z\ — > 0, 
z 2 — > 1, and z 3 — > oo we find the usual three-point function behaviour 

r 3 

/ dp Yl e ^ c,N = X ~ El <V 7rl(JV+1/2) E < °i\ Zl2 \ 2 te-°i-^) l^l^i-^-as) \ Zl3 \2{*2-<n-*3) i( au(T2} a3 ) 
J i=i 

(4.49) 

with 



I{a 1 ,a 2 ,a 3 )= / ^— ^— (4.50) 

|/3| 2 + \6\ 2 j (\a + /3| 2 + | 7 + 5| 2 ) (|a| 2 + | 7 |2 
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The result of this integral was quoted in [7] , but many steps were omitted and the full 
evaluation is quite sophisticated. For completeness we have included a full derivation 
in Appendix F. The result is 



Comparing this with the DOZZ asymptotics (4.44) we find complete agreement, with 
the saddle points included depending on the sign of ImQ^fli — V^ 2 )- We also see 
that apparently A(b) = i7r~ 3 fe~ 3 e~ 2TE , which would be interesting to check by explictly 
treating the measure. That it is imaginary is unsurpising given the complex integration 
cycle. 

4.4 Summary 

This concludes our argument that the analytic continuation of the DOZZ formula 
in Regions I and II is described by the Liouville path integral evaluated on a complex 
integration cycle that changes as we cross Stokes lines. The behaviour is completely 
analogous to that of the Gamma function as described in Appendix C. This has a 
qualitative explanation that was explained in section 1.1.1. The integral representation 
of the Gamma function is the zero mode part of the Liouville path integral, and the 
complex saddle points that we studied for Regions I and II differed only by shifting the 
zero mode. What we learned in this section is that in Regions I and II there are no 
additional subtleties in the analytic continuation in rji beyond those that are already 
apparent in the zero mode. 

5 Four-Point Functions and the Interpretation of Singular Sad- 
dle Points 

We now confront the issue first raised in section 3.3: for most complex values of the 
r)i, there are no nonsingular solutions of Liouville's equation with the desired boundary 





x 



r(o-i + o- 2 - o~3)r(o-i + Q- 3 - o- 2 )r(o- 2 + 0-3 - o~i)r(o-i + o- 2 + q~3 - 1) 

r(2 ( T 1 )r(2 ( x 2 )r(2a3) 



(4.52) 
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conditions. The candidate solution (3.25) fails to be a solution because of zeroes of the 
denominator function 

f{z,z) = A{z-z ) + B{z-z ) + .... (5.1) 

At such a zero, <p c = —2 log / — log A is singular, and perhaps more seriously, it is also 
generically multivalued. Around a zero of / with winding number k, <p c changes by 
—4mk. 

This seems to raise a serious challenge to any attempt to interpret the full analytic 
continuation of the DOZZ formula in terms of conventional path integrals. In this 
section we will study this further. We will make three arguments that even when C is 
multivalued, the expression (3.25) still makes some sense and controls the asymptotic 
behaviour of the DOZZ formula. We will first show that there is a minor redefinition of 
the action which agrees with the formula (2.25) when there are no singularities but is 
finite even in the presence of zeroes of the denominator. Moreover it correctly produces 
the analytic continuation of (2.25). We will then show that the presence of singularities 
actually allows the full multivaluedness of the action (4.31) to be realized by analytic 
continuation of the "solutions." Finally we will probe the saddle points that dominate 
the three-point function by including a fourth light operator. For the case that we are 
able to implement this test - the case that the light operator is degenerate - we will 
find agreement with the (3.25) for all values of the rji. We will close by commenting on 
the implications for general four-point functions. 

5.1 Finiteness of the "Action" 

We begin by observing that in Region I defined at the beginning of subsection 
4.2, we included a restriction on the imaginary parts of the r^'s to ensure that the 
denominator in (3.25) did not vanish away from the operator insertions. However, the 
formulas that followed seemed to know nothing about this additional restriction; the 
multivaluedness in the expressions for Cj and Sl cannot be activated without violating 
the conditions ReQ^T^) < 1 or Rer/j < |, regardless of the imaginary parts of the 
rji's. Moreover the expression (4.31) for the action can easily be continued to values 
of rji where the denominator vanishes, and its value is perfectly finite there. This is 
perhaps unexpected because near a zero of the denominator, one has 

<f> c (z,z) « -2log[A{z-z ) + B{z-z )), (5.2) 

which has a logarithmic singularity as well as a branch cut discontinuity. 28 With such 
discontinuous behavior, the kinetic term in the Liouville action J cP£<9 a c <9 a c certainly 

28 We consider the case that \A\ ^ \B\, which is generically true for complex ?y's. When the rji's are 
real and a\ is also real then we can have \A\ — \B\, we will comment on this below. 
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diverges. The finite analytic continuation of the action therefore cannot be computed 
by naive application of (2.25). 

We begin by observing that for solutions with no additional singularities we can 
rewrite (2.25) as 



As before, the c?j are small discs centered around Z{. We propose that even in the 
presence of zeroes of the denominator of e^ c , this is still the correct form of the action, 
with the integral defined by removing a small disc of radius e centered around each 
zero and then taking e — > 0. The divergence from the discontinuity in C is avoided 
since / is continuous, but we still need to show that there is no divergence as e — > 0. In 
particular near a zero at z — Zq, we have the expansion (5.1), so we can approximate 
the contribution to the integral from the vicinity of Zq as 



The radial integral is logarithmically divergent, but as long as \A\ ^ \B\ the angular 
integral is zero! The higher order corrections to / will produce manifestly finite correc- 
tions to the action, and in fact one can show that this definition of the action is invariant 
under coordinate transformations of the form z — > z + 0(z 2 ). This is thus analogous 
to the principal value prescription for computing the integral of 1/x across x = 0. We 
claim that the action computed this way agrees with what one gets by analytic contin- 
uation in rji. To justify this, we need to show that we can continue to use the trick of 
differentiating with respect to rji to calculate the action. This requires a demonstration 
that a multivalued "solution" is a stationary point of the improved action. To show 
this we can compute the variation of the improved action under /—>•/ + 6f with 5f 
continuous; most terms are clearly zero when evaluated on a multivalued "solution," 
but a potentially nontrivial boundary term is generated by the integration by parts: 



For intuition, we observe that this boundary term is also present near each of the 
operator insertions. Near the operator at Zj, we have / ~ r 2r,t , and the boundary term 
produces a nontrivial variation —-$%5f. This variation is cancelled by the variation 
— 2n it) 2 ' * d0(fi c of the regulated operator. The point however is that for / obeying (5.1), 
this boundary term is automatically zero by itself since the angular integral vanishes. 
So in this sense, a multivalued "solution" is a stationary point of the action. 




d 2 £ [dfdf/f 2 + l/f 2 ] + boundary terms. 



(5.3) 




(5.4) 




(5.5) 
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For orientation, perhaps we should mention that a singlevalued <p c with singular- 
ities away from the operator insertions can never be such a stationary point. Indeed, 
generalities about elliptic differential equations ensure that a solution of the complex 
Liouville equations is smooth away from operator insertions. In the singular case, it is 
only because <ft c is multivalued that it may be, in some sense, a stationary point of the 
action. 

With this explanation of what the action means in the presence of singularities, 
we may drop the conditions on the imaginary parts of rji from both Regions I and II 
and the story of the previous section goes through unchanged. This argument does fail 
in the special cases where \A\ = \B\, for which higher order terms near the singularity 
are important and the singularity may be non-isolated. We will view this just as a 
degenerate limit of the more general situation. In particular we can continue from 
Region I to anywhere else in the ^-plane without passing through a configuration with 
a singularity with \A\ = \B\, so this subtlety should not affect our picture of the analytic 
continuation of (2.28). 

Before we move on, we observe that there are two different kinds of multivaluedness 
being discussed in this section. One is with respect to rji, and the other is with respect 
to z,z. For convenience we summarize the multivaluedness of various quantities in the 
following table: 





z, z behaviour at fixed rji 


rji behaviour at fixed z, z 


g</>c 


singlevalued 


singlevalued 


b 2 S L 


trivial 


defined up to addition of 

27ri(Y2iVi m i~\~ n ) with rrii all 
even or all odd 


Q 


trivial 


defined up to addition of 2ni 




trivial 


defined up to multiplication 
by a sign 


a 1 /a 2 


trivial 


singlevalued 


f 


singlevalued 


defined up to multiplication 
by a 2, ^-independent sign 


4> c 


possibly singlevalued, possi- 
bly monodromy of addition 
of 4Tri about points where 

/ = o 


defined up to addition of 2ni 
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5.2 Multivaluedness of the Action 



We saw in subsection 4.2 that the action (4.26) is highly multivalued as a function 
of the rji, with the multivaluedness arising from the function F{rf) defined in (4.27). 

We can now interpret this multivaluedness of the action as a consequence of the 
multivaluedness in z,z that <ft c acquires in the presence of zeroes of /. This multi- 
valuedness does not affect the kinetic and potential terms of the action as defined in 
section 5.1, since they depend only on /, which is singlevalued as a function of z, z. 
But the terms — Y2i ^ d6(p c that come from the regulated operator insertions are 
sensitive to this multivaluedness. Their contribution to the action is 



where C, is the constant term in C near the operator insertion. Using the formulas 
(4.22)-(4.24) for Cj, we see that continuing along a closed path in the parameter space 
of the rji can shift Cj by an integer multiples of 2Tri, hence shifting the action by an 
integer linear combination of the quantities 27rii]i. We can see the same effect in the 
formula (4.26) for the action; the same processes that cause a shift in the Ci cause an 
equivalent shift in the function F in this formula, leading to the same multivaluedness. 
For example, on a path on which f]i + r] 2 — r] 3 circles around an integer value, shifting 
C% and C*2 by 2iri and C% by —2iri, there is a corresponding shift in the action from 
F(ih -h 772 — 773)- 

It is important to note that it is only because <p c can be multivalued as a function 
of z, z that we can realize the full multivaluedness of the action in rj. We argued below 
equation (4.24) that any continuation in rji that passes only through continuous c 's 
cannot produce monodromy for the difference of any two Cj's because of continuity. 
But once we allow paths in rji that pass through multivalued (and thus discontinuous) 
c 's, these differences can have the nontrivial monodromy necessary to produce the 
full set of branches of the action. Thus the multivaluedness of the action in rji has a 
natural interpretation once we allow solutions of the complex Liouville equations that 
are multivalued in z, z. 

5.3 Comparison With The DOZZ Formula 

We are finally ready to consider in general the semiclassical asymptotics of the 
DOZZ formula (2.28). The DOZZ formula is constructed from the function Y^/fe), 
where 77 is a linear combination of the t]i. In all, seven Tf, functions appear in the 
numerator or denominator of the DOZZ formula. To evaluate the small b asymptotics 
of this formula, one needs the small b asymptotics of the functions. This is given in 




(5.6) 
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(A. 7) for i] in a certain strip in the complex plane; it can be determined in general by 
using the recursion relations (A. 3) to map rj into the desired strip. In the process, the 
recursion relation generates a function that can be expanded as a sum of exponentials, 
as in (4.14); we interpret this as a sum over different complex critical points. 

For generic rji, when evaluating the asymptotics of the DOZZ formula using the 
asymptotic formula (A. 7), we will need to apply the recursion relations to all of the 
T&'s. There is just one crucial difference from the derivation of eqn. (4.37). The final 
factor in that formula has an expansion in positive or negative powers of exp(27r? Yli r li)i 
where ^ rji entered because in that derivation, we had to apply the recursion relation 
only to one of the Y& functions, namely T^Q^^/ft). In general, we have to allow for 
the possibility that the argument of any one of the seven Y& functions in the DOZZ 
formula may leave the favored strip. So Y2i Vi ma ¥ t> e replaced by the equivalent 
expression appearing in any one of the other Y& functions, namely 2r)i, 771+772 — ?73, or 
any permutation thereof. 

In the process, it is not quite true that the action can be shifted by 2ni £\ for 
arbitrary integers TTij. Rather, the rrii are either all even or all odd. This holds because 
similarly the Y& functions in the DOZZ formulas are all functions of Yli c iVi/b, where 
the Cj are all even (the factors in the numerator of the DOZZ formula) or all odd (the 
factors in the denominator). 

5.3.1 A Further Comment 

One interesting point about this is that for some values of the rji, singlevalued complex 
solutions of Liouville's equations do exist. But even in such regions, we may need to use 
the recursion relations to compute the asymptotics of the DOZZ formula, and hence 
we seem to need the full multivaluedness of the action, even though from the present 
point of view this multivaluedness seems natural only when the classical solutions are 
themselves multivalued. The reason that this happens is that in continuing in rji from 
Region I to these regions we necessarily pass through regions where <p c is multivalued in 
z. When we arrive at the region of interest it is then possible that although a continuous 
single-valued solution exist we have actually landed on a discontinuous one. 29 The 

29 For a simple example of this phenomenon, consider the function 



For r\ real and positive we can define the branch of the logarithm so that h is a continuous function 
with an ambiguity of an overall additive factor of 2niN . But if we choose such a branch and then at 
each point x continue in r\ around a circle containing rj = 0, this will produce a shift of 2iri near x = 1 
but not near x = 0; the resulting function will thus be discontinuous even though a continuous choice 
of branch exists. 
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locations and strengths of these discontinuities will depend on the path in 77. This 
allows the full multivaluedness of the action to be realized, since the discontinuities will 
not affect the kinetic term when written in terms of / but they will allow independent 
shifts of 4> c by 2mN near the operator insertions and infinity. 

These discontinuities are admittedly unsettling so we note here that in section 6, 
we explain a different point of view in which the full multivaluedness of the action is 
equally natural for any values of the rji. 

5.4 Degenerate Four-Point Function as a Probe 

The previous two arguments for the role of multivalued "solutions" in the Liouville 
path integral were rather indirect. We give here a more direct argument. In section 
2.4, we reviewed Teschner's formula (2.44) for the exact four-point function of a light 
degenerate field VLj/2 with three generic operators V ai . This expression is meromorphic 
in «j, and choosing all three cVs to scale like 1/6 we can study its semiclassical limit 
for any values of the fy. Moreover we can compare this to (2.26), which says that in 
the semiclassical limit this correlator can be evaluated by replacing the operator 
by the function exp(— b(f>) = exp(— </> c /2), where 4> c is the saddle point determined by 
the three heavy operators. If there are several relevant saddle points 4> c ,n, iVeT, with 
action Sl,n, then (2.26) gives 30 

(V m / b (z 4 , Z4)V m/b (z 3 , Z 3 )V^ b / 2 (z2, Z2)V m/h {zi, Zi)) 

~ e ~(l>c,N{z2,Z2)/2 e -S L , N 

N 

Using the definitions (2.27) and (2.32) and also (3.2), this implies 

G l234 (x,x) « Vx j Zu jj Zu j 2 ^2f N (z2,z 2 )e-^. (5.7) 



N 



Si34,n is a branch of (4.26) without its position-dependent terms, with the branch 
labelled by N, and with the replacement 772 — > 774. Explicitly: 

b 2 S 134 ,n = {r)i + 773 + 774 — 1) log A + F(t)i + 774 - 7/3) + F(?7i + 773 - 774) + F(r) 3 + r] A - rjx) 
+ F( Vl + 773 + 774 - 1) - F(2 Vl ) - F(2 m ) - F(2t7 4 ) - F(0) 
+ 2m{n + mirji + m 3 r] 3 + 777,4774). (5.8) 



30 As discussed below (2.26), we have omitted Z2-independent factors that are O(b ) in the exponent. 
These come from the functional determinant and corrections to the action (2.25). These factors will 
cancel between the two sides in (5.14) below. 
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Here n, mj are integers determined by the branch N. We saw in section 3.3 that 
1/ fx is uniquely determined (independent of N), so is uniquely determined 



up to sign. (The sign comes from the choice of square root in defining a\.) By comparing 
this with the semiclassical limit of Teschner's formula (2.44), we can thus explicitly 
check the position dependence of the saddle point f(z, z)\ 

Rewriting Teschner's proposal (2.44) with a condensed notation, we have 



G 



1234 1^, X ) 



12^34- 



J r _(x)J r _(x) 



C + 12C34+ 
C~ 12C34- 



3F + (x)3F + (x) 



(5.9) 



Teschner's recursion relation can be rewritten as 
C*34+C + 12 1 



C34-C Yl 

X 



1 - 6(2ai - b)Y 

7 2 (26(2cn - l))7(6(a 3 + a 4 - <*i - b/2)) 



7(6(«i + a 4 - a 3 ))7(6(ai + a 3 - a 4 - b/2))j(b(ai + a 3 + a 4 - Q - b/2)) ' 
and using (3.22) and (3.23) and taking the semiclassical limit this becomes 

C34+C" 



(5.10) 



12 



Cm- C 



— >■ 



12 



a 2 
ai 



(5.11) 



Here a>\ and a>2 are the constants in the semiclassical solution (3.25), with the replace- 
ment r] 2 — > rj^. In the same limit, we can see from (3.21) that 



F+(x) P l - rn {x) 
T- (x) P Vl (x). 



(5.12) 



In checking this, it is useful to recall that we can send 773 — )■ 1 — 773 in the definition 
of P 1 ~ m (x) without changing the function since this is one of Kummar's permutations 
from Appendix B. We thus find that in the semiclassical limit we have 



G 



1234 (.^j x ) 



C 



34- 



P r,1 (x)P 



m 



[x) - — P 1 - r > 1 (x)P 1 - r > 1 (x) + 0(b) 
ai 



(5.13) 



With the help of (3.25), we find that this will agree with (5.7) if 



31 



-S34_ 



\z\a\ 


^34 






Zl3 





Sl34,N+0(b) 



(5.14) 

31 In deriving this formula, we neglected O(b ) terms in the exponent of C34-. These are the same 
terms that we previously neglected on the right-hand side of (5.7), since the difference between rji and 
rji — b/2 affects them only at subleading order. So this equation really needs to be true to order O(b ) 
in the exponent. 
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Beginning with this equation we explictly include the branch dependence of a\ for the 
rest of the section. Semiclassically the structure constants C134 and C34- are in the 
same region of the r\i plane since their r/j's differ by something that is 0(b 2 ), so we 
can assume they are both a sum over the same set of branches N. This justifies our 
equating the sums term by term in (5.14). Using (5.8), we see that: 



SlU,N ~ S 34 -,N =- log A + i-KTTli + - 



log 7(771 + 774 - 773) + log 7(7/1 + 7/3 - 7/4) 
+ log 7(7/1 + 7/3 + 7/4 - 1) - log 7(7/3 + 7/4 - 7/i) - 2 log 7(2?7i] 

+ 0(6). (5.15) 

Comparing with (3.24), we see that (5.14) is clearly satisfied up to an overall branch- 
dependent sign. 

To see that this sign works out, we need to give a more careful argument. First we 
can define 

log 7(771 + 7/4 - 7/3) + log 7(7/1 + 7/3 - 7/4) + log 7(7/1 + 7/3 + 7/4 - 1) 

- log 7(7/3 + 7/4 - 7/i) - 2 log 7(27/i) + 77T7711 . (5.16) 





Zl3 






^34 



oijv = 1 — n — r ex P 



The logarithms are defined by continuation from real 77's in Region I along a specific 
path, which gives an unambiguous meaning to rh\? 2 The signs will match in (5.14) if 
mi = 777i. To demonstrate this, recall that near z\ we may write 

(p c ,N = -47/ilog \z - zi| + C 1>N , (5.17) 

with 

Ci N = ~ 27T77711 - log A - (1 - 27/i) log log 7(7/1 + 7/4 - 7/ 3 ) 

P34r 

- log 7(7/1 + 7/3 - 7/4) - log 7(7/1 + 7/3 + 7/4-I) 

+ log 7(7/3 + 7/ 4 - 7/1) + 2 log 7(27/1). (5.18) 

Here the logarithms are defined by analytic continuation along the same path as in 
defining a^jv- Since d ^ L '* = — Ci,/y, we are justified using mi in this formula. Finally 
near z = Z\ we have 

so in (5.7) we should choose the branch of fjy, and thus of a^jv, with rh\ = mi. 

This completes our demonstration of (5.14). We consider this to be very strong 
evidence that at least for the case of the degenerate four-point function, the Liouville 
path integral is controlled by singular "solutions" throughout the full rji three-plane. 



:!2 



It does not matter what the path is, but we need to choose one. 
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5.5 Four-Point Function with a General Light Operator 

The discussion of the previous section showed that a certain type of four-point func- 
tion is semiclassically described by singular "solutions" of Liouville's equation. More 
specifically, the nontrivial position dependence of the correlator (5.7) was captured by 
the function fwfa, z 2 )- The effect of the singularities is rather benign, however; the 
correlator simply has nontrivial zeros as a function of the position of the light operator. 
As argued at the end of section 3, the zeros of Jn are generically stable under quantum 
corrections and thus are actually zeros of the exact four-point function (2.44). There is 
nothing inherently wrong with such zeros, but this observation is troubling nonetheless. 
The reason is that these zeros are smooth only because the light operator is exactly 
degenerate. If instead of the operator e~^ c ^ 2 we had considered a more general light 
operator e CT ^ c , then a semiclassical computation based on equation (2.26) (the other 
three operators are still heavy) would have given 

\ 7 1 4cr I Via 

G 12U (x,x) « G \-° \ 24 J y>(^ 2 r 2j e-W (5.20) 

1^34 | | ^14 | ^ 

Here Go is a O(b ) factor from the fluctuation determinant and the corrections to 
the action, both of which we expect to be independent of z 2 , and Si^n is given by 
(5.8). The problem however is that in the vicinity of a point Zq where fN(z 2 ,z 2 ) ~ 
A(z 2 — Zq) + B(z 2 — Zq), this correlator is generically singular and discontinuous! 33 We 
can quantify the nature of these singularities by using the winding number introduced 
at the end of section 3, and we find that the semiclassical correlator has winding number 
— 2cr around Zq if \A\ > \B\ and winding number 2a if \A\ < \B\. The winding number is 
not an integer because the function is discontinuous. It cannot be changed significantly 
by small corrections, and since it is generically nonzero we are tempted to conclude 
that the exact four-point function must also be discontinuous as a function of the light 
operator position at finite but sufficiently small 6! 34 

This situation would not be entirely without precedent; in the SL(2,M.) WZNW 
model appropriate for studying strings in AdS$ [38-41] it was shown in [42, 43] that 

33 0nc might hope that the discontinuity could cancel in the sum over the different branches N, but 
this will not work because for any given generic values of r\\ , r\2 , V3 > a there will be a single dominant 
saddlepoint that is parametrically larger as b — > 0. 

34 In general it is of course possible for a smooth function to have a semiclassical approximation 
which is discontinuous, a simple example is r(x/\) > which has a nne of zeros turn into a branch cut as 
A — >• 0. A more sophisticated example that we have been studying extensively in this paper is Tb(x/b), 
which exhibits the same phenomenon. That this does not happen for the four-point function under 
consideration is a special consequence of the semiclassical formula (5.20) for the correlator, where the 
nontrivial ^-dependence is all in a factor that is finite as b — > and the factor that goes like e~ x l h 
is independent of z%. 
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the exact 4-point function of certain operators has singularities when all four operators 
are at distinct positions. This could be seen semiclassically from stringy instantons 
going "on-shell" and was reproduced exactly using the machinery of the Knizhnik- 
Zamolodchikov equation [44] . In that situation however the singularities were localized 
to isolated points and the correlator was continuous away from those points. In the 
remainder of this section we will give an argument that in Liouville there are in fact no 
singularities, isolated or otherwise, in the exact four-point function when the operator 
positions do not coincide. We will then close the section with some speculation about 
where our semiclassical argument goes wrong. We caution however that we will use 
some plausible pieces of lore that have not strictly been proven, so our argument is 
slightly heuristic. 

We take as a starting point the exact formula (2.39) for the Liouville 4-point func- 
tion, which we reproduce here for convenience 



This formula is strictly true only when Re(«i + c^) > Q/2 and Re(«3 + 04) > Q/2. 
Away from this region, which we will certainly be with three heavy operators obey- 
ing the Seiberg bound and one light operator, there are additional discrete terms that 
are residues of the finite number of poles that have crossed the contour of integration. 
Looking at this expression, we see that there are only two possible sources of singu- 
larities in x, x. The first is singularities of the Virasoro conformal blocks J-1234 as a 
function of x, and the second is possible divergence of the integral over P for particular 
values of x. We will address each of these issues, beginning with possible singularities 
of the conformal blocks. 

The conformal blocks are expected to have branch points at x — 0,1, 00, which 
correspond to the UV singularities of the correlator when the operator at Z2 approaches 
the operators at z\, Z3, or Z4. The singularity at x = is manifest from the definition 
(2.40), and the singularity at x = 1 arises from the nonconvergence of the series in 
(2.40) when \x\ = 1. When all operator weights are real and positive the fact that 
the radius of convergence of this series is indeed one follows from the convergence of 
inserting a complete set of states in unitary quantum mechanics. The convergence 
for generic complex operator weights has actually never been proven in the literature, 
although it was conjectured to be true in [45] and discussed more recently in [30, 46]. 
In [46] it was proven that if the radius of convergence is indeed one, then there are no 
other singularities with \x\ > 1 except for the singularity at infinity. We will also not 
be able to prove this convergence, but we give two pieces of evidence in favor of it. 




C(ai, a 2 , Q/2 + iP)C{a 3 , a 4 , Q/2 - iP) 
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First we note that the c — > oo limit of the conformal block, which turns out to mean 
including only descendants of the form (L_i) n |Q/2 + iP) in the sum in (2.40), can be 
evaluated explictly from the definition and gives [45, 47] 

lim ^1234 (Ai, A P , x) = x Ap - Al - A2 F(A P + A 2 - A 1; A P + A 3 - A 4 , 2A P , x). (5.22) 

c— >oo 

As discussed in Appendix B, this hypergeometric function is singular only at x — 
0, 1, oo. So any additional singularities of J-1234 would have to disappear in the c — > 00 
limit, which seems unnatural. When c is finite but one of the external legs is degenerate 
we can again compute the conformal block, with result (2.48). Again it only has 
singularities in the expected places. In 27 years of studying these functions as far as we 
know no evidence has emerged for singularities at any other points in x, so from now 
on we assume that they do not exist. 

The other possible source of singularities in the four point function (5.21) is diver- 
gence of the integral over P. To study this further, we need large P expressions both 
for the structure constants and the conformal blocks. The appropriate asymptotics for 
Tfe are quoted (with some minor typos we correct here) as equation 14 in [30]: 

3 17V 

logTb(x) = x 2 logx — -x 2 —x 2 + O(xlogx) Imx — )■ ±00. (5.23) 

We will not derive this, but it isn't hard to get these terms from our expression (4.34) 
for the semiclassical limit of T& with x scaling like 1/6. 35 Using this in (2.28), we find 
that at large real P we have 

C(a u a 2 , Q/2 ± iP) = iQ- p2 +v( p ^ p ) . (5.24) 

The structure of the conformal blocks at large P was studied by Al. B. Zamolodchikov 
in a series of papers [48, 49] , 36 ; he obtained the following remarkable result: 

J*i23 4 (Ai, A, x) =(16g) A -^a;^- Al - A2 (l - x)^- A2 ~ A3 

X 6 3 (q) £ f i -4(Ai+A2+A 3 +A 4 ) ^ + / A )) _ (535) 

35 To do this, we observe that for large 77 we can use Stirling's formula to approximate log 7(2:) inside 
the integral expression (4.27) for F(r/). This is not quite the same as a full finite-^ derivation since in 
principle there could be subleading terms in b that become important for sufficiently large 77, but we 
have checked this formula numerically at finite b with excellent agreement so apparently this does not 
happen. The formula (4.34) was valid only for 77 in a certain region, but using the recursion relations 
to get to other regions will not affect things to the order we are working in (5.23) so (5.23) is valid for 
arbitrary Re(x). 

36 English translations are availiable online but hard to find. This also especially the case for reference 
[45], which gives a beautiful exposition of the general formalism of [17] that is more complete than 
anything else in the literature. The most accessible place to find the formula quoted here seems to be 
in section 7 of [7] , but beware of a notational difference in that our conventions are related to theirs 
by 1 <-> 2. 
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Here 6 3 (q) = X]^L-oo ^ an< ^ 1 = ex P [ — TfK{\ — x)/K(x)], with 

Ki - X) ~ 2 / (5 ' 26) 

This q can be interpreted as exp(i7rr), where r is the usual modular parameter of the 
elliptic curve y 2 = t(l — t)(l — xt). So in particular Imr is always positive and one 
always has \q\ < 1 when the elliptic curve is smooth (that is, for x ^ 0,1, oo). For 
fun we note that, like most things in this paper, K(x) is actually a hypergeometric 
function: equation (B.18) gives K(x) = § F(l/2, 1/2, 1, x). 

The derivation of (5.25) uses certain reasonable assumptions about the semiclassical 
limits of correlation functions; we will be explicit about them in Appendix D, where 
for convenience we review the origin of the leading behaviour 

F 12U (Ai,A P ,x) ~ (16g) P \ (5.27) 

This will be sufficient for our study of the integral in (5.21); combining it with (5.24) 
we see that the integral will converge, for i ^ 0,1, oo, given the fact that \q\ < 1. Thus 
the integral cannot generate any new singularities. This completes our argument that 
the Liouville four-point function (5.21) cannot have any new singularities in x. 

So what is wrong with our semiclassical argument for such singularities in the 
beginning of the section? To really understand this we would have to compute the 
semiclassical limit of (5.21) and compare it to our formula (5.20). For the moment, this 
is beyond our ability. We may guess however that the problem lies in our assumption 
that the factor Go is independent of z-x- This was true for the degenerate computation 
in the previous section, but the singularity we discovered here perhaps suggests that 
more sophisticated renormalization of the nondegenerate light operator is required in 
the vicinity of any singular points of the "solution". It is at first unsettling that the 
renormalization of the operator should depend on the positions and strengths of the 
other heavy operators, but we already saw in section 5.1 that even the "principal value" 
prescription for evaluating the action depended on these things at distances arbitrarily 
close to the singular point. Thus we expect that once an appropriate renormalization 
is performed, the semiclassical singularity in (5.20) will be smoothed out. It would be 
good to be more explicit about what this renormalization is, but we will not try to do 
so here. 

A different perspective on this four-point function is provided by the Chern-Simons 
formulation of Liouville theory, which we will introduce momentarily. In this formu- 
lation it seems clear that there are conventional nonsingular solutions that exist for 
any r\i and which can be used to study the semiclassical limit of this correlator; in this 
version of things it seems apparent that no singularity can emerge. 
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6 Interpretation In Chern- Simons Theory 



6.1 Liouville Solutions And Flat Connections 

In section 3.1, to a solution of Liouville's equations we associated a holomorphic 
differential equation 

' r +W(z))f = (6.1) 



( + W(z) ) / = 0. (6.2) 



dz 2 

and also an antiholomorphic differential equation 

/ d 2 
\dz 2 

Locally, (6.1) has a two-dimensional space of holomorphic solutions, and (6.2) has 
a two-dimensional space of antiholomorphic solutions. We constructed a solution of 

Liouville's equation from a basis ( U ) of holomorphic solutions of (6.1) along with 

W 

a basis \^~J °f antiholomorphic solutions of (6.2). This construction applies on any 

Riemann surface S, though we have considered only S 2 in the present paper. 

Globally, in passing around a noncontractible loop in E, or around a point at which 

there is a singularity due to insertion of a heavy operator, the pair ( | has in general 

W 

non-trivial monodromy. The monodromy maps this pair to another basis of the same 
two-dimensional space of solutions, so it takes the form 

:MIH(:). <«» 

where M is a constant 2x2 matrix. Actually, the determinant of M is 1, so M 
takes values in SL(2, C). One way to prove this is to use the fact that the Wronskian 
udv — vdu is independent of z, so it must have the same value whether computed in the 
basis u,v or the basis u, v. This condition leads to det M = 1. Alternatively, we may 
observe that the differential equation (6.1) may be expressed in terms of an SX(2,C) 
flat connection. We introduce the complex gauge field A defined by 

Since these 2x2 matrices are traceless, we can think of A as a connection with gauge 
group SL(2, C). 37 On the other hand, a short calculation shows that the condition for 



37 0ur convention for non-abelian gauge theory is that — + A^, so in particular 
[Dp, D v ] = d^Ay — <9„A M + [Ap, A v \ and the gauge transformation is — > gD^g^ 1 , with g 6 G. 
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a pair to be covariantly constant with respect to this connection is equivalent to 

w 

requiring that / is a holomorphic solution of the equation (6.1) while g = df/dz. Thus 
parallel transport of this doublet around a loop, which we accomplish by multiplying 
by U = p e -f A * dz , is the same as analytic continuation around the same loop. In 

particular if we define the matrix S — [ ™ _ ) , then we have US = SM T . Taking 

\ou ov J 

the determinant of this equation, we find that M G SL(2, C). 

Similarly, the antiholomorphic differential equation (6.2) has monodromies valued 
in SL(2, C). This may be proved either by considering the Wronskian or by introducing 
the corresponding flat connnection A, defined by 



A z = 0, A z =\ \ VV ^\. (6.5) 




(It is sometimes convenient to take the transpose in exchanging z and z, and we have 
done so, though this will not be important in the present paper.) 

The connections A and A have singularities near points with heavy operator in- 
sertions. The monodromy around these singularities can be inferred from the local be- 
havior of the solutions of the differential equation. For example, the solutions of (6.1) 
behave as z v and z l ~ v near an operator insertion at z = with Liouville momentum 
a = rj/b. The monodromies of these functions under a circuit in the counterclockwise 
direction around z = are exp(±27rz^). An invariant way of describing this, without 
picking a particular basis of solutions, is to say that 

TrM = 2cos(27rr]). (6.6) 

Similarly, the behavior of the local solutions near z = implies that the monodromy of 
the antiholomorphic equation (6.2) around z = has the same eigenvalues, and hence 
again obeys (6.6). 

More generally, the two flat connections A and A are actually gauge-equivalent 38 
and have conjugate monodromies around all cycles, including noncontractible cycles 
on E (if its genus is positive) as well as cycles of the sort just considered. This is 
guaranteed by the fact that / = uu — vv has no monodromy, since the Liouville field 
C is e*° = l/Xf 2 . 

So a solution of Liouville's equations - real or complex - gives us a flat SL(2, C) 
connection over E that can be put in the gauge (6.4) and can also be put in the 



:5N 



The explicit gauge transformation between them is g = 



I dudu + dvdv —udu — vdv 
\ —udu — vdv uu + vv 
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gauge (6.5). The basic idea of the present section is that, by a complex solution of 
Liouville's equations, we should mean in general a flat SL(2, C) connection, up to 
gauge transformation, which can be gauge-transformed to either of those two forms. 
We do not worry about what sort of expression it has in terms of a Liouville field. 

The attentive reader may notice that we have cut some corners in this explanation, 
because in section 3 the reference metric was chosen to be flat in deriving the holomor- 
phic differential equations. This is not possible globally if £ has genus greater than 
1, and even for S = S 2 , it involves introducing an unnatural singularity at infinity. A 
more precise description is to say that A is a flat connection that locally, after picking a 
local coordinate z, can be put in the form (6.5), in such a way that in the intersection of 
coordinate patches, the gauge transformation required to compare the two descriptions 
is lower triangular 



A flat connection with this property is known as an oper. This notion is explained 
in section 3 of [21], but we will not need that degree of detail here. The global char- 
acterization of A has the same form (with upper triangular matrices replacing lower 
triangular ones, given the choice we made in (6.5)). Our proposal then is that a classi- 
cal solution of Liouville theory is a flat connection whose holomorphic structure is that 
of an oper, while its antiholomorphic structure is also that of an oper. 

6.2 Some Practice 

A few elementary observations may give us some practice with these ideas. Let us 
first consider the main example of this paper, namely S 2 with insertions of three heavy 
operators of Liouville momenta a« = r]i/b, at positions z = z%. The monodromies Mj 
around the three points will have to obey 




(6.7) 





(6.8) 



In addition, the product of the three monodromies must equal 1: 



M 1 M 2 M 3 = 1. 



(6.9) 



Equivalently 



M 1 M 2 = M 3 



-i 



(6.10) 



from which it follows that 
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And of course we are only interested in a flat bundle up to conjugacy 



Mi^gM ig -\ geSL(2,C)- (6.12) 

To start with, let us just ignore the oper condition and ask how many choices of 
the Mi there are, up to conjugacy, that obey the conditions in the last paragraph. We 
can partially fix the gauge invariance by setting 

Mi = r n JL)- ( 6 - 13 ) 



e ~ 2niri \ 

The remaining freedom consists of diagonal gauge transformations 



»=[J A -.). < 6 - 14 > 



where we only care about A up to sign, since a gauge transformation by g = — 1 acts 
trivially on all gauge fields and monodromies. In general, we can take 

M 2 = fj (6.15) 

If we look for a solution with q = 0, we soon find that, for generic values of the rji, once 
we adjust p and r to get the right values of TrM 2 and det M 2 , we cannot also satisfy 
(6.11). So we take q ^ 0, in which case A in (6.14) can be chosen uniquely, up to sign, 
to set q — 1. Then, imposing det M 2 = 1, we get 

M 2 =( P X V (6.16) 

xPQ - 1 qj 

Now the conditions TrM 2 = 2 cos(27r?7 2 ), Tr M\M 2 = 2cos(27r?73) give two linear equa- 
tions for p and q which generically have a unique solution. So M 2 and therefore 
M3 = M 2 ~ 1 M 1 _1 are uniquely determined. 

The conclusion is that a flat bundle on the three-punctured sphere with prescribed 
conjugacy classes of the monodromies Mj is unique, up to gauge equivalence, even if 
we do not require the oper conditions. 39 The unique SL{2, C) flat bundle with these 
monodromies can be realized by a holomorphic differential equation and also by an 



This remains true for non-generic values of the rji where the derivation in the last paragraph does 
not quite apply, unless the r\i equal (0,0,0) or a permutation of (0, 1/2, 1/2). To verify this requires 
only one subtlety: if for some i, TrMj = ±2, so that the eigenvalues of Mi are equal, then one should 

not assume that Mi can be diagonalized; its Jordan canonical form may be ± ( ^ 
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antiholomorphic one. The proof of this statement is simply that functions W and W 
with the right singularities do exist, as in (3.18). 

Since it can be realized by both a holomorphic differential equation and an anti- 
holomorphic one, the unique SL(2, C) flat bundle on the three-punctured sphere with 
monodromies in the conjugacy classes determined by the rji is a complex solution of 
Liouville's equations in the sense considered in the present section. What one would 
mean by its action and why this action is multivalued will be explained in section 6.3. 

It is instructive to consider a more generic case with s > 3 heavy operators, with 
parameters rji, inserted at points Zj G S 2 . Now there are s monodromies Mj, i = 1, . . . , s. 
They are 2x2 matrices, constrained by 

Tr M { = 2cos(2vr?7 i ), det M; = 1. (6.17) 

We also require 

M 1 M 2 ...M S = 1, (6.18) 
and we are only interested in the Mj up to conjugacy 

Mi^gMig- 1 . (6.19) 

A simple parameter count shows that the moduli space M. of flat bundles over the s- 
punctured sphere that obey these conditions has complex dimension 2(s — 3). Instead, 
let us ask about the subspace V C M. consisting of flat bundles that can be realized 
by a holomorphic differential equation (6.2). We already know that the potential W 
that appears in the holomorphic differential equation is unique for s = 3. When one 
increases s by 1, adding a new singularity at z = Z{ for some i, one adds another 
double pole to W, giving a new contribution AW — c/(z — Zi) 2 + c'/ (z — Zj). But c is 
determined to get the right monodromy near z±, so only d is a new parameter (usually 
called the accessory parameter). Hence the dimension of V is s — 3; V is a middle- 
dimensional subspace of M.. (For a more complete account of this standard result, see 
section 8 of [21].) Similarly, the subspace V of flat bundles that can be realized by an 
antiholomorphic differential equation is middle-dimensional. 

A complex solution of Liouville theory in the sense that we consider in the present 
section corresponds to an intersection point of V and V. As V and V are both middle- 
dimensional, it is plausible that their intersection generically consists of finitely many 
points, or possibly even that it always consists of just one point. Unfortunately we do 
not know if this is the case. All we really know is that for any s, if the ^ are real and 
obey the Seiberg and Gauss-Bonnet bounds, then there is a real solution of Liouville's 
equations, and this corresponds to an intersection point of V and V. 
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6.3 Interpretation In Chern-Simons Theory 

To explain what one would mean in this language by the action of a classical solu- 
tion, and why it is multivalued, the main idea is to relate Liouville theory on a Riemann 
surface £ to Chern-Simons theory on E x I. The basic reason that there is such a rela- 
tion is that Virasoro conformal blocks (which can be understood as building blocks of 
Liouville theory) can be viewed as physical states in three-dimensional Chern-Simons 
theory. This was first argued in [18] and has been reconsidered much more recently 
[21]- 

We start with an SL(2, C) connection A with Chern-Simons action 

Scs = 4^ J A Tl ( AAdA+ t AAAAA ) ( 6 ' 2 °) 

on a three-manifold M. Scs is invariant under gauge transformations that are con- 
tinuously connected to the identity, but not under homotopically non-trivial gauge 
transformations. For M = S 3 , the homotopically non-trivial gauge transformations are 
parametrized by tts(SL(2, C)) = Z. The integer invariant of a gauge transformation is 
often called winding number. In defining S'cs, we have picked a convenient normaliza- 
tion, such that under a homotopically nontrivial gauge transformation on S 3 of winding 
number n, Scs transforms by 

Scs -> S cs + -p-. (6.21) 

With b understood as the Liouville coupling parameter, this matches the multival- 
uedness of Liouville theory that comes from the trivial symmetry <p c — » <p c + 2i\i. 
Conventionally, the Chern-Simons action is normalized to make S'cs singlevalued mod 
2-7riZ and the homotopically nontrivial gauge transformations are regarded as symme- 
tries [50]. For our purposes, this would be far too restrictive (since we do not want 
to assume that b 2 is the inverse of an integer). Rather, in the path integral, we con- 
sider integration cycles that are not invariant under homotopically non-trivial gauge 
transformations, and we do not view homotopically nontrivial gauge transformations 
as symmetries of the theory. In other words, we adopt the perspective of [11]. The inte- 
gration cycles are middle-dimensional in the space of SL(2, C)-valued flat connections. 
A basis of the possible integration cycles is given by the cycles that arise by steepest 
descent from a critical point of the action. 

The Yang-Mills field strength is defined, as usual, by J 7 = dA+AAA. The classical 
equations of motion of Chern-Simons theory are simply 

7 = 0. (6.22) 
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We take the three-manifold M to be simply M = £ x I, where / is a unit interval and £ 
is the Riemann surface on which we want to do Liouville theory. The fundamental group 
of M is therefore the same as that of E, and so a solution of the classical equation of 
motion (6.22) is just an SL(2, C) flat connection on £. This being so, one may wonder 
what we have gained by introducing a third dimension. 

The answer to this question is that to do Liouville theory, we need more than an 
SX(2, C) flat connection on E. It must obey two conditions: (i) it can be described 
by a holomorphic differential equation, and (ii) it can also be described by an anti- 
holomorphic differential equation. It is possible to pick boundary conditions at the two 
ends of I so that condition (i) is imposed at one end and condition (ii) at the other 
end. Such boundary conditions were introduced in Chern-Simons theory in [18] and 
used to relate that theory to Virasoro conformal blocks; a variant related to Nahm's 
equations and other topics in mathematical physics has been described in [21]. 

We may call these oper or Nahm pole boundary conditions. For the very schematic 
purposes of the present paper, almost all that we really need to know about them is 
that they completely break the SL(2, C) gauge symmetry down to the center ±1 of the 
gauge group. 

Now let us consider the topological classification of gauge transformations on E x I. 
The gauge transformation is described by a map g from £ x / to SL(2, C). At the 
left end of E X I, g must equal 1 or —1. For the present paper, an overall gauge 
transformation by the center of SL(2, C) will be of no interest, since it acts trivially 
on all gauge fields, so we can assume that at the left end, g — 1. But then there are 
two choices at the right end, namely g = 1 and g = —1. After we make this choice, 
the remaining freedom in describing g topologically is given by tt 3 (SL(2,C)) = Z. 
The homotopy classification of gauge transformations is by Z x Z 2 , where the Z 2 factor 
classifies the relative value of g at the two ends and Z classifies the twist by ii 3 (SL(2, C)). 

This last statement is not completely trivial; we must verify that the homotopy 
classification is by a simple product Z x Z 2 and not by a nontrivial extension — > 
Z — 7- T — > Z 2 — > 0. Concretely, the question is whether a gauge transformation with 
g = — 1 on the right end has integer or half- integer winding number. In fact, the 
winding number is always an integer. To see this, it suffices to exhibit an example of a 
gauge transformation with g = — 1 on the right end and with integer winding number. 
We can simply pick g to be a function on £ x / that only depends on the second factor; 
such a map can be constructed from a path from 1 to —1 in the group SL(2, C). (The 
conclusion just stated remains valid when monodromy defects are included, as we do 
momentarily. For this, it suffices to note that by continuity the question is independent 
of the values of the rji, while if one of the rji vanishes, we can forget the corresponding 
defect.) 
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Figure 2. Monodromy defects, depicted as horizontal dotted lines, in E x /. In the example 
shown, S = S 2 and the number of monodromy defects is 3. 

6.4 Liouville Primary Fields and Monodromy Defects 

Some basic points about how to interpret various Liouville fields in the Chern-Simons 
description have been explained in [21]. The main point of concern to us is how to 
incorporate a primary field of Liouville momentum a = r]/b at a point p G £. The 
answer is simply that in the Chern-Simons description, the gauge field A should have 
an appropriate monodromy around the codimension two locus px/cExJ. This can 
be achieved by requiring A to have a suitable singularity along £. If z = re %e is a local 
coordinate that vanishes at p, then we require that the singular behavior of A should 



where the ellipsis refers to less singular terms. This singular behavior has been chosen 
so that the eigenvalues of the monodromy are exp(±27iir]) . We call a singularity of 
this kind in the gauge field a monodromy defect. (What happens when a monodromy 
defect meets the Nahm pole singularity at the ends of / has been analyzed in section 
3.6 of [51]. The details are not important here.) If in the Liouville description, there 
are several primary fields, say at points pt given by z = z iy then in the Chern-Simons 
description we include monodromy defects on pi x I for each i (Fig. 2). 

Now we want to work out the topological classification of gauge transformations in 
the presence of monodromy defects. The ansatz (6.23) is only invariant under diagonal 
gauge transformations along p. So now a gauge transformation g is constrained as 
follows: at the ends of £ x J, it equals ±1 (and we assume it to equal +1 on the left 
end), while along the monodromy defects it is diagonal, and away from the boundary 
and the monodromy defects, it takes arbitrary values in SL(2,C). 

Let us look at what happens along a particular monodromy defect. A diagonal 



be 




(6.23) 
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gauge transformation can be written 



9 = 




e i§ 
p- x e 




(6.24) 



where p is positive and d is real. Let us parametrize the interval / by a parameter y 
that equals, say, at the left end of I and 1 at the right end. We constrain d to vanish 
at y — (where g — 1), and to equal nm at y — 1, where m is even if g = 1 on the 
right end of £ x /, and m is odd if g — — 1 there. We call m the winding number along 
the defect. It is sometimes useful to write it as 



We have normalized m so that wrapping once around the maximal torus of SO (3) corre- 
sponds to m = 1, while wrapping once around the maximal torus of SU{2) corresponds 
to m — 2. 

So in total, with s monodromy defects, the topological classification of gauge trans- 
formations is by integers (n, mi, m?,-, ■ ■ ■ > n^s), where n is the bulk winding number and 
rrii, i = 1, . . . , s are winding numbers defined just along the monodromy defects. The 
rrii are either all even or all odd, since their oddness or evenness is determined by the 
behavior of g. 

In the presence of a monodromy defect, it is necessary to add one more term to 
the action. The reason is that in the presence of a monodromy defect, we would like a 
flat bundle on the complement of the monodromy defect that has the singularity (6.23) 
along the defect to be a classical solution. Such a flat bundle obeys 



Here we write K for the support of the monodromy defect, and 5k is a two-form delta 
function supported on K. But in the presence of the Chern-Simons action (6.20) only, 
the equation of motion is simply J 7 = rather than (6.26). To get the equation we 
want, we must add to the action a term 




(6.25) 




(6.26) 




(6.27) 



Finally, we can determine how the action transforms under a gauge transformation 
in the presence of a monodromy defect. (See section 4.2.6 of [11] for an alternative 
explanation.) We have already discussed the behavior of the Chern-Simons term under 



gauge transformation, so what remains is to understand what happens to the new 
interaction Sk- For K = p x I, only A y , the component of A in the y direction, 
appears in (6.26). Under a diagonal gauge transformation (6.24), the diagonal matrix 
elements of A y are shifted by ^fd\og(pe^)/dy. Taking the trace and integrating over 
y, we find that Sk is shifted by 



(There is no contribution involving <ilog p, since J Q dy{d\og p/dy) = 0, as p = 1 at both 
endpoints.) 

More generally, let us go back to the case of s heavy operator insertions, with 
Liouville parameters rji, i — 1, . . . , s, inserted at points Pi e S. In Chern-Simons theory, 
they correspond to monodromy defects, supported on Ki = p^ x /. In classifying gauge 
transformations, we introduce a winding number associated to each monodromy 
defect. There is also a bulk winding number n. The shift in the total action S = 
$cs + Si under a gauge transformation is 



6.5 Interpretation 

The moral of the story is that in the Chern-Simons description, critical points are 
flat connections on E x /, with prescribed behavior near the ends and near monodromy 
defects, modulo topologically trivial gauge transformations. Topologically nontrivial 
gauge transformations cannot be regarded as symmetries because they do not leave the 
action invariant. Instead, they generate new critical points from old ones. 

For the main example of this paper - three heavy operators on S 2 - all critical 
points are related to each other by topologically nontrivial gauge transformations. This 
means that there is a simple way to compare the path integrals over cycles associated 
to different critical points. 

In fact, let A be a connection that represents a critical point p. Suppose a gauge 
transformation g with winding numbers n and m 1; . . . , m s acts on A to produce a new 
critical point A'. Let Z p and Z p > be the path integrals over integration cycles associated 
to A and to A', respectively. Z p i and Z p are not equal, since the gauge transformation 
g does not preserve the action. But since g transforms the action by a simple additive 
c-number (6.29), there is a simple exact formula that expresses the relation between 
Z p i and Z p : 






(6.29) 




(6.30) 
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We will discuss the interpretation of this formula in Liouville theory in section 7.3 
below. 



7 Timelike Liouville Theory 

So far in this paper we have analytically continued in a but not in b. From the point 
of view advocated in the introduction, this is somewhat artificial; we should allow our- 
selves to consider the path integral with arbitrary complex values of all parameters and 
then study which integration cycles to use to reproduce the analytic continuation from 
the physical region. For complex b's with positive real part, we can simply continue 
the DOZZ formula and the machinery of the preceding sections is essentially unmodi- 
fied. Indeed numerical results for complex b were given in [7], confirming the crossing 
symmetry of the four-point function based on the DOZZ formula. As mentioned in 
the introduction, in various cosmological settings it is desireable to define a version of 
Liouville that has real central charge that is large and negative. The most obvious way 
to try to do this is to continue the DOZZ formula all the way to purely imaginary b, 
since the formula (2.5) will then be in the desired range [22]. This has been shown to 
fail rather dramatically [23], as we will discuss in the following subsection. However, 
we first introduce some conventional redefinitions to simplify future formulas when b is 
imaginary. We begin with 

b=-ib (7.1) 
<f> = i$ (7.2) 

Q = iQ-6j=iQ, (7.3) 
after which the action (2.1) becomes 

Sl = ^LJ [- d ^ d b$9 ab ~ QR$+ 47r//e 2 ^] . (7.4) 

The theory with this action is conventionally referred to as "timelike" Liouville theory, 
since the kinetic term has the wrong sign. In this equation, superscripts are procreating 
at an alarming rate, so we pause to remind the reader that g is the reference metric 
and does not undergo analytic continuation. We will use "hat" , as in b, exclusively to 
refer to the timelike analogues of standard Liouville quantities. The central charge is 
now 

c L = l- 6Q 2 , (7.5) 



- 71 - 



which for small real b accomplishes our goal of large negative central charge. The 
physical metric becomes g a b = e Q 9ab, so the boundary condition on <fi at infinity is 

$(z,z) = -2Q log \z\ + 0{l). (7.6) 

To talk about exponential operators, it is convenient to make one final definition 

a = ia, (7.7) 

which gives conformal weights 

A (e~ 2 ^ = A (e" 2S? ) = a(a - Q). (7.8) 

In the presence of heavy operators 5?j = rji/b, the generalized action (2.25) for the 
rescaled field C = 2b(p = 260 with a flat reference metric is 




Here A = TTfib 2 = —A, and in fact other than an overall sign change this is the only 
difference from the expression of this action in terms of the "unhatted" variables. Note 
that (p c and rji do not need to be "hatted" since they are the same before and after the 
analytic continuation. The equation of motion is now 

dd<f) c = -2A> - 27T^r^ 2 (£ - (7.10) 

i 

which for positive /i is just the equation of motion for constant positive curvature with 
conical deficits at the heavy operators. When b and rji are real and rji is in Region II, 
described by (4.19), this equation has a real solution. As discussed below (4.19), this 
solution can be constructed from spherical triangles. In the FRW/CFT application of 
timelike Liouville, this real saddle point is identified with the asymptotic metric in a 
Coleman-de Luccia bubble [3, 4]. 

7.1 The Timelike DOZZ Formula 

The redefinitions of the previous section make clear that at the classical level the 
relationship between spacelike and timelike Liouville is straightforward. Much less 
clear is the question of the appropriate integration cycle for the path integral when b 
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is imaginary. One way to attempt to specify a cycle is to try to continue the DOZZ 
formula from real b. As just mentioned, this does not work. We can see why by 
considering more carefully the analytic properties of Tb(x) in b [23]. From (A.l) we see 
that the defining integral for Tj,(x) does not converge for any x when b is imaginary, 
which is already a sign of trouble, but this could possibly be avoided by deforming the 
contour. A more sophisticated argument from [23] is as follows: consider the function 



where for the moment we take b to have positive real part and negative imaginary part 
to ensure that both T's can be defined by the integral (A.l). This function is entire 
and has simple zeros everywhere on the lattice generated by b and 1/6, as illustrated 
in Figure 3. Using the recursion relations (A. 3) we can show that H b obeys: 



H b (x) = T b (x)T ib (-ix + ib), 



(7.11) 



H b (x + b)= e^ 2bx -^H b (x) 
H b (x + 1/b) = e^-^H h {x) 



(7.12) 



It is convenient to here introduce a Jacobi ^-function 



10 



Mb 



Figure 3. Zeros of H b (x). The solid circles come from the zeros of T b (x) while the empty 
circles come from zeros of Ti b (—ix + ib). 



9 1 {z,t) =i J2 (-l)V wr(n - 1/2)2+2mz(n - 1/2) lmr>0, (7.13) 

n=— oo 

which is entire in z for any Imr > and obeys 

6 1 (z + l,T) = e-^e 1 (z,r) 

B x { z + r, r) = e^-^O^z, r). (7.14) 

By cancelling the terms n — 1, 2, ... with n, = 0, —1, ... we see that it has a zero at z — 0, 
and by applying these recursion relations we see that it has zeros for all z = m + nr 
with m,n G Z. In fact these zeros are simple and are the only zeros, which follows from 
the standard product representation of the theta function. This function is useful for 
us because we can now observe that 



e 2 



if (x 2 +x/b- 



xh) e x {xibMb 2 ) (7.15) 



obeys the same recursion relations and has the same zeros as Hb(x). Their ratio is 
doubly periodic and entire in x, and must therefore be a function only of b. We can 
determine this function by setting x = | + h and recalling that T b (Q/2) = 1. The 
result is that 41 



9 1 (x/b,l/b 2 ^ 



We can now use this formula to study the behaviour of T& near imaginary b; since 
Yib(— ix + ib) = i we move ^ U P towards the postive real axis then will ap- 
proach the region of interest. But (7.16) reveals that doing this continuation requires 
Q\ to approach the real r-axis. This is actually a natural boundary of analytic contin- 
uation for 9i, with a nonlocal and extremely violent singularity running all along the 
real r-axis. The detailed form of the approach to the singularity depends strongly on z, 
so there is no possibility of cancellation between the two 8is in Hb except for at special 
values of z. This shows that for generic values of x, Tb simply cannot be continued 
to generic imaginary b [23]. This is the origin of the failure of [22] to make sense of 
timelike Liouville theory in this way. 

What then are we to do? One possibility is to restrict to special values of b and 
a where the continuation can still be nonsingular; this is explored in [27]. We are 
interested however in generic complex values of the parameters so this will not work 
for us. A very interesting proposal was made by Al. B. Zamolodchikov [23], and also 



40 Conventions for Jacobi ^-functions are rather inconsistent, so we note that this definition is im- 
plemented in Mathematica as EllipticTheta[l, nz, e l7rr ] . 



41 



This formula was also derived in [27] 
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independently by Kostov and Petkova [24-26]. A key observation is that although 
we cannot continue the DOZZ formula to imaginary b, we can continue Teschner's 
recursion relations. For real b the essentially unique solution of these recursion relations 
is the DOZZ formula, but for generic complex b the solution is not unique since we 
can engineer a-dependent combinations of ^-functions by which we can multiply any 
solution of the recursion relations to produce a new solution. But when we get to 
imaginary b, it turns out that there is again an (almost) unique solution, which is 
not given by analytic continuation of the DOZZ formula. 42 This solution is not quite 
unique because one can multiply it by an S?-independent arbitrary function of b without 
affecting the recursion relations. Fixing this normalization in a way that we explain 
momentarily, and which slightly differs from the choice in [23], the solution is: 



e -^(z2i a i-Q)/ b 



d(a ll a 2l a 3 ) = S \-n f i 1 (-P)b 2+2p 
b L 

T^(5?i + a 2 -a 3 + &)T^(«i + a 3 - a 2 + b)T^(a 2 + a 3 - 5?i + &)T^(Si + a 2 + a 3 - Q + b) 

T ? (6)Tr~(2ai + b)T~ b (2a 2 + 6)T ? (2a 3 +% 

(7.17) 

We will refer to this as the timelike DOZZ formula. The power of 7r/i7(— b 2 ) differs 
slightly from CIO in [23], but the choice we have made here is given by a scaling ar- 
gument in the path integral and is required for our interpretation of timelike Liouville 
as being a different integration cycle of ordinary Liouville. We have also divided C.10 
from [23] by a factor of §-7(1 — b 2 )^{2 — l/b 2 ); as just mentioned these choices do not 
affect the recursion relations and can be interpreted as an ambiguity in the normaliza- 
tion of the operators, but we will see below in section 7.2.3 that the choice we have 
made here is supported by semiclassical computation. Moreover in section 7.3 we will 
see exactly that it is the natural choice for our interpretation of the timelike Liouville 
path integral. 

We can also write down an exact two-point function. Since unlike the three-point 
function the two-point function (2.38) does have a good analytic continuation to imag- 
inary b, it is natural to choose the timelike 2-point function to agree with this analytic 



42 The reason that purely real and purely imaginary b have essentially unique solutions of the recur- 
sion relations is that the lattice generated by b and l/b becomes degenerate in these two cases and 
functions with two real periodicities are highly constrained. The freedom involving multiplying by 
^-functions of the a's goes away since in these cases some of the 9 functions are always evaluated at 
t = 0. The original DOZZ formula (2.28) and the formula (7.17) below are related at complex b by 
such a factor, which is why they can not be continued into each other. For more details see [23, 33]. 
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continuation. This gives 



G(a) = -—C(0,a,a) 

Z7T 



1 

6 2 L 



-It 



(2a-Q)/b 



-iT(2S-Q)/6 T (2S6+6 2 ) 7 



1 

fe 2 



2a 



6 

(7.18) 

Note that for real a, this expression is not positive-definite, as expected from the 
wrong-sign kinetic term. Its relation to the three-point function is somewhat arbitrary, 
unlike in spacelike Liouville where there was a clear rationale for the formula (2.37). 
In particular setting one of the a's to zero in the timelike DOZZ formula does NOT 
produce a 5-function. Indeed the timelike DOZZ formula has a finite and nonzero 
limit even when Si — > 0, «2 ^ S3. In [23], this was observed as part of a larger 
issue whereby the degenerate fusion rules mentioned below equation (2.43) are not 
automatically satisfied by the timelike DOZZ formula. In [28], this was interpreted 
as the two-point function being genuinely non-diagonal in the operator weights. We 
will not be able to explain this in a completely satisfactory way, but we will suggest a 
possible resolution below in section 7.4. 

7.2 Semiclassical Tests of the Timelike DOZZ formula 

In this section we will show in three different cases, analogous to the three cases 
studied above for the spacelike DOZZ formula, that the semiclassical limits of (7.17) and 
(7.18) are consistent with our claim that they are produced by the usual Liouville path 
integral on a different integration cycle. This task is greatly simplified by observing 
that we can trivially reuse all of our old solutions (or "solutions") and expressions 
for the action. Say that we have a solution Cj jv(?7i, X,b, z,z) of the original Liouville 
equation of motion (2.15). It is easy to check that 



<Pc,n{vu A, b, z, z) = c ,at(^, A, b, z, z) - in 



(7.19) 



obeys (7.10). We can also compute the action by noting that if we define the original 
modified action (2.25) as Sl [rji, A, b, Zi, z$\ 4> c ,N{j)i, A)], then we have 



Si 



r)i, -A, -ib, z u z^, (/)c,N(Vh ^) 



S, 



r]i, \,b, z h z^ (j) C)N (r]i, X) 



T)i, X,b, Z h Zi, (f> c ,N{Vh ^) 



b 2 



(7.20) 



The left hand side of this is just (7.9), so we can thus compute the action for timelike 
Liouville theory by simple modification of our previous results. 
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7.2.1 Two-point Function 

Using (7.20) and (4.8), we find that the timelike version of the saddlepoint (4.3) has 
timelike action 



S 7 



1 

b 2 



2mN(l — 277) + (277 — 1)A + 2 ( (1 — 277) log(l - 2rj) - (1 - 2t?) 



+ 2(77 — ?7 2 ) log |^i 2 | 



(7.21) 



The semiclassical limit of (7.18) with a heavy is: 



G( v ) (e 2 ^ 1 - 2 ^ - l) 



ljexpji [-(l-2rj)logA + 



2((1 - 2ri)log(l -2rj) -(l-2r7)) 

(7.22) 



which is matched by a sum over the two saddle points N = and N — 1 with actions 
given by (7.21). Note that the integral over the moduli would again produce a diver- 
gence, but that unlike in the DOZZ case this divergence did not seem to be produced 
by the limit a\ — » 0, a 2 = «3- Note also that there is now no Stokes line in the 77 
plane, there are always only two saddle points that contribute. This is analogous to 
the integral representation of 1/T(z) as discussed in appendix C. 

7.2.2 Three-point Function with Heavy Operators 

Similarly for three heavy operators in Region I, the timelike version of (4.21) has 
timelike action 



Sr. = ~; 



b 2 



1 - y^Vi) lo S^ - @i + ^2 - 63) log \z 12 \ 2 - (5 X + 5 3 - 8 2 ) log I ^13 1 2 

(8 2 + 83- 81) log \z 2 3\ 2 + F(t]i + 772 - Vs) + F(Vi +V3- m) + F{r]2 + V3- Vi) 
F(r]i + r]2 + 773 - 1) - F(2 Vl ) - F{2r l2 ) - F(2 m ) - F(0) 



2m(N- 1/2)(1 



(7.23) 



where we have defined 8{ = r)i(r)i — 1), consistent with (7.8). This clearly has the right 
position dependence for a timelike three-point function. The action in other regions is 
always an analytic continuation of this action along some path, but to be definite we 
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also give the timelike action from Region II as well: 



S, 



1 

b 2 



1 - lo S^ - (Si + 8 2 - 8 3 ) log \z 12 \ 2 - (?i + ? 3 - ? 2 ) log ki 3 | 2 

(8 2 + 83- 81) log I223 1 2 + F(r]i +rj 2 - r] 3 ) + Ffa + r] 3 - r/ 2 ) + F(?7 2 +r] 3 - 771) 
Ffa + % + Vs) - F(2 Vl ) - F(2 m ) - F{2r l3 ) - F(0) 



2 \ (1 - j2 vt) iog(i -E^-(!-E^) + 2?riiV ( 1 - E 



(7.24) 



To compare these with the timelike DOZZ formula, we can again make use of the 
asymptotic formula (4.34). The terms in (7.17) that don't involve approach 

e l=Soi (i7r+ 2iog6-iogA)+o(i/6) j / 725 n 

and using (4.34) we find that in Region I the T^'s combine with this to give 
1 r 



C(Vi/b) ~ exp 



b 2 



- f i - E^l log ^ + + 712 ~ % ) + F ^ ni + ^ - 

+ i 7 ^ +V3- Vi) + F(vi +V2 + V3-I)- F(2r ]l ) - F(2r] 2 ) 



F(2 V3 )-F(0)+iir(l-J2vi) 



(7.26) 



Comparing with (7.23), we see that only the saddle point with N = contributes. In 
Region II as before (4.37) we need to use (A. 4) to shift one of the Tg's before using the 
asymptotic formula (4.34), giving: 



Ti 



2 &) ~ 7 ((^^-l)/6 2 ) _1 6p( 2(1 -^^- ( ^^- 1/2)2 )eP F( ^^. (7.27) 



The result in Region II is 
1 



C(r)i/b) ~ exp 



b 2 



- f 1 - E^l logA + F(r]i + r? 2 - 773) + ^(?7i + V3- r} 2 ) 

+ F{ V2 +tj 3 - Vl ) + F{ Vl + V2 + 773 - 1) - F{2i ll ) - F(2r7 2 ) 

- F(2 % ) - F(0) +2(1-J2 Vi) log(l - - 2 ^ - E ^) 



+ Z7r(l - J^T?; 



^(i-EiW)/^ 2 _ e -i»r(i-Ei •»)/&" 



(7.28) 
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Comparing this with (7.24), we see it matches a sum over two saddle points with N = 
and N — 1. Unlike the spacelike DOZZ formula there are no Stokes walls in Region II, 
in complete analogy with the situation for 1/T(z) explained in appendix C. 

7.2.3 Three-point Function with Light Operators 

As a final semiclassical check of the timelike DOZZ formula (7.17), we will calculate 
its b limit when all three operators are light and compare to a semiclassical 
computation analogous to that from section 4.3. We will define Oi = y = — y, which 

gives A — > a as b — > 0. Manipulations similar to those leading up to (4.44) now give 
C(-a$, -<t$, - a 3 b) = _ 27 rz6- 3 A 1 -^^- 1 /6 2 e -2/6 2 -27 £ +o(biog6) ^(E^-m/^) _ 1 

r(i - 2<7i)r(i - 2<7 2 )r(i - 2<r 3 ) 



X 



T(l +<ti-<t 2 - o"3)r(l + a 2 -a 1 - (T 3 )r(l + a 3 - a x - cr 2 )r(2 - £\ a { ) ' 

(7.29) 

From the structure of this formula, it appears that we will be able to interpret as a 
sum over two complex saddle points as with Region II in the previous section. There 
is a subtlety however in that to produce the T-functions that will emerge from the 
modular integral in our imminent semiclassical computation, we need to apply the 
Euler reflection formula T(x)T(l — x) = ir/ sin7rx to each of them. Anticipating this 
result, we write: 

di-aj, - a 2 b, -a 3 b) = 6-3^1-^,^-1/6^-2/^-2^+0(6^6) ^ME^-i+i/S") _ 
T(a 1 + a 2 - cr 3 )T(ai + o- 3 - ^2)^(^2 + ^3 ~ gjjTQ^ Qj - 1) 

r(2a 1 )r(2 ( 7 2 )r(2 ( 73) 

X ^g47ri(Ti — 1) ( e 4vrio- 2 — 1) (e 47 ™°"3 — 1) 

(7.30) 

The structure of the terms in the third line show that a much more complicated set of 
saddlepoints are needed to explain this result than in the spacelike case (4.44). At the 
end of this section we will explain why this happens. 

The semiclassical formula analogous to (4.46) for this correlation function is 



r 3 

NGT J i=l 

(7.31) 

Here we have assumed that the fluctuation determinant and Jacobian parametrized by 
A(—ib) are just the analytic continuations of their spacelike counterparts. This should 



- 79 - 



be true if our path integral interpretation is correct, and we will see momentarily that 
this works out. <p c ^ is the timelike "solution" with branch choice N, related to the 
usual spacelike "solution" by (7.19). Explicitly 

$ c (z,z) = 2iriN(z,z) - log A - 21og(|az + /3| 2 + \-?z + 5\ 2 ). (7.32) 

Based on (7.30), we have allowed iV to vary with position to allow the different branches 
of the action to be realized. This is one of the situations discussed in section 5.3.1 where 
discontinuous "solutions" must be included even though single- valued solutions exist. 
Computing the action (7.4) and simplifying the modular integral as in section 4.3 we 
find 

d{-a x % -a 2 b, -a 3 b) « AH^-^-^V^/K a 2 , a 3 ) ^ e -^i^ i+ n/P) 

7VeT 

(7.33) 

Here — n is the value of N at oo and — is its value near the various insertions. 
Using (4.51) and comparing this with (7.30), we find complete agreement provided 
that A(—ib) = 5~ 3 7r~ 3 e~ 27E . Recalling that at the end of section 4.3 we found A(b) = 
ib~ 3 7T~ 3 e~ 2,yE , this indeed works out as expected. 

The set T of included branches is now rather complex; it can be read off from 
(7.30) but we will not try to characterize it more precisely. We observe however that 
many branches that correspond to discontinuous "solutions" are now definitely needed. 
This is different than what we found for spacelike Liouville in section 4.3, where the 
contributing saddle points were single-valued and continuous and just the same as in 
Region II for heavy operators. The reason for this distinction is that, as explained 
in appendix F, the modular integral over SL(2, C) converges only when the ex's obey 
certain inequalities (F.20). In spacelike Liouville with the ex's in Region II the integral 
is convergent, so we can evaluate it without any contour deformation. In timelike 
Liouville, when the cr's are in Region II many of the inequalities are violated and 
the integral must be defined by analytic continuation. This continuation results in 
additional Stokes phenomena, which changes the contributing saddle-points. 

7.3 An Exact Check 

The checks of the previous section were semiclassical, but we will now give an exact 
argument that the timelike DOZZ formula (7.17) is produced by evaluating the usual 
Liouville path integral on a new integration cycle. We will show that the ratio of the 
spacelike and timelike DOZZ formulae must have a specific form and then demonstrate 
that it does. 
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We begin by defining: 

Z p (oci,Zi,Zi) =/ V<p c V ai {z 1 ,z 1 )...V an {z n ,z n )e~ SL , (7.34) 

JC P 

where here p is a critical point of the action with heavy operators as sources, and the 
path integral is evaluated on the steepest descent cycle C p that passes through p. As 
discussed in the introduction, this quantity is not in general equal to the Liouville 
correlator; we need to sum over such cycles with integer coefficients a p as in (1.6). We 
will now argue however that the p-dependence of Z p is quite simple. First recall the 
exact version of the action (2.10) 43 

Sl = ^2 / ^ H 9 ^) 2 + 16Ae0c ] + + &2 ) / ^ d6 +^( 1 + 2b2 + &4 ) lo g^ 

(7.35) 

We note that under the transformation <p c — > <p c +2niN, we have Sl — > Sl + ^^-(1+o 2 ). 
Semiclassically the operator V a defined by (2.24) transforms as V a — >■ V a e 2ma ^ b under 
the same transformation. Since the Seiberg bound ensures that the renormalization 
needed to define this operator precisely is the same as in free field theory, this is 
actually the exact transformation of V a . Moreover the path-integral measure T><p c is 
invariant under the shift. This means that if two p's differ only by adding 2iriN, then 
with a slight abuse of notation we have the simple relation 

Z p+2mN = e^ N ^' b '^Z p . (7.36) 

This result is exact, and more generally it shows that the result of integrating over a 
sum of integration cycles of this type can be factored out from the correlator: 

oo oo 

Z= aP +2mN Z p+2mN = Z p a^ 2mN e 2mN ^^l b ~ l ' b2 \ (7.37) 

N=-oo N=-oo 

Thus in general the ratio of Z's which are computed on different cycles, both of the 
form Y^n=-oo aP+2mN Cp+2TriN, will be expressible as a ratio of Laurent expansions in 
e 2TTi(J2iOi t /b~i/b 2 ) integer coefficients. This is a rather nontrivial constraint; for 

example it implies that the ratio is invariant under shifting any particular aij by — > 
on + b. There is also a more subtle invariance of the form b — >■ ^=p and ctj — > ^=p- 
Unfortunately as discussed in section 5, to understand the DOZZ formula in the 
full range of ctj's it is not sufficient to only consider integration cycles that differ by 



43 When we discuss discontinuous "solutions" momentarily, the kinetic term should really be under- 
stood to be expressed in terms of / in equation (5.3). 
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a global addition of 2iriN. We found semiclassically in (5.2) that to fully explain the 
DOZZ formula it was necessary to consider discontinuous "solutions" that differ by 
different multiples of 2ni at the different operator insertions. To proceed further we 
need to assume that we can apply the machinery of the previous paragraph to these 
"solutions" and their associated "integration cycles of steepest descent." The idea 
is that the action (7.35), with the kinetic term expressed in terms of / as in (5.3), 
changes only by an overall c-number if we shift the field configuration by 2niN with 
a position-dependent N G Z. The change in the action depends on the value of iV at 
infinity, and the contributions of operator insertions also shift in a way that depends 
on the value of N in their vicinity. For the reader who is uncomfortable with this, 
we note that in the Chern-Simons interpetation espoused in section 6, these additional 
"solutions" were just as valid and conventional as the usual ones. So one could in 
principle rephrase what follows in Chern-Simons language, which would perhaps make 
it sound more plausible. We will henceforth assume that the relationship between Z p 
and Z p i is given by a formula analogous to (6.30) in the Chern-Simons version: 

Z p , = Z p e-^( n+ ^ miaib ) (7.38) 

Here n and rrii are the differences in N at infinity and near the various operator inser- 
tions, and rrii are either all even or all odd. 

We will now compute the ratio of the spacelike and timelike DOZZ formulas for 
a region of b where both make sense, with the goal being to check that their ratio is 
consistent with this result. Using (7.17) expressed in terms of the "unhatted" variables 
as well as (2.28), (7.16), and (7.11), we find: 

C(-iai,-ta2,-ia 3 ) = _ 2 Il lim T &( e ) ^g/b-bMQ-T, m) 
C(ai,a 2 ,a 3 ) b ^oT H b (e) 

Hb(J2j oii - Q)H b {ai + a 2 — otz)H b {ai + a 3 - a 2 )H b (a 2 + a 3 - «i) 
X H b {2a 1 )H b (2a 2 )H b {2a 3 ) 
= -27ue-W(^Mi+fc 2 )/2) 

fligjFz, ^giggg^ ^Wii^f^, ^fliQgff^ j) , 

ei(O,^)0 1 (2a 1 /6,^)0 2 (2a 1 /6,^)0 3 (2ai/6,^) ' 

Here 9[(z,t) = ^(z,t). We can simplify this a bit by using (7.14) to shift the 
argument of one of the ^-functions: 

C(-ia 1 , -ia 2 , -ia 3 ) _ Vi{— ^— , p)Vi{ b , p)Vi\ b , pJ y n b > 



C{a u a 2 , a 3 ) " #1(0, ^)0 1 (2a 1 /6 J ^)# 2 (2a 1 /6, ±)6 3 {2a x /b, 

(7.40) 
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We'd now like to express this as a ratio of sums of terms of the form e iF\ n+ ^>i miaib ) 
with integer coefficients. To facilitate this, we define 



e 1 {z,T) = -ie- iwT ' A+i ' Kz e l {z ) T)= J2 



^iiTTn(n— l)+2nizn 



n=— oo 
oo 



9 (r) = -^-e- i7TT/4 9' 1 (0,T)= V (-l) n ne*" Tn ^ 1 >, (7.41) 

Z7T * — ' 

n=— oo 

in terms of which we have: 

„'™ „'™ \ ft I Si 1 ^fl f ai +02-0:3 1 Nfl f 01+03-0:2 1 Nfl / 02+03-01 1 \ 



C(a u a 2 , a 3 ) e Q {^)6 l {2a l /b, £)0 2 (2ai/6, &)9 3 (2 ai /b, ±) 

(7.42) 

From (7.41), we see that the right hand side of this equation now explicitly is a ratio of 
the desired form. This completes our demonstration that the timelike DOZZ formula 
is given by the ordinary Liouville path integral evaluated on a different integration 
cycle. Note in particular that the ratio of the two is bad both for purely real and 
purely imaginary b, which illustrates the failure to directly continue between the two. 
This argument also confirms our choice of prefactor in the timelike DOZZ formula, 
since other choices, including the one made in CIO from [23], would have polluted this 
result. 

7.4 Is Timelike Liouville a Conformal Field Theory? 

Unlike most sections which are titled by questions, in this case our answer will be an 
optimistic "maybe". We have established that the timelike DOZZ formula is computed 
by evaluating the Liouville path integral on a particular choice of cycle, which means 
that its correlation functions will necessarily obey the usual conformal Ward identities. 
So in the sense that any local path integral which computes correlators that obey the 
conformal Ward identities is a conformal field theory, it is clear that timelike Liouville 
theory fits the bill. For example as a consequence of this our semiclassical computations 
confirmed the usual position dependence of the two- and three-point functions. But the 
real meat of this question is understanding to what extent timelike Liouville theory fits 
into the standard conformal field theory framework of [17]. At least one thing that seems 
to work is that the derivation of the timelike DOZZ formula from the recursion relations 
confirms that the four-point function with a single degenerate operator constructed in 
the standard way is crossing symmetric. There has however been justifiable concern in 
the literature [22, 23, 28] about the fact that the timelike DOZZ formula does not obey 
the degenerate fusion rules when its arguments are specialized to degenerate values. 
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The simplest manifestation of this is the nonvanishing of C(0, o?i, a 2 ) when Si 7^ a 2 , 
as discussed below (7.18). 

The reason that this is troubling is that semiclassically it seems obvious that 
lim a _K) e 2Q ^ = 1. If this were really true as an operator equation, it would imply that 
the timelike Liouville two-point function is nondiagonal in the operator dimensions. 
Since the diagonal nature of this function is a consequence only of the Ward identities, 
and we know the Ward identities are satisfied just from the path integral, something 
has to give. What the Timelike DOZZ formula seems to tell us is that sending a — > in 
the three-point function does not produce the identity operator, but instead produces 
another operator of weight zero that does not obey the degenerate fusion rule. 44 The 
existence of such an operator is usually forbidden by unitarity, but timelike Liouville 
is necessarily nonunitary so this does not contradict anything sacred. We believe that 
this is the correct interpretation. 45 

As evidence for this proposal, we consider the differential equation obeyed by u, v 
for the semiclassical three-point function with three heavy operators: 

d 2 u + W(z)u = 

with 

771.(1 - 771) 2:122:13 772(1-772)^21^23 . ?73(1 - ?73)^3i^32l I 



W(z) 



Z — Z\ z — Z2 z — 23 



(z - zi)(z - z 2 )(z - z 3 )' 



We observe that if 77! — > 0, there is still a regular singular point at z — Z\ that only 
cancels if we also have 772 = 773. 46 When 772 7^ 773, the solution will generically have a 
logarithmic singularity at z = z\. In this limit the standard semiclassical solution (3.25) 
that we reviewed previously breaks down, and a new solution needs to be constructed. 
We interpret this as the three-point function of a new nontrivial operator of weight zero 
with two conventional Liouville operators. 47 In the spacelike case this also could have 
happened, but since for real b spacelike Liouville is unitary such an operator cannot 



44 Recall that even in spacelike Liouville there was a subtlety with computing degenerate correlators 
by specializing the general correlators to degenerate values, as discussed below (2.43). 

45 We thank V. Petkova for useful correspondence on this point. He points out that this nondecou- 
pling happens in the c < 1 Coulomb gas formalism, which is closely related to the Timelike DOZZ 
formula evaluated at degenerate points. 

46 We are here assuming the Seiberg bound Re(rji) < 1/2. 

47 Thc monodromy matrix M of the differential equation about z\ in this limit has in some basis the 
form ( \\i with A some function of 772 and 773. This matrix has one eigenvector with eigenvalue 1, but 



,A 1, 

is not diagonalizeable. In the Chern-Simons interpretation this means that there is still a monodromy 
defect in the gauge field even after we send r\\ — > 0. 
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exist and the O(b ) corrections to the saddlepoint conspire to set the correlator to zero. 
In timelike Liouville there is no reason for this conspiracy to happen, and indeed from 
the timelike DOZZ formula we see that it does not. We take the fact that this extra 
singularity disappears only when rj2 = rjz as evidence that, contrary to the worries 
expressed in [22, 23, 28], the real two-point function of timelike Liouville theory is 
indeed diagonal. 

Perhaps a natural framework to discuss a CFT that includes an extra operator of 
dimension that does not decouple is "logarithmic" CFT. Something which remains 
mysterious about this interpretation however is that there does not seem to be any 
candidate primary operator we can express in terms of the Liouville field to play this 
role. We leave this unresolved for future work, but we note that in the Chern-Simons 
formulation it is straightforward to describe the nondegenerate primary of dimension 
0; it corresponds to a monodromy defect with unipotent monodromy as explained in 
the footnote on the previous page. 

The main open question that would allow a more systematic understanding of 
timelike Liouville as a CFT is to identify the set of states on which we should factor- 
ize correlation functions. In spacelike Liouville theory this question was answered by 
Seiberg [8], and is formalized by the expression (2.39) for the four-point function. In 
that case the key insight came from study of minisuperspace and the analogy to scat- 
tering off of an exponential potential. A similar analysis for timelike Liouville theory 
was initiated in [29] and studied further in [28], but the Hamiltonian is non-hermitian 
and subtle functional analysis seems to be called for. We have not tried to extend the 
minisuperspace analysis of [28, 29] to the full timelike Liouville theory, but it seems 
that this would be the key missing step in establishing the appropriate basis of states 
to factorize on. This would then allow construction of four-point functions as in (2.39) 
for spacelike Liouville, and one could check numerically if they are crossing symmetric. 
Since in the end of the day we know that the path integral does produce crossing- 
symmetric four-point functions that obey the Ward identities, it seems certain that 
such a construction is possible; it would be good to understand it explicitly. 

8 Conclusion 

In this conclusion we summarize our main results and suggest a few directions for future 
work. We began by trying to assign a path integral interpretation to the full analytic 
continuation of the DOZZ formula (2.28) for the three-point function of Liouville pri- 
mary operators. Our technique was to compare the semiclassical limit of the DOZZ 
formula and various other correlators to the classical actions of complex solutions of 
Liouville's equation. We found that for certain regions of the parameters the analytic 
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continuation is well described by the machinery of Stokes walls and complex saddle 
points, and in particular we showed that the old transition to the "fixed-area" region 
[8] can be reinterpreted in this manifestly-local language. The main surprise was that 
in order to properly account for the full analytic continuation it was necessary to in- 
clude multivalued/discontinuous "solutions", whose actions were defined according to a 
simple prescription in section 5.1. In section 5.5,we saw that these singularities naively 
suggested singularities in the four-point function, but argued that they were in fact 
resolved by quantum corrections. One is tempted to declare this an example of the 
quantum resolution of two-dimensional gravitational singularities. 

Two situations come to mind where these ideas may be relevant. In [52] a statistical 
model of bubble collisions in three-dimensional de Sitter space was constructed in which 
the 4-point function had singularities when the operators were not coincident, similar to 
the naive result in section 5.5. This theory has a yet to be well-understood relationship 
to dS/CFT in three dimensions, which is expected to have a Liouville sector coupled to 
a nonunitary CFT [5] . Perhaps in a more refined version of this theory the singularity 
could be resolved as in section 5.5. Secondly, in three-dimensional Euclidean quantum 
gravity with negative cosmological constant, it was found in [53] that including only real 
smooth solutions in the path integral produces a partition function that does not have 
the correct form to come from a CFT computation. Perhaps other complex "solutions" 
need to be included? 48 More generally the use of singular "solutions" seems to be a 
new phenomenon in field theory and we wonder where else it could appear. 

We then discussed how the question of analytic continuation could be reformulated 
in the Chern-Simons description of Liouville theory, where we found that the picture 
of analytic continuation in terms of Stokes phenomenon is more conventional and all 
relevant solutions seem to be nonsingular. It would be interesting to get a more precise 
picture of these solutions; since explicit formulas exist on the Liouville side it seems 
plausible that they may also be achievable on the Chern-Simons side. This would allow 
a more concrete realization of the ideas suggested in section 6. 

Finally we used the tools developed in the previous sections to discuss an expression 
(7.17), proposed in [23], for an exact three-point function in timelike Liouville theory. 
We found that we could interpret this formula as being the result of performing the usual 
Liouville path integral on a different integration cycle, which we demonstrated both 
semiclassically and exactly. We also discussed the extent to which timelike Liouville 
theory can be understood as a conformal field theory, arguing that it probably can but 
that the spectrum of states to factorize on needs to be understood before the question 

48 For quantum gravity in higher dimensions it is unlikely that the path integral makes sense beyond 
the semiclassical expansion about any particular background, so in particular it is probably not well- 
defined enough for us to ask about Stokes phenomenon. 
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can be decisively settled. Even before this question is addressed however, we already 
consider our results sufficient motivation to begin using the formula of [23] to study 
the various proposed applications of timelike Liouville theory to closed string tachyon 
condensation [22] and FRW/CFT duality [3]. 
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A Properties of the T5 Function 

The function has now become standard in the literature on Liouville theory, but 
for convenience we here sketch derivations of its key properties. The function can be 
defined by 



Here Q = b + \. The definition reveals that Tb(Q — x) = Yf,(x). When x = the 
second term in the integral diverges logarithmically at large t, and at small but finite 
x it behaves like logx. Yj, therefore has a simple zero at x = as well as x = Q. 
To extend the function over the whole x-plane, we can use the identity 49 




sinh 2 ((g/2-x)|) 
sinh ^ sinh ^ 



< Re(x) < Re(Q). 



(A.l) 
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(A.2) 
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This identity is derived in Appendix E. 
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to show that in its range of definition T& obeys 



Where: 



T b (x + 6) = 7 (te)6 1 - 2ba; T b (x) 
T b (x + l/b) = 7 (x/6)6^- 1 T ;) (x). 

, , r(x) 

7(rc) 



(A.3) 



r(i -x) 

These recursion relations are the crucial property of T b from the point of view of 
Liouville theory, among other things they are what allow a solution of Teschner's re- 
cursion relations to be expressed in terms of T&. The recursion relations also show 
that the simple zeros at x = 0, Q induce more simple zeros at x = —mb — n/b and 
x = {m! + 1)6 + (n' + l)/b, with m,m' and n,n' all non-negative integers. It it is also 
useful to record the inverse recursion relations: 



T b (x - b) = 7 (fcr - 6 2 )- 1 b 2bx - l - 2b2 T b (x) 
T b {x - 1/6) = 7 (x/6 - I/ft 2 )" 1 b 1+ ^-^T b (x). 



(A.4) 



We will also need various semiclassical limits of Y^. 50 Rescaling t by b and using 
the identity 



\ogx 



J j [<r* - e- xt ] Re(x) > 0, 



we see that 

b 2 \ogT b {r,/b + b/2) 

+ 



00 dt 

T 



n 



log b 



t 



2U4 



1 - 



~24 



+ 



sinh 2 [(r) - \)t/2] 



sinh 



(A.5) 



When < Re(r/) < 1, the subleading terms in the series 1 + t 2 b A + . . . can be integrated 
term by term, with only the 1 contributing to nonvanishing order in b. From the 
identity (A. 2), we can find 



F(rf) = I \og / y(x)dx 

1/2 



dt 

T 



(V - 1/2)V* - 



2sinh 2 ((r7-l/2)|) 



sinh ( |) 



< Re(r/) < 1, 



so using this we find the asymptotic formula: 

T b ( v /b + 6/2) = e U-^/^o Sb +F( v)+ o^)] 



< Re(^) < 1. 



(A.6) 



51) 



We thank A. Zamolodchikov for suggesting the use of (A. 2) in the first of these derivations. 
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In particular if we choose rj to be constant as b — » only caring about the leading 
terms, then we have 

T b {r}/b) = e ^[^)-(^i/2) 2 iog^(Mo gfe) ] < Re{r]) < 1? (A ?) 

which is useful for our heavy operator calculations in section 4. 

For light operator calculations we will also be interested in the situation where the 
argument of scales like b. Looking at the b — y limit of the first recursion relation 
in (A. 3) we find 

T b ((a + l)b)^\r b (ab). (A.8) 

(TO 

One solution to this relation is 

T b (ab) « ^h(b), (A.9) 

where h(b) is independent of a. Unfortunately this solution is not unique since we can 
multiply it by any periodic function of a with period one and still obey the recursion 
relation. We see however that it already has all of the correct zeros at a = 0, —1, —2, ... 
to match the T b function, so we might expect that this periodic function is a constant. 
This periodic function in any case is nonvanishing and has no poles, so it must be the 
exponential of an entire function. If the entire function is nonconstant then it must 
grow as a — > oo, which seems to be inconsistent with the nice analytic properties of T b . 
In particular (A. 6) shows no sign of such singularities in rj as rj — > 0. We can derive 
h(b) analytically, up to a b- independent constant which we determine numerically. The 
manipulations are sketched momentarily in a footnote, the result is 



^b(vb) = exp 
I cr 



-^log6 + F(0)/6 2 + O(6 2 log6) 



(A.10) 



The numerical agreement of this formula with the asymptotics of the integral (A.l) is 
excellent;in particular we find C = 2. 50663. 51 The constant C will cancel out of all 



51 To do this numerical comparison, it is very convenient to first note that for Re(cr6) > 0, we have 
log T b ((or + \)b) = -(i - obf logfo + ±F(ab 2 ) + I(a, b), with 

2 



I(a,b) 



dt (2 1 \ sinh 2 [(25? - 30*/2] 



t \t sinh(£/2)y sinh 



; 

2b 2 



This integral approaches a finite limit as b — > 0, which makes it easy to extract the leading terms 
in (A.10) and also to do the numerical comparison with (A.l). Although this derivation required 
restrictions on cr, the final result does not since it can be continued throughout the a plane using the 
recursion relation. Of course the asymptotic series is only useful when a is O(b ). 
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of our computations since we are always computing ratios of equal numbers of TVs. 
This precise numerical agreement also confirms our somewhat vague argument for the 
absence of an additional periodic function in a. As an application of this formula we 
can find the asymptotics of Tq from the DOZZ formula: 



C 



~ log b + F(0)/b 2 + 0(b 2 log b) 



(A.11) 



B Theory of Hypergeometric Functions 



This appendix will derive the results we need about hypergeometric and P-functions 
[54]. No prior exposure to either is assumed. Our initial approach is rather pedestrian; 
it is aimed at producing concrete formulas (B.10-B.13) which illustrate the monodromy 
properties of various solutions of Riemann's hypergeometric differential equation and 
the "connection coefficients" relating them. This "toolbox" approach is convenient for 
practical computations, but the disadvantage is that it involves complicated expressions 
that obscure some of the underlying symmetry. In section B.5 we give a more elegant 
general formulation in terms of the integral representation, which illustrates the basic 
logic of the previous sections in a simpler way but is less explicit about the details. It 
also allows us to recast the three-point solutions of section 3.3 in an interesting way. 

B.l Hypergeometric Series 

We begin by studying the series 

F(a,b,c,z) = f2^^z". (B.l) 
(c) n n. 

n=0 

Here (x) n = x(x + l) • • • (x + n — 1) = F p^"- > . It is easy to see using the ratio test that if 
c is not a negative integer, then for any complex a and b the series converges absolutely 
for | z | < 1, diverges for \z\ > 1, and is conditional for \z\ = 1. It is also symmetric in a 
and b, and we can observe that if either a or b is a nonpositive integer then the series 
terminates at some finite n. One special case which is easy to evaluate is 

F( M ,2,,) = -Mlz£), 

which shows that the analytic continuation outside of the unit disk is not necessarily 
singlevalued. 

We will also need the value of the series at z — 1. We will derive this from the 
integral representation in section B.5: 

n«AcA) = *f V{c ~*~ h l Re(c-a-b)>0. (B.2) 
r(c — a)T(c — b) 
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B.2 Hypergeometric Differential Equation 

By direct substitution one can check that the function F(a, b, c, z) obeys the follow- 
ing differential equation: 

z{\ - z)f" + (c-(a + b + l)z)f - abf = 0. (B.3) 

This second-order equation has three regular singular points, at 0,1, and oo. Since 
F(a, b, c, z) is manifestly nonsingular at z = 0, its analytic continuation has potential 
singularities only at 1 and oo. There is the possibility of a branch cut running between 
1 and oo. We saw this in the special case we evaluated above, and it is standard 
to choose this branch cut to lie on the real axis. We can determine the monodromy 
structure of a general solution of this equation by studying its asymptotic behaviour in 
the vicinity of the singular points. By using a power-law ansatz it is easy to see that 
any solution generically takes the form 52 

A Q (z) + z 1 - c B Q {z) as z ->■ 

/(:) - { z-MooCl/z) + z^B^il/z) as z ->• oo , (B.4) 

Ai(l - z) + (1 - z) c - a - b B x (l - z) as 2^1 

with Ai(-),Bi(-) being holomorphic functions in a neighborhood of their argument being 
zero. The solution of (B.3) defined by the series (B.l) is a case of (B.4), with Aq(z) = 
F(a, b, c, z) and Bq = 0. We will determine the A{ and Bi at the other two singular 
points later. These expressions confirm that a solution of (B.3) will generically have 
branch points at 0, 1 and oo. 

B.3 Riemann's Differential Equation 

It will be very convenient for our work on Liouville to make use of Riemann's 
hypergeometric equation, of which (B.3) is a special case. This more general differential 
equation is: 



+ 1 1 -g- g' + 1-/3-/3' + 1 - 7 ~ Y ) ft+ 

\ Z - Z\ z- z 2 z- z 3 J 

aa'z 12 z 13 (3f3'z 21 z 23 Jj'z 3 iZ 3 2 } f 



(B.5) 



Z — Zi z-z 2 z-z 3 J (z - Zi){z - z 2 ){z - z 3 ) 

along with a constraint: 

a + o / + /3 + /3 , + 7 + 7 / = 1. (B.6) 

52 When the two exponents given here become equal at one or more of the singular points, there are 
additional asymptotic solutions involving logs. We will not treat this special case, although it appears 
for the particular choice a = b = 1/2, c = 1 in appendix D. 
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Here z^ = z^ — Zj, the parameters a, (3, 7, a', (3', 7' are complex numbers, and the 
constraint is imposed to make the equation nonsingular at infinity. The points 2, are 
regular singular points. This is in fact the most general second-order linear differential 
equation with three regular singular points and no singularity at infinity. To see that 
this reduces to the hypergeometric equation (B.3), one can set z\ = 0, z 2 = 00, z 3 = 1, 
a — 7 = 0, /3 = a, /3' — b, and a' = 1 — c. 

Solutions to Riemann's equation can always be written in terms of solutions of the 
hypergeometric equation; this is accomplished by first doing an SL(2, C) transformation 
to send the three singular points to 0, 1, and 00, followed by a nontrivial rescaling. To 
see this explicitly, say that g(a,b,c,z) is a solution of the differential equation (B.3), 
not necessarily the solution given by (B.l). Then a somewhat tedious calculation shows 
that 

\z-z 2 ) \z-z 2 ) V z 13 (z-z 2 )J 

is a solution of the differential equation (B.5). 

Near the singular points any solution behaves as 53 

Ax{z - z x ) a + Bi(z - z x ) a ' as z -» z x 
/( : ) ~ { A 2 (z - z 2 f + B 2 (z- z 2 f asz^z 2 (B.8) 
A 3 (z - z 3 )^ + B 3 (z - z 3 y as z -> z 3 , 

so the monodromies are simply expressed in terms of a, a', /3, . . . 

B.4 Particular Solutions of Riemann's Equation 

We now construct explicit solutions of Riemann's equation that have simple mon- 
odromy at the three singular points in terms of the hypergeometric function (B.l). 
Given equation (B.7), the most obvious solution we can write down is 

f^\z)=( Z ^Y ( Z -^Y F(a + P + 1 ,a + P' + 1 ,l + a-a' 



z-z 2 ) \z-z 2 ) \ z l3 {z-z 2 )J 

We denote it f( a > because the holomorphy of the series (B.l) at ensures that any 
nontrivial monodromy near z\ comes only from the explicit factor (z — Z\) a . The 
differential equation is invariant under interchanging a f> a', so we can easily write 



53 As before, there is a caveat that when the two exponents are equal at one or more of the singular 
points, there are additional solutions involving logs. 
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down another solution that is linearly independent with the first (assuming that a 7^ 
a'): 

{z) = ( Z -^) ^ ( Z -^) 7 F (a + /3 + 7, a' + + 7, 1 + a' — a, Z J^A) . 
\z-z 2 J \z-z 2 J \ z 13 (z-z 2 )J 

This solution has the alternate monodromy around z — 0. 

The differential equation is also invariant under /?•<-)■ /3' and 7 H 7': the former 
leaves the solutions f^ a \ f^ a '^ invariant and can be ignored, but the latter apparently 
generates two additional solutions. We can find even more solutions by simultaneously 
permuting {z%,a, a'} -H- {z 2 ,f3,f3'} {^3,7,7'}, so combining these permutations we 
find a total of 4x6=24 solutions, known as "Rummer's Solutions" . 54 Since these are all 
solutions of the same 2nd-order linear differential equation, any three of them must be 
linearly dependent. 

To pin down this redundancy, it is convenient to define a particular set of six 
solutions, each of which has simple monodromy about one of the singular points. The 
definition is somewhat arbitrary as one can change the normalization at will as well as 
move around the various branch cuts. We will choose expressions that are simple when 
all ^-dependence is folded into the harmonic ratio 

g s ~ *0 (B 9) 

Zl3{Z ~ Z 2 ) 

54 This derivation of Kummer's Solutions using the symmetric equation (B.5) is quite straightforward, 
but if we had used the less symmetric equation (B.3) then they would seem quite mysterious. 
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Our explicit definitions are the following: 

p a ( x ) = z «(i _ x ) 7 F(a + /3 + 7, a + (3' + 7, 1 + a - a', x) 

= x a (l - x)- a - p F {^a + p + 7, a + p + 7', 1 + a - a', ^37) 

p a \ x ) = x a '(i - x y'F(a' + f3 + 7', a' + P' + 7', 1 + ol - a, x) 

= x a \l - x)- a '-^F (a' + (3 + 7, a' + /3 + 7', 1 + a' - a, 

py( x ) = x <*(\ - x ) 7 F(a + p + 7, a + p' + 7, 1 + 7 - 7', 1 - z) 

= - x) 7 F(a' + /5 + 7, a' + /?' + 7, 1 + 7 - 7', 1 - z) 
pV( x ) = x «(i _ x y'F(a + p + 7', a + p' + 7', 1 + 7' - 7, 1 - x) 
= x a \l - xy'F(a' + p + 7', a' + p' + V, 1 + 7' - 7, 1 - a) 

P^O) = - x)"""^ fa + /3 + 7, a + /3 + y, 1 + /3 - p', ^—^J 
= x a '(l - x)- a '-?F (a' + P + j,a' + P + i,l + P- P', — ?-j-J 

PP'(x) = x a (l - x)- a ~P'F (a + /?' + 7, a + /?' + 7', 1 + /?' - P, 

= x a '(l - x)- a '-?'F fa' + P' + 1) a' + P' + 7', l + P'-P, ^—^J . (B.IO) 

These formulas are somewhat intimidating, but they follow from the simple permuta- 
tions just described. For convenience in the following derivation we give two equivalent 
forms of each. More symmetric integral expressions for them will be described in section 
B.5. 

Since only two of these can be linearly independent, there must exist coefficients 
dij such that 

p<*( x ) = a Q7 P 7 (x) + a a yP 7 '(x) 

P a '(x) = a a ^P^(x) + a aV P 7 '(x). (B.ll) 
These coefficients are called connection coefficients. To determine them we can evaluate 
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these two equations at x = and x = 1, which gives 

r(l + a-a') r (y-7) 



T(a + + y)r(a + P' + f) 

r(i + a - a')r(7 - i) 
' r(« + /3 + 7 )r(a + /3' + 7 ) 

r(i + a' - a)r(V - 7) 
r(a' + /3 + y)r(a' + /3' + y) 
r(i + a' -a)r(7-y) 



(B.12) 



r(a' + /3 + 7 )r(a' + /3' + 7)' 
In solving these equations one uses (B.2). 55 Similarly one can find: 

r(l + a - a')Y{p' - p) 



r(a + /3' + 7)r(a + /3' + 7') 

r(l + a - a')T{P - p') 
" T(a + P + 7 )r(a + p + 7') 



r(i + a' - a)r(/3' - 
^ " r(a' + /3' + 7 )r(a' + /?' + 7') 1 J 

r(i + a' - a)r(/3 - /?') 



r(a' + p + 7)r(a / + /? + 7') 



Finally we note that our expressions for the connection coefficients allow us to 
derive some beautiful facts about the original hypergeometric function F(a,b,c, z). 
First making the replacements mentioned below (B.5), we see that (B.ll) gives: 

, \ r(c)r(c — a — b) „ , , , . 

Fa, 6, c, z = -(F a, 6, 1 + a + 6 - c, 1 - z 

r(c-o)r(c-6) 

r( c )r(a + &-<:),, ^ c _ a _ 6 ; ' 



r(a)r(&) 



(1 - z) c - a - b F(c - a, c - 6, 1 + c - a - 6, 1 - z). 



This gives explicit expressions for Ai(l — z) and B\{1 — z) for F(a,b,c, z), as 
promised above. We can also set a = 0, a' = 1 — c, /3 = 0, P' = c — a — 6, 7 = a, 7' = 6, 
Zl = 0, Z2 = 1, and z 3 = 00, in which cases B.ll gives: 

F(a, 6, c, z) = ^ C)r(6 ~^ (-z)- a F(a, 1 - c + a, 1 - 6 + a, z~ l ) 

T(c-a)T(b) (Rl5) 

+ r! C)r t"fV z )^ 1 - c + 6, 1 + 6 - a, z" 1 ). 
1 (c — o)l (a) 



55 Two identities which are useful are sin(x) sin(y) = sm(x + y — z) sin(z) + sin(z — x) sin(z — y) and 
r(x)T(l-x) = — — 7- — \ . 
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This expression gives A^l/z) and B^iA/z) for F(a, b, c, z), and in fact it gives the 
full analytic continuation of the series (B.l) in the region \z\ > 1, since the hypergeo- 
metric series on the right hand side converge in this region. We can thus observe that 
indeed the only singular behaviour of the function F(a, b, c, z) is a branch cut running 
from one to infinity. 

B.5 Integral Representations of Hypergeometric Functions 

We now consider the integral representations of the hypergeometric and P functions. 
We begin by defining 



I c (a,b,c,z)= I dss a - c (s-l) c - b -\s- z)~ a , (B.16) 
c 



where C is some contour to be specified in the s-plane. If we insert this integral into 
the hypergeometric differential equation (B.3) we get 



/ ds i~ [s a ~ c+1 {s - \y-\s - zy a - 1 } = 0, 

Jc ds 



(B.17) 



so this function will be a solution of the equation as long as C has the same and initial 
and final values for the quantity in square brackets. As an example we can choose C to 
run from from one to infinity, which is allowed when Re c > Re b > 0. The monodromy 
of the integrand as z circles zero is trivial everywhere on the contour, so we expect this 
solution to be proportional to the original hypergeometric function (B.l). Indeed we 
have 

T(r) f°° 

F(a,b,c,z) = ^y-^ ^ I dss a - c (s-l) c - b -\s- z)- a Rec>Re6>0. 

(B.18) 



r(6)r( c -6) h 

To establish this we can use the binomial expansion 

T(a + n) /z\ n 



1 " = Y tv t, } S - (B.19) 

^ T(a)T(n + 1) Vs' 



n=0 



to expand the integrand. We then change variables to t — 1/s and use Euler's integral 
for the Beta function 

p(x,y) = T ^ T ^ = [ dtt x ~\l-t) y - 1 Rex > 0, Rey > 0. (B.20) 
i[x + y) J 

This representation allows an easy determination of the value of the hypergeometric 
function at z — 1. Changing variables s — 1/t and using the /3 function integral we 
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Figure 4. The Pochhammer contour. 



find: 56 

F(a, b, c, 1) = ^ C)r( ° ~ ° ~ *| Re(c-a-b)>0. (B.21) 
1 (c — a)L (c — o) 

The main point however is that by integrating on other contours it is possible to 
get other solutions of the hypergeometric differential equation in a straightforward way. 
For example if we integrate from zero to z, which requires 2 > 1 + Re 5 > Rec> 0, 
then by changing variables to w = z/s it is easy to see that we get 

~1— CT71/1 i L. „ 1 I „ „ O „ ~\ T(2 C) /" t_ cf \ C — a— 1 



F(l+6-c, l+a-c2-c, s) = - fe v _ ' _ jf da^l-a) — ^z-a)" 6 , 

(B.22) 

which is the other linearly independent solution of the hypergeometric differential equa- 
tion with simple monodromy at z = 0. More generally there are four singular points 
of the integrand, and placing the contour between any two gives six different solutions 
which correspond to the six solutions that have simple monodromy at 0, 1, oo. 

Unfortunately these simple contours require strange inequalities on a, b, c to be 
satisfied which we certainly do not expect to hold for the general solutions we are 
considering in Liouville theory. To find contours that produce solutions for arbitrary 
a, b, c is more subtle. The trick is to use a closed contour that winds around both 
points of interest twice, but in such a way that all branch cuts are crossed a net zero 
number of times. This is called a Pochhammer contour, it is illustrated in Figure 4. If 
the inequalities we've been assuming are satisfied then we can neglect the parts of the 
contour that circle the endpoints and it collapses to the one that runs between the two 
points times a simple factor that depends on the choice of branch of the integrand. In 
general, if the inequalities are not satisfied, we can just use the Pochhammer contour. 
For example we can write 



56 0nc might guess that this formula should also require Rec > Re& > 0, but these conditions are 
a relic of our simple choice of contour. The Pochhammer contour we discuss below will remove these 
extra conditions, but it cannot remove the condition Re (c — a — b) > since from (B.4) we see that 
the function actually diverges at z = 1 if this is violated. 
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(B.23) 

where C is a Pochhammer contour involving 1 and oo and the extra factor j^^awz^ii 
cancels the sum over four traversals out to infinity and back. 57 This expression gives 
the full analytic continuation of the hypergeometric function in all of its parameters. 

We can also construct an integral representation of an arbitrary solution of Rie- 
mann's more general differential equation. We begin by defining 

p c =( z - Zl Y( z - Z2 f( z -z,y 

x f ds(s - Zl ) a,+ ^- x {s - z 2 ) a+ P'+^- l (s - z 3 ) a+ ^'- 1 (s - z)- a -^. (B.24) 
Jc 

By explicit substition into Riemann's equation and rewriting as a total derivative as in 
(B.17), we see that see that this integral will be a solution as long as the quantity 

V = (s - Zl ) a ' + ^(s - z 2 ) a+ P' + ^(s - z 3 ) a+ ^'(s - (B.25) 

has the same value at both ends of the contour. For a closed contour this means 
the contour must cross each branch cut of the integrand a net zero number of times. 
We may thus choose Pochhammer type contours involving any two of the four branch 
points Zi, Z2-, Z3, z. There are six such choices, and the integral evaluated on these six 
different choices is proportional to our P a ,P a defined by (B.10). The detailed 
proportionality depends on the branch choices and we will not work it out here. 

So far what we have gained in elegance over our previous formulation in terms of 
the series (B.l) we have arguably lost in the sophistication of the contours and the 
branch choices. But for 77's in the physical region this presentation allow a very nice 
repackaging of the Liouville solution (3.25). The quantity f(z,z) = uu — vv defined 
by (3.2) has an enlightening expression in terms of the types of integrals we have been 
considering so far. Our claim is that 

/ =fo\z - z 2 \ 2 \z - Zl \ 2a \z - z 2 \ 2 ?\z - z 3 | 27 

x J d 2 s\s - Zl \ 2 ^' + ^-% - z 2 \ 2 ^' + ^\s - z 3 \ 2( - a+ ^'-V\s - z|~ 2 (<*+^), 

(B.26) 

with the parameters a, a', . . . given by equation (3.19) and the integral being taken over 
the full s-plane. As long as this integral converges it is clearly monodromy invariant, 

57 The precise form of this factor depends on a choice of branch for the integrand. A different choice 
would multiply the integral by a z-indepcndent constant. 
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and it solves the holomorphic and antiholomorphic differential equations (3.4,3.5) by 
the same argument as just given for the integral (B.24). We saw in section 3.3 that these 
two properties were sufficient to uniquely determine / up to an overall normalization, so 
this establishes our claim. The conditions for the convergence of this integral, expressed 
in terms of the r/j's, are 

Re{ Vl + V2 - V3 ) <1 (B.27) 
Re(r/i + 773 - r/ 2 ) <1 
Re(r) 2 + V3- Vi) <! 

Re (r/i +T)2 + T) 3 ) >1. 

These are certainly obeyed in the "physical" region in Liouville. Are they equivalent to 
it, or more precisely to Region I from section 4.2? Actually, they imply < Rer/j < 1 
for all i, while in Region I we would have had < Rer/j < 1/2. But for operators 
obeying the Seiberg bound, the integral converges only in Region I. For more general 
rji, such an expression would require a more sophisticated type of integral. For real 77's 
in Region I, this expression is manifestly positive and it shows that there cannot be any 
zeros of /, something that was not clear from our old expression (3.25). We wonder if 
Liouville solutions for correlators with more than three heavy operators can be written 
in terms of generalizations of this integral. 



C Gamma Functions and Stokes Phenomena 
C.l Generalities 

In this appendix, we review the Stokes phenomena that occur for the Gamma 
function and its reciprocal, as they are closely related to the zero mode integrals of 
spacelike and timelike Liouville theory. 58 T(z) has the following integral representation 
[56] for Re z > 0: 

poo POO 

T{z) = / f-^dt = z z e-'^-^dcj), (C.l) 

(With a slightly different change of coordinates by t = rather then t = ze^, we could 
have put this in the form of the Liouville zero mode integral, as in eqn. (1.11).) 
The exponent X = — z[e^ — <fi) in (C.l) and (C.14) has critical points at 

<f) n = 2mn, neZ. (C.2) 

58 For a much more detailed introduction to Stokes phenomena, see section 2 of [55]. For a treat- 
ment of the Gamma function along lines similar to what follows, see [12]; see [13, 14] for previous 
mathematical work. 
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To each such critical point n , one attaches an integration contour C n . This is a contour 
that passes through the critical point <p n and along which the exponent X has stationary 
phase while ReX has a local maximum. More briefly, we call this a stationary phase 
contour. Alternatively, 59 the contour C n can be defined as a contour of steepest descent 
for h = Re (— z(e^ — <p)). For a steepest descent contour C n , it is straightforward to 
determine the large z behavior of the integral 



The maximum of ReX along the cycle C n is, by the steepest descent condition, at 4> n . 
For large z, the integral can be approximated by the contribution of a neighborhood 
of the critical point. In our case, the value of X at a critical point is — z(l — 2nin), so 
asymptotically 



(times a subleading factor that comes from approximating the integral near the critical 
point). 

Now let us consider the integral (C.l, initially assuming that z is real and positive. 
The Gamma function is then defined by the integral (C.l), with the integration cycle 
C being the real <fi axis. On the real axis, there is a unique critical point at = 0. 
Moreover, for real z, X is real on the real axis, and the contour of steepest descent from 
= is simply the real axis. Thus, if z is real and positive, the integration cycle C 
in the definition of the Gamma function is the same as steepest descent cycle Co, on 
which the asymptotics are given by (C.4). So we get the asymptotic behavior of the 
Gamma function on the real axis: 



This is essentially Stirling's formula (the factor l/\ / 2nz in Stirling's formula comes 
from a Gaussian approximation to the integral (C.4) near its critical point). 

Now let us vary z away from the positive z axis. The Gamma function is still defined 
by the integral (C.l), taken along the real (f) axis, as long as Rez > 0. As soon as z is 
not real, it is no longer true that the steepest descent contour Co coincides with the real 
axis. However, as long as Re z > (we explain this condition momentarily), the steepest 

59 In complex dimension 1, the stationary phase contour through a critical point coincides with the 
steepest descent contour, but in higher dimension, the stationary phase condition is not enough to 
determine C n and one must use the steepest descent condition. For much more on such matters, see 
section 2 of [55]. Notice that, in our case, because the function X has opposite sign in (C.14) relative 
to (C.l), the cycle C n is different in the two cases, though we do not indicate this in the notation. 




(C.3) 




(C4) 




(C.5) 
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descent contour C is equivalent to the real axis, modulo a contour deformation that is 
allowed by Cauchy's theorem. Hence, Stirling's formula remains valid throughout the 
half-plane Rez > 0. 

If we want to analytically continue the Gamma function as a function of z, in general 
we will have to vary the integration contour C away from the real axis. To analytically 
continue beyond the region Rez > 0, we can let the integration contour C move away 
from the real <fi axis, so that the integral still converges and varies analytically with z. 
In the case of the Gamma function, there is some restriction on the ability to do this, 
since the Gamma function actually has poles at z = and along the negative z axis. 

Now we come to the essential subtlety that leads to Stokes phenomena. As one 
varies the parameters in an integral such as (C.l), the steepest descent contours C n 
generically vary smoothly, but along certain "Stokes lines" (or Stokes walls in a problem 
with more variables), they jump. In our case, the only relevant parameter is z, so the 
Stokes lines will be defined in the z-plane. For generic values of z, the C n are copies of 
K (topologically) with both ends at infinity in the complex <p plane. For example, for z 
real and positive, C n is defined by Im z = 2rm and actually is an ordinary straight line 
in the <fi plane. However, for special values of z, steepest descent from one critical point 
p leads (in one direction) to another critical point p'. Whether this occurs depends 
only on the argument of z, so it occurs on rays through the origin in the z plane; these 
rays are the Stokes lines. As one varies z across a Stokes line £, the steepest descent 
contour from p will jump (on one side of i, it passes by p' on one side; on the other 
side of £, it passes by p' on the other side and then heads off in a different direction). 

For the Gamma function, we can easily find the Stokes lines. Since the steepest 
descent cycles have stationary phase, they can connect one critical point p to another 
critical point p' only if the phase of I is the same at p and at p' . For the critical point 
at = 2'Kin, the value of ImX is c n = Im (—2(1 — 2irin)). So c n = c n > for n ^ n' if and 
only if Re z = 0. We really should remove from this discussion the point z = where 
our integral is ill-defined for any noncompact contour (and the Gamma function has a 
pole), so there are two Stokes lines in this problem, namely the positive and negative 
imaginary z axis. 

There is one more basic fact about this subject. Away from Stokes lines, the 
steepest descent contour C n are a basis for the possible integration cycles (on which 
the integral of interest converges) modulo the sort of contour deformations that are 
permitted by Cauchy's theorem. So any integration contour C - such as the one for the 
analytically continued Gamma function - always has an expansion 




(C.6) 



n 
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where the a n are integers, and the relation holds modulo contour deformations that 
are allowed by Cauchy's theorem. Since the integral over any of the C n always has the 
simple asymptotics (C.4), the asymptotics of the integral over C are known if one knows 
the coefficients a n . As one varies z in the complex plane, C will vary continuously, but 
the C n jump upon crossings Stokes lines. So the asymptotic behavior of the integral 
for large z will jump in crossing a Stokes line. The well-behaved problem of large z 
asymptotics is therefore to fix an angular sector in the complex z-plane between two 
Stokes lines and consider the behavior as z — > oo in the given angular sector. Actually, 
this would be a full picture if there were only finitely many critical points. In the case 
of the Gamma function, it will turn out that the sum (C.6) is an infinite sum if Re z < 
and moreover this infinite sum can diverge if z is real and negative. To get a simple 
problem of large z asymptotics, one must keep away from the negative z axis, where 
the Gamma function has its poles, as well as from the Stokes lines. 

C.2 Analysis Of The Gamma Function 

Now let us make all this concrete. The integral (C.l) converges along contours 
that begin and end in regions where Re(— z(e^ — 0)) — > — oo. These regions have been 
shaded in Fig, 5 for the case that z is real and positive. The C n are the horizontal 
lines Im0 = 2im. One can see by hand in this example that any integration cycle 
that begins and ends in the shaded regions is a linear combination of the C n , as in eqn. 
(C.6). The real </> axis - which is the integration contour Co in (C.l) - coincides with 
Cq, and is indicated in the figure as the Relevant Contour. 

In Fig, 6, we have sketched how the steepest descent contours C n are deformed 
when z is no longer real but still has positive real part. In passing from Fig, 5 to Fig, 6, 
the C n evolve continuously and it remains true that the integration contour C defining 
the Gamma function is just Co, modulo a deformation allowed by Cauchy's theorem. 

However, the C n jump upon crossing the Stokes lines at Rez = 0. This is shown 
in Fig. C.2 for the case Imz > 0. While in Fig, 6, the steepest descent contours for 
the Gamma function have one end in the upper left and one end to the right, in Fig. 
C.2, they end to the right in both directions. As a result, although nothing happens 
to the contour C that defines the Gamma function in going from Fig, 6 to Fig. C.2, to 
express C as a linear combination of the C n 's, we must in Fig, C.2 take an infinite sum 




(C.7) 



n>0 
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Figure 5. For the Gamma function integral (C.l) to converge, the integration contour must 
begin and end in regions of the eft plane with Re(— z(e^ — <fi)) — > — oo. These regions are 
shaded here for the case that z is real and positive. In addition, we show the critical points 
at (j> n = 2-irin (represented in the figure by dots) and the steepest descent cycles C n , which 
are the horizontal lines Imcf) = 2im. 



The dependence on n of the Gamma function integral over C n is very simple: 

/ d0exp(-z(e* - <p)) = exp{2mnz) / d<p exp{-z(e* - <p)). (C.8) 

This is because a shift — > <f> + 2nin maps Co to C n and shifts X by X — > X + 2mnz. 
This formula is an analog for the Gamma function of eqn. (6.30) for Liouville theory. 
Using (C.7) and (C.8), we find that in the quadrant Rez < 0, Imz > 0, the Gamma 
function is 

F(z) =z z ! d0exp(-z(e* - <p)) = z z J2 \ d(f)exp(-z(e^ - 0)) (C.9) 

JC n >Q JCn 

OO „ 

^exp(27rmz) / rf0exp(-2(e^ - <p)). (CIO) 

n J Cn 



-Z 

n=0 



From (C.9), we can read off the asymptotic behavior of the Gamma function in the 
quadrant in question. Approximating the integral over Cq by the value at the maximum, 
and performing the sum over n, we get 

T{z)~z z e- Z - l — — . (C.ll) 

1 — exp{2mz) 
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Figure 6. This figure is equivalent to Fig. 5, except that now z is complex but still with 
Rez > 0. (In drawing the figure, we have taken the case Imz > 0.) The critical points are 
unchanged from the case that z is real, but the steepest descent contours are changed. The 
phrase Relevant Contour labels the contour Cq that controls the asymptotics of the Gamma 
function in this region. 

(Again, a prefactor analogous to the \j\/2ixz in Stirling's formula can be found by 
evaluating the integral over Cq more accurately.) From this point of view, the poles 
of the Gamma function at negative integers arise not because of a problem with the 
integral over C but because of a divergence of the geometric series. The factor 1/(1 — 
exp(2niz)) in this formula is important only near the negative real axis. 

For use in the main text of the paper we note that a similar analysis of the case 
where Kez < 0, Imz > gives 

T(z) ~ z z e~ z - -. — - — — , (C.12) 

and that these formula can all be combined to give 

_ je zlogz - z+0< - losz ^ Re(z) > 

r (*) = i i re *iog(-*)-*+o(io g (-*)) Re(z)<0, ,(JO; 



alTTZ o — ITVZ ' 



where the logarithms are always evaluated on the principal branch. 
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Figure 7. This is the analog of Fig, 6, but now for the case that Hez < 0, Imz > 0. (For 
Hez < 0, Imz < 0, just turn the figure upside down.) A Stokes phenomenon has occurred, 
relative to Fig, 6. In Fig, 6, each steepest descent curve for the Gamma function connects 
the shaded region in the upper left to one of the shaded regions on the right. This is also the 
behavior of the contour C that defines the Gamma function. However, for Re z. Imz < 0, the 
steepest descent contours connect two adjacent shaded regions on the right. To construct the 
contour C that controls the Gamma function, one must take an infinite sum ^2 n>0 C n . 

C.3 The Inverse Gamma Function 

We can play the same game for the inverse of the Gamma function, starting with 
the integral representation 

J- = -L / t- z e*dt = —z- z+1 { e*( e ~* + *>e-+d<l>. (C.14) 
T(z) 2mJ Ct 2m J c V K > 

In (C.14), Ct starts at real — oo — ie, encircles the branch cut along the negative real t 
axis, and ends up at — oo + ie. 

To arrive at the right hand side of (C.14), we have made the coordinate change 
t = ze~^, which differs slightly from the transformation t = ze^ used in deriving 
(C.l). The critical points are still at <ft n = 2nin. Once again the shaded regions in 
Fig, C.3 are the ones in which the integral is convergent, as ReX — > — oo, where now 
X = z(e~^ + 0). The steepest descent contours are shown in Fig, C.3 and connect 
adjacent shaded regions on the left of the figure. For z real and positive, the image 
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Figure 8. For the inverse Gamma function integral, with z real and positive, the critical 
points and steepest descent contours are as shown here. Regions in which the integral is 
convergent are again shaded. The function 1/T(z) is defined by an integral over a contour C 
that coincides with the steepest descent contour Co associated to the critical point at z = 0. 
This contour connects two adjacent shaded regions on the left of the figure. 



in the plane of the contour C t in the t plane is the steepest descent contour Co that 
passes through the critical point at = 0. Since C is a steepest descent cycle, the 
asymptotic behavior of the integral in this region is just e x , with X evaluated at the 
critical point. So in this region, 

1 



T(z) 



~ z- z+l e z . 



(C.15) 



(As always, subleading factors, including powers of z, can be determined by evaluating 
the integral more accurately near the critical point.) 



When we vary z, the regions in the <p plane in which the integral converges move 
up and down slightly. The integration contour C defining the inverse Gamma function 
varies smoothly and continues to connect two adjacent shaded regions. However, when 
Re z < 0, the qualitative behavior of the steepest descent contours C n changes. As 
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Figure 9. For Rez < and Imz > 0, the steepest descent contours C n connect a region on 
the left to the upper right, as shown here. On the other hand, after varying z, the integration 
contour C that defines the inverse Gamma function still connects two adjacent regions on the 
left. (It is not drawn here.) It is therefore no longer true that C equals one of the steepest 
descent contours; rather, C is a difference Co — C\. 



sketched in Fig, 9, for such values of z, each C n connects a shaded region in the left of 
the (f) plane to the upper right. The formula expressing C in terms of steepest descent 
contours is now C = C — C\. (Here C — C\ starts in a shaded region on the left, heads 
to the upper right, and then doubles back to an adjacent shaded region on the left - 
thus reproducing C.) Accordingly, (C.15) is modified to 

' (1 - e 27Tiz )z~ z+1 e z . (C.16) 



T(z) 



Of course, this result for the asymptotic behavior of 1/T(z) is equivalent to the 
result that we had earlier for the asymptotic behavior of T(z). However, seeing this 
behavior directly from the Stokes phenomena that affect an integral (rather than the 
inverse of an integral) is useful background for the body of this paper. 

Just as in the relation between spacelike and timelike Liouville theory as studied 
in this paper, the Gamma function and the inverse Gamma function are essentially 
given by the same integral, evaluated on different integration contours. This fact is 



- 107- 



also related to the functional equation obeyed by the Gamma function. We have 



2tu 



r(i-z) 



**-Vtft. (C.17) 



The integrand on the right hand side can be converted to the integrand of the Gamma 
function integral (C.l) if we substitute t — > —t. This maps Ct to another integration 
contour C' t and gives 



2iri 



r(i-z) 



exp(m(z - 1)) / t'^e^dt. (C.18) 



So the inverse Gamma function, apart from some elementary factors and a substitution 
z — > 1 — z, is given by the same integral as the Gamma function, but with a different 
integration contour. 



D Semiclassical Conformal Blocks 

In this appendix we give a derivation, based on [48] , of an asymptotic formula for the 
Virasoro conformal block at large intermediate operator weight. The original argument 
was somewhat terse and implicitly involved certain assumptions about the semiclassical 
limit of correlators, so in our view its validity has not been established completely 
rigorously from the definition (2.40). It has however survived stringent numerical tests 
[7], and in our discussion in section 5.5 we were comfortable assuming it to be true. 
We will expand out the argument here and try to be clear about what the assumptions 
are. The asymptotic formula is 

^"1234(^4, A, x) ~ (16g) A , (D.l) 

with 

q = exp[-nK(l-x)/K(x)] (D.2) 

and 

K( - 1 f 1 dt 
{X) ~2J ~Jt(l-t)(l-xt)' ' ) " >1 

The idea is to study a five-point function of a light degenerate primary operator 
with four primary scalar fields of generic operator weight in the semiclassical limit 
c 3> 1. The external conformal weights Aj are taken to be of order c, and the internal 
weight A p is initially also taken to be of this order; it will eventually be taken much 
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larger than c. 60 If we parameterize the central charge as c = 1 + 6(6 + 1/b) 2 , then 
the conforma 
correlator as 



the conformal weight of the light degenerate operator is — \ — We will write the 



(0 4 (Z4, z 4 )0 3 (z 3 ,z 3 )*(z, z)0 2 (z 2 , z 2 )O x { Zl , Zl )), (D.4) 

where ^ is the degenerate operator and we choose |z 4 | > \z 3 \ > \z\ > \z 2 \ > \zi\- This 
correlator obeys equation (2.43), which here gives 

4 

6 2 



{0 A 0^l0 2 O x ) = 0. (D.5) 



61 



We can also expand the correlator using the operator product expansion [17] 

2 (z 2l z 2 )0 1 (z 1 ,z 1 ) = J2c P 21 \z2i\ 2{Ap - A '- A2) Y,^) k ^f^^°l k 

p k,k 

(D.6) 

which allows us to extract all dependence on z 2 . 

(0 4 3 *0 2 0i) =^c 2 p 1 |^ 21 | 2 ( A ^ Ai - A2 ) 

V 

x ^(^i) fc (^ 2 i)^/3 2 f (0403*0? •*>>. (D.7) 

In these formulae, as discussed below equation (2.40) the sum over k is only heuristic 
and more precisely includes a sum over all descendants at a given level k. The operators 
are the Virasoro descendants of the primary O p . Now say that we define 

*p(*> *>«>**)= (o 4 o 3 o p ). (D - 8) 

In the semiclassical limit we can think of the function ty p as the classical expectation 
value of the degenerate operator in the presence of the other operators, and since the 
light operator has weight of order c° we expect it to have a finite limit at c — > oo. In 
Liouville theory, this is the statement that light operators just produce O(b ) factors 
in the correlation function as in equation (2.26). 



60 This argument could be applied to any CFT with a c that can be large and primary operators with 
the desired weights. "Light" in this context just means that the operator weight of the degenerate 
field scales like c°. Liouville is a theory that fits the bill, but we will use general CFT language to 
avoid the subtleties of Liouville factorization. 

61 We assume here that as in Liouville, the only primaries that appear in the operator product 
expansion are scalars. We could drop this assumption at the cost of slightly more complicated formulas, 
the result (D.l) would be the same. 
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So far we have written only exact formulas, but we now come to the first ap- 
proximation: in the semiclassical limit we claim that the same formula holds also for 
descendants of O p with the same function ty p : 

(0 4 3 *Oj, k &) « * P {z,z; zt,zi)(0 4 3 0^). (D.9) 

The justification for this is that the correlator (O^Os^Op*'^) can be written in terms 
of a series of differential operators acting on [O^O^O p ). The differential operators 
are of the form [17] 

(m - l)Aj 1 



C = V 



-d- 



(D.10) 



i=3,4,2 

Here A 2 = — | — ^p. Similarly the correlator (0^0 3 O p k,k ^) can be written in terms of a 
similar series of differential operators acting on {0±0 3 O p ), but with the sum in (D.10) 
being only over 3,4. The point however is that in the semiclassical limit we expect 

(0 A 3 O p ) ~ e~i^ ( D . n) 

for some S c i, and since we have taken A 3)4 ~ c while ty p and A 2 are both O(c ), we see 
that the i = z term in (D.10) becomes unimportant, and also that in the % = 3, 4 terms 
we can neglect derivatives acting on This establishes (D.9). In the semiclassical 
approximation we thus have 

(0 4 3 ^f0 2 1 ) « ^ Cf^xl 2 ^-^-^^, Z 4 , Z 4 , *3, *3, *1, *l) 

x ^(^1)^^21)^^(0403^^). (D.12) 

k,k 

We can make this formula more elegant by defining 

F 12 34(A,, A p , Zi ) = (z 21 )^^ D^ ^? ' (D.13) 
after which we get 

(0 4 3 ^0 2 0i) ~ ^VpCS^OtOaOJF^Ai, A p , Zi )F 1234 (A h A p ,^)- (D.14) 
p 

We note for future convenience that .F1234 becomes the Virasoro conformal block J-1234 
after sending z 4 — > 00, z 3 — > l,z 2 —tx,Zi —> 0. Based on its definition, we can guess 
that -F1234 has a semiclassical limit [48] of the form 

F 123 4 ~e-f /d , (D.15) 
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where f c \ is called the "semiclassical conformal block". 62 

We will now study the implications of the expressions (D.14, D.15) for the differen- 
tial equation (D.5). It will be convenient to view z and z as independent. We observe 
that for generically different A p 's, the various terms in the sum over p in (D.14) have 
different monodromy as z 2 circles z\. For each p the different terms in the differential 
equation (D.5) have the same monodromy, so in order for the equation to be solved by 
(D.14) it seems reasonable to expect that it must actually be solved separately for each 
p. In the semiclassical limit the action of the derivatives with respect to z\ on is 
suppressed as in our discussion below (D.10), so we find that the differential equation 
can be converted into an ordinary differential equation just involving ty p : 

* p = 0, (D.16) 
with 

C i = d i {S d + f dl ). (D.17) 

Here 5i = 6 2 Aj. In this type of differential equation the parameters Cj are referred to 
as "accessory parameters"; clearly if we can learn something about them then we are 
learning about the semiclassical conformal block. In [48] it was shown that for A>c, 
a combination of symmetry and the WKB approximation allows a determination of all 
Cj, and thus of the semiclassical conformal block in that limit. 

This argument begins with the observation that there are three linear relations on 
the accessory parameters which come from demanding that the term in round brackets 
in (D.16), which is related to the semiclassical limit of the stress tensor by the Ward 
identity (T(z)0 4 3 2 1 ) = £ti + ^ft) (0 4 3 2 1 ), vanishes like z" 4 at 

infinity: 

i 

{CiZi - 5,) = 

i 

(Cizl - 2S iZi ) = 0. (D.18) 



62 In fact this exponentiation has never actually been proven directly from the definition (2.40), 
although it has been checked to high order numerically. We thank A. B. Zamolodchikov for a discussion 
of this point, and for a summary of his unpublished numerical work. 




- Ill - 



There is thus only one independent accessory parameter, which we take to be C 2 . If 
we then take the limit z± — >• oo, Z3 — > 1, Z2 — > x, Z\ — > 0, equation (D.16) becomes 



o 5i S 2 5 3 61 + S 2 + 5 3 - 5 4 C 2 x(l - x) 



dt + ^ + — + -rr^-, + 



(D.19) 



z 2 (z — x) 2 (1 — z) 2 z(l — z) z(z — x)(l — z) 

In the same limit (D.17) simplifies to 

c 2 = d x f cl . (D.20) 

One way to parametrize the effect of C 2 is to study the monodromy of ty p as z 
circles both x and z\. We will work out this relationship below, but we first note 
that this will give us what we want because we can also determine this monodromy 
from the definition of The reason is that as discussed in section 2.4, the four- 
point function (O^O^Op) receives contributions only from two intermediate conformal 
weights. If we parametrize conformal weights as A = a(b + 1/b — a) then these are 
A± = a±(b + 1/b — a±), with a± = a p ± b/2. These contributions behave near z = 
like Z A ±- A p- A * ) which semiclassically becomes z 2 ( 1± "\/ 1 ~ 46 A p) ? so their monodromy 
matrix in this basis is 




^=1 q (D.21) 

So far everything we have said is valid for A p ~ c, but we now observe that if we take 
A p 3> c then we can solve (D.19) in the WKB approximation, where we include only 
the first and last terms. This gives approximate solution: 



\I/p ~ exp 



±y/x(l-x)C 2 



dz' 



X) 



y/z'{l- Z '){z! 

Comparison with the mondromy matrix (D.21) in the same limit gives 



(D.22) 



C2 ~~x(l-x)K(xy [DM} 

Finally we can integrate this by observing that K(x)d x K(l — x) — K{\ — x)d x K(x) = 
— ix ^_ x ^ , which follows from the Wronskian of the hypergeometric differential equation 
obeyed by K(x) = | F(l, |, 1, x), with the normalization determined by expanding near 
x = 0. This gives 

K(l — x) 

f d « 7r& 2 A„ v - - ; + constant. (D.24) 

AX 
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We can determine the constant to be — b 2 A log 16 by matching the series expansion 
near x = to the normalization of the leading term in the conformal block 



J- 1234 (A,, A p , x) = x A -- A ^ (1 + 0{x)) , (D.25) 

which at last gives 

^1234 ~ (16g) A . (D.26) 



E An Integral Expression for logT(z] 
In this appendix we derive the identity 



iogr(z) 



dt 
T 



(z - l)e" 



-zt 



Re\z] > 0. 



(E.l) 



We by differentiating the definition (C.l) of the Gamma function with respect to 



r'(z) 



dj™t z - l e- l dt 



dz 



dtt z ~ x e- 1 



t* -1 log |t|e -t tft 



(E.2) 



ds , 



lim 

e-S>0 



dtt z - x e' 1 



e - e 



(is 



-stl 



lim 

e->0 



lim 

e-5-0 



dtt z - l e- 1 



ds , 



e - e 



-stl 



dtt 



2-1 



ds 
s 



(E.3) 
(E.4) 



Here we have used integral identity log b = J °° ^f(e~ x — e~ bx ). While the integral 
currently is convergent, we have to reexpress the lower limit in terms of e in order to 
break the integral up into divergent pieces. We do this to allow us to enact coordinate 
changes on the second term. The divergences will cancel each other out. We make the 
coordinate change p = t(l + s) on the second term. The change is linear in t. We 
exchange the p and s integrals and evaluate the t and p integrals giving 



lim 

<E-S>0 



dtt^e- 1 



T(z) lim 

£->0 



ds 



lim 

e-S>0 



ds 
s 



dtt z - x e- l e' st 



o 



r(z) 



lim 

e->0 



3 ds 
— e 
s 

00 ds 



lim 

e-S-0 



00 ds 
s 

oo 



dp 



P 



z-l 



— e " — lim 

s <^o 



1 + s \s + 1 
ds ( 1 



8(1 + 8) \S + 1 



(E.5) 
(E.6) 
(E.7) 



We see that the Gamma function factors out of the integral representation of V(z). 
To continue, we need to make one more coordinate change before putting the integrals 
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back together 1 + s = e s which yields ds = e s ds. Since e — > 0, when s ~ e then 
l + s = e s = l + s' + C^s 2 ) ~ 1 + e and s ~ e. The lower limit remains unchanged in 
the second integral. (E.7) then becomes 



= r 


» 


lim 


/. " e 


— lim 










= r 




lim 


f> 


— lim 










e->0 


= r 




lim 




— lim 













s(l + s) Vs + 1 



z-l 



(e s - l)e s 



e 



S\Z— 1 



(E.8) 
(E.9) 
(E.10) 



(1-e-J 

Now we can put the integrals back together and take the limit e — > 0, resulting in 



T(z) 



ds 



1-e- 



Integrating (E.ll) with respect to z yields (E.l 

S'(z) r ds 



\ogT(z) 



z-le" 



(E.ll) 



(E.12) 



F Integral over the SL(2,C) Moduli of the Saddle Point for 
Three Light Operators 



In section 4.3, we encountered the following integral: 

dfi(a, (3, 7, 5) 



I(ai, 02,03) = J 



2(72 / \ 2(73 ' 

\a\ 2 + I7I 2 



(F.l) 



'W + W) ai (|« + /3| 2 + |7 + ^| 2 ) 
In the text we claimed this integral is given by 

3 r(cri + cr 2 - cr 3 )r(cr 1 + 0-3 - o-2)r(cr 2 + 03 - 0"i)r(oi + 2 + 3 ~ 1) 



/(01, 02,03) 



r(2o!)r(2o 2 )r(2o 3 ) 



(F.2) 

and in this appendix we will show it. We will see along the way that the integral is 
divergent unless certain inequalities involving the o^'s are satisfied, so we will assume 
them as we come to them and then in the end define the integral away from those 
regions by analytic continuation. 

5> ^ 2 ~~ 7+5' 



Following [7], we begin by performing the coordinate change £1 = 4 , £2 = ttt, and 



£3 = -. The measure becomes dfi(a, p, 7, 8) 



l(«i-6)(6-6)(?3-Ci)l 2 ' 



and the integral 
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becomes 

x (1 + Ki| 2 )- 2CT1 (1 + !6lT 2ff2 (l + le 3 | 2 )- 2<T3 (F.3) 

Here v\ = <j\ — o~2 — 03, v<i = 02 — <J\ — 03, z/3 = 0-3 — o~i — 02 ■ This expression is invariant 
under the SU(2) subgroup of SL(2,C) given by & ->■ -^±^, with |/| 2 + |#| 2 = 1. We 
can use this to send £ 3 — > oc: 



IK<7 2 ,<7 3 )=7T | rf 2 ^ 2 6Kl-6r 2 ^ 3 (l + Kl| 2 )- 2,71 (l + l6| 2 )- 2ff2 . (F.4) 

From here the result is quoted in [7] without further explanation, we will fill in the 
steps. The evaluation will involve repeated use of the defining integrals of the T and (3 
functions, which we reproduce for convenience: 



POO 

T(x) = / t x - l e- l dt Re [a] > 
Jo 



(F.5) 



P(x,y) = ^ g)r( ^ = C \ x -\l-\) y - 1 d\= dt , ^ ' (F.6) 
^ K ,y) T(x + y) 7o Jo (l + ty+y 1 J 

Re [a] > 0, Re[y] > 0. 

We will also need a lesser-known version of the Feynman parameters which is used in 
closed string theory. 63 

1 1 f°° 

— = -—— / dttAP^e-W with Re A >0. (F.7) 
\z\ A T(A/2)J 1 ; 

Now to business. Starting with integral (F.4), we convert it into four real integrals 
defined by the coordinate change 

£1 = x + iy £ 2 = u + iv. (F.8) 

Noting that d 2 ^i = ^{d^i A = ^(dx + idy) A (dx — idy) = dx Ady with a similar 
identity for d 2 ^2, we find equation (F.4) becomes 



I = 7i I dxdydudv((x-u) 2 + (y-v) 2 )- 1 - i ~ as - ai -^\l + x 2 + y 2 )- 2ai (l+u 2 + v 2 )- 2a2 . 

(F.9) 



This identity can be easily derived by changing variables to t = t\z\ 
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We now use the identity (F.7) three times with A = 4<Ti,4<T2,2 + 2z/ 3 , which requires 
Re <7i > 0, Re cx 2 > 0, Re (a% + a 2 — 0-3) < 1, to get: 

dxdydudv / difjdxdK—— (F.10) 

x exp [-{if)[(x - w) 2 + (y - v) 2 } + x (l + x 2 + y 2 ) + re(l + w 2 + t> 2 )}]. 
Collecting powers of x,y,u,v we have 

dxdydudv / dipdxdn—— x Wo v Wo ; (F.ll) 

Jo r(i + cr 3 - 0-1 - o- 2 )r(2cri)r(2cr 2 ) 

x exp [-{(x 2 + y 2 )(^ + x) + ( u2 + *> 2 ) W> + «) - 2*l>(ux + yv) + x + «}]■ 
Completing the square for x and y and factoring out (if) + x) gives 

/oo /-oo ^0-3-0-2-0-1^,20-1-1^.20-2-1 
dxdydudv / difjdxdn—— v Wo s (F.12) 
Jo r(i + 0-3 - o x - «7 2 )r(2«7i)r(2«7 2 ) 



x exp [-(V> + X){(z - 



2/ 



^ + v V' + X 



x exp -(if) + x)i (V + v )-y (m 2 + f 2 ) 



+ « / 2 , 2a/ ^ \ , X + « 



ip + X \(^P + x) 2J ^ + x- 

The x, y integral is now a straightforward Gaussian integral. Changing variables by the 



linear shift s = x — j^- , t = y — f^- and performing the integral over s and t we find 



^03 — 02 — 01^,201 — 1^.202 — 1 

I = 7T 2 / d-ucfc / dibdydn— — — — — — — (F.13) 



x exp [— < (m 2 + w 2 ) 



V> 2 

■0 + K - 



*P + X 

We can also do the u and v integrals, which give 64 



+ X + K 



7j-3 /-oo ^,03-02-01 y2<ri— 1^.2(72— lg—x—« 

dx dn 



r(i + a 3 -a 1 - C r 2 )r(2 ( T 1 )r(2 ( T 2 )7 r + + - V] 

(F.14) 

So far the required manipulations have been fairly obvious, but to proceed further 
we need to perform a rather nontrivial coordinate change on if), x, an d We will 



64 For the u, v integrals in (F.13) to converge the prefactor of the u and v exponetial terms in (F.13): 
[■0 + k— ^p^\ > 0. It can be shown that this is so, by putting everthing over a common denominator, 
and noting that the term becomes a sum of positive quantities. 
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motivate it by answer analysis of (F.2). Summing the exponents of the ip, x, an d k 
terms in the numerator of (F.14), we get one minus the argument of the first Gamma 
function in (F.2). This implies that in the new coordinates, ip, Xi an d « must have some 
common factor p in order to produce this first Gamma function. This also means p must 
equal x + K due to definition (F.5). Now \ and k are linearly independent, therefore we 
propose the coordinate change: x = pcos 2 9 and k = psin 2 9. To determine ip we see 
that if we take the arguments of last two Gamma functions in the numerator of (F.2), 
they add up to 2cr 3 . That means the factor r ( CT3+CT2 ~^|)j^ 3+(Jl ~' j2 ) — f3{a 3 + cr 2 — o\ , 03 + 
<j\ — o<i) must be a factor in the integral. This Beta function will involve factors the of 
sin 2 9 and cos 2 9 from x an d k. A quick glance at the factors of x an d k in (F.14) shows 
that they are not correct. However, if ip instead included a factor of cos 2 (6) sin 2 (8), 
we would get the proper factors for the Beta function in terms of the integral over 9. 
ip must then have one factor of p and one factor of cos 2 (6 l ) sin 2 (6 l ). Lastly ip must be 
linearly independent from x an d k since we want a one to one mapping, so ip must 
have a third factor which we will call (. The coordinate change is then 

ip = pCcos 2 (#)sin 2 (#) x = pcos 2 (9) n = psm 2 (9). (F.15) 

The Jacobian for this is 

dipdxdn = 2p 2 cos 3 (fl) sm 3 (9)dpd9d(. (F.16) 

After the performing the transformation and collecting common terms, (F.14) be- 
comes 

2tj-3 roo roc r-rr/2 CT1 -|_ CT2 4. CT3 _l_l£l_|_ CT 3_ CT2 _ (Tl _l 

d( dpd9 



T(l + a 3 -a 1 -a 2 )T(2a 1 )T(2a 2 ) J J J r (1 + C) 

x (cos0) 2CT3+2CT1 - 2CT2 - 1 (sin0) 2CT3+2<T2 - 2ffl - 1 exp[-p]. (F.17) 

We now change the 9 coordinate to A = sin 2 9 which gives d9 = — , dX and thus 

^3 POO POO pi ^cri+0- 2 +CT3-l-l ( Al+CT3-(T2-0-l-l 

I= T{l + a,-a 1 -a 2 )T{2a l )T{2a 2 )j J J o ^ d P dX 

x (1 - A) ,T3 + ,71 - ,J2 - 1 A <T3+,J2 - <T1 - 1 exp [-p]. (F.18) 

These transformations have completely factorized the integral, which we can at last 
evaluate using (F.5) and (F.6). The result is: 

j _ ^ 3 r(ori + a 2 + q- 3 - i)r((xi + cr 2 - q~3)r(a 2 + - q~i)r(o- 3 + - q 2 ) 

r(2a 1 )r(2a 2 )r(2a 3 ) ' 



- 117- 



In addition to the inequalities mentioned below (F.9), in doing these final integrals 
we had to assume that 

Re(a 1 + a 2 - <7 3 ) >0 (F.20) 
Re(<Ji + 03 — a 2 ) >0 
Re((T 2 + cr 3 - 0-j ) >0 
Re (01 + 2 + 03) >1- 

These inequalities are easy to understand from equation (F.4); they come from conver- 
gence of the integral when £1 — > £ 2 and when either £1 or £ 2 or both go to infinity. The 
inequalities Re(0i) > 0, Re(0 2 ) > that we assumed in deriving (F.9) automatically 
follow from these, but the third inequality, Re(0i + 2 — 03) < 1, does not. This final 
inequality is somewhat mysterious becuse it breaks the symmetry between 01,02,03 
and does not follow in any obvious way from the convergence of (F.4). In fact it is 
a relic of our method of evaluating the integral; in deriving (F.10) we used (F.7) to 
introduce a factor of t^- — r which then cancelled out of the integral in the end. 

r(l+<T3— 01— <X2) to 

We could have avoided this inequality by deforming the contour in (F.7). For another 
way to see that the inequality is spurious, we observe that if we had used the symmetry 
to send £1 or £ 2 to infinity instead of £ 3 in deriving (F.4) then a different inequality 
related to this one by symmetry would have appeared that we did not need to use in 
our evaluation. Thus the only conditions for convergence of the integral (F.l) are given 
by (F.20). 
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